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Abstract 

This thesis is based on three different projects, all of them are directly linked to 
the classical general theory of relativity, but they might have consequences for 
quantum gravity as well. 

The first chapter deals with pseudo-Finsler geometric extensions of the classi- 
cal theory, these being ways of naturally representing high-energy Lorentz sym- 
metry violations. In this chapter we prove a certain type of "no-go" result for 
significant number of theories. This seems to have important consequences for 
the question of whether some weaker formulation of Einstein's equivalence prin- 
ciple is sustainable, if (at least) certain types of Lorentz violations occur. 

The second chapter deals with the problem of highly damped quasi-normal 
modes related to different types of black hole spacetimes. First, we apply to this 
problem the technique of approximation by analytically solvable potentials. We 
use the Schwarzschild black hole as a consistency check for our method and de- 
rive many new and interesting results for the Schwarzschild-de Sitter (S-dS) black 
hole. One of the most important results is the equivalence between having a ratio- 
nal ratio of horizon surface gravities and periodicity of quasi-normal modes. By 
analysing the complementary set of analytic results derived by the use of mon- 
odromy techniques we prove that all our theorems almost completely generalize 
to all the known analytic results. This relates to all the types of black holes for 
which quasi-normal mode results are currently known. 

The third chapter is related to the topic of multiplication of tensorial distri- 
butions. We focus on an alternative approach to the ones presently known. The 
new approach is fully based on the Colombeau equivalence relation, but techni- 
cally avoids the Colombeau algebra construction. The advantage of this approach 
is that it naturally generalizes the covariant derivative operator into the general- 



ized tensor algebra. It also operates with much more general concept of piecewise 
smooth manifold, which is in our opinion natural to the language of distributions. 
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Introduction 



This thesis is based on three different projects. All of these three projects lie some- 
where on the boundary between the classical general theory of relativity and the 
so far unknown quantum theory of gravity All of them deal with the phenom- 
ena that can be fully understood by classical (non-quantum) language, but at the 
same time describe, or at least strongly indicate, high energy modifications of the 
classical theory, mostly due to quantum gravity 

Ever since Einstein's general theory of relativity was established as a highly 
successful theory of gravity, physicists have also realized its unavoidable limits. 
The first set of problems is due to the theory itself. As is proven by rigorous math- 
ematical theorems, physically reasonable situations lead in general relativity to 
solutions containing singularities. The other set of problems is due to the fact that 
the theory is classical. To be able to consistently describe the full interaction be- 
tween quantum fields and gravity, one naturally needs to go beyond the classical 
language and somehow quantize gravity. These two sets of problems are, how- 
ever, closely related. Singularities appear in the situations where we do not expect 
the classical description to be relevant. 

There are presently many ideas, suggestions and conjectures about the quan- 
tum theory of gravity. They are related to the many different approaches to the 
problem that have appeared over the last decades, some of them having quite 
different backgrounds. (These are approaches such as string theory, loop quan- 
tum gravity, non-commutative geometries, etc.) Despite the fact that the basic ap- 
proaches are very different, there is a wide agreement about the results obtained 
within the semi-classical regime. (For example, the effect of Hawking radiation.) 



Surprisingly, the results of semi-classical approach already seem to give us con- 
siderable information about some of the objects of full quantum gravity, such as 
black holes. These results (such as the relation between horizon area and black 
hole entropy) must be recovered by every quantum theory of gravity, that has 
ambitions to be correct. 

First two chapters in my thesis present results obtained in close collaboration 
with my supervisor Prof. Matt Visser. They were published in various journals 
(see the list of publications at the end of the thesis). The last chapter is related to 
my own work and its significantly shortened version will be prepared for journal 
submission in the immediate future. 

The first chapter of the thesis deals with some possible classical geometric ex- 
tensions of the general theory of relativity. Such geometric extensions are related 
to possible high-energy Lorentz symmetry violations, by many physicists being 
assumed to be one of the possible effects of the future quantum gravity theory. To 
search for a link between some generalized geometry and possible high-energy 
Lorentz violations seems to be, with respect to the spirit of Einstein's equivalence 
principle, a very natural approach. 

The second chapter deals with a semi-analytic approach to the highly damp- 
ed quasi-normal modes of various classical black hole space-times. It is in fact 
a topic from classical general relativity, but it is also of particular interest of the 
quantum gravity community. This is because the asymptotic quasi-normal modes 
behaviour is suspected |77llll0| to be connected to the area spectrum of the quan- 
tum black holes. 

The third chapter represents a topic from the field of mathematical physics 
closely related to the general theory of relativity. It offers new ideas about how 
to fully generalize the language of differential geometry into the distributional 
framework. It is again a conceptual extension which has direct relevance for the 
classical theory, but it might not be unreasonable to assume that it can have im- 
portant consequences for quantizing gravity as well. 

At the end of the thesis we summarize our results. This is followed by ap- 
pendices containing some of the most common physics and mathematics results 



relevant for the calculations in this thesis. 



Chapter 1 

Pseudo-Finsler extensions to gravity 



1.1 Introduction 

Possible low-energy manifold-like limits of quantum gravity One of the most 
significant problems in the recent history of theoretical physics is the problem 
of quantization of gravity. Most of the candidates for a quantum gravity theory 
suggest that the description of spacetime is at the fundamental level far from the 
traditional concept of manifold. 

Despite this fact, one might be still interested in whether these theories have a 
deeper manifold-like low-energy limit; one more subtle than the ordinary pseudo- 
Riemannian geometry. If they have such limit, it must be definitely an extension 
of pseudo-Riemannian geometry, but in the same time it must be "close to" the 
highly successful concept of pseudo-Riemannian geometry. 

Ultra-high energy violations of Lorentz symmetry The theoretical physics com- 
munity has recently exhibited increasing interest in the possibility of ultra-high- 
energy violations of Lorentz invariance (^HS5HS5HS7HIQ8HIMII11HIS|. Specif- 
ically, recent speculations regarding Lorentz symmetry breaking and /or funda- 
mental anisotropics and /or multi-refringence arise separately in the many and 
various approaches to quantum gravity. 
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Such phenomena arise in loop quantum gravity 11251 1591 , string models H96l 
1112111131 , and causal dynamical triangulations HI 12, and are also part and parcel 
of the "analogue spacetime" programme |15|, and of many attempts at devel- 
oping "emergent gravity" l|T6ll24ll97| . Recently, the ultra-high energy breaking 
of Lorentz invariance has been central to the Horava-Lifshitz models ||78l |79l |80l 
11541 11551 11751 - Of course not all models of quantum gravity lead to high-energy 
Lorentz symmetry breaking, and the comments below should be viewed as ex- 
ploring one particular class of interesting models. 

The connection between Lorentz violations and geometry The extensions of 
pseudo-Riemannian geometry can typically modify dispersion relations, so one 
can be interested in seeing if such modified dispersions relations can be natu- 
rally embedded in some extension of pseudo-Riemannian geometry These exten- 
sions could be then naturally viewed as a low energy manifold-like non-pseudo- 
Riemannian limit of a given quantum gravity approach. 

In the other direction, if we wish to follow the spirit of Einstein's equivalence 
principle and develop a geometric spacetime framework for representing Lorentz 
symmetry breaking, either due to spacetime anisotropics or multi-refringence, 
then it certainly cannot be standard pseudo-Riemannian geometry. 

This strongly suggests that carefully thought out extensions and modifications 
of pseudo-Riemannian geometry might be of real interest to both the general rel- 
ativity and high-energy communities. 

Why focus on the light cone structure? In particular, when attempting to gen- 
eralize pseudo-Riemannian geometry, the interplay between the "signal cones" of 
a multi-refringent theory and the generalized spacetime geometry is an issue of 
considerable interest: 

• In multi-refringent situations it is quite easy to unify all the signal cones in 
one single Fresnel equation that simultaneously describes all polarization 
modes on an equal footing. 
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• In a standard manifold setting, where we retain the usual commutative co- 
ordinates, we shall see that it is natural to demand that each polarization 
mode can be assigned a specific geometric object. This object is in fact a 
Lorentzian analogue of what mathematicians know as Finsler norm (see the 
next section for details). 

• In standard general relativity the (single, unique) signal cone almost com- 
pletely specifies the spacetime geometry — one needs only supplement the 
signal cone structure with one extra degree of freedom at each point in 
spacetime, an overall conformal factor, in order to completely specify the 
spacetime metric, and thereby completely specify the geometry. This is ul- 
timately due to the fact that in standard pseudo-Riemannian geometry the 
scalar product is a simple bi-linear operation. Unfortunately in the more 
general pseudo-Finsler geometry^! life is more difficult, but one might still 
have a hope that following the guideline of simplicity the light-cone behav- 
ior could hold a crucial piece of information about the overall geometry. 

Bi-refringent crystal and beyond Considerable insight into such Finsler-like 
models can be provided by considering the "analogue spacetime" programme, 
where analogue models of curved spacetime emerge at some level from well un- 
derstood physical systems IIT5I . In particular, the physics of bi-axial bi-refringent 
crystals [29] provides a particularly simple physical analogue model for the math- 
ematical object introduced some 155 years ago by Bernhard Riemann |132[ , and 
now known as Finsler distances (again see the next section for more details). We 
shall soon see that this mathematical object can reasonably easily be extended to 
a Lorentzian signature pseudo-Finsler spacetime, with an appropriate pseudo- 
Finsler norm. 

Note that we are not particularly interested in the properties of bi-refringent 



1 For the details of what we mean by "pseudo-Finsler geometry" see the next section. 

2 It must be emphasized that, despite many misapprehensions to the contrary uni-axial bire- 
fringent crystals are relatively uninteresting in this regard; they do not lead to Finsler 3-spaces, 
but "merely" yield bi-metric Riemannian 3-geometries. 
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crystals per se, we use them only as an exemplar of Finsler 3-space and Finsler 
space-time, as a guidepost to more complicated things that may happen in Fins- 
lerian extensions to general relativity As a result of this fact in the next step we 
show that all our observations from the particular bi-refringent crystal case hold 
in general bi-metric situationsj. This is because our ultimate goal is to be able 
to say something about the (presumed) low-energy manifold-like limit of what- 
ever quantum theory (or class of quantum theories) is leading to a bi-metric / 
bi-refringent theory approximately reproducing Einstein gravity. 

The no-go result While the use of Finsler 3-spaces to describe crystal optics is 
reasonably common knowledge within the community of mathematicians and 
physicists studying Finsler spaces, it is very difficult to get a clear and concise 
explanation of exactly what is going on when one generalizes to Lorentzian sig- 
nature space-time. In particular the fact that any relativistic formulation of Finsler 
space needs to work in Lorentzian signature (- +++), instead of the Euclidean 
signature (+ +++) more typically used by the mathematical community, leads to 
many technical subtleties (and can sometimes completely invalidate naive con- 
clusions). Moreover unlike the spacetime Finsler norm, defining a spacetime Finsler 
metric is fraught with technical problems. These problems seem to be fundamen- 
tal and hold in arbitrary bi-metric situations. 
So the basic things we assert are: 

• What is exceedingly difficult, and we shall argue is in fact outright impossi- 
ble within this framework, is to construct a unified and still simple formal- 
ism that moves "off-shell" (off the signal cones). 

• This is a negative result, a "no-go theorem", which we hope will focus at- 
tention on what can and cannot be accomplished in any natural way when 
dealing with multi-refringent anisotropic Finsler-like extensions to general 
relativistic spacetime. 



3 For the difference between bi-refringence and bi-metricity see for example 111761 . 
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To this end, our "no-go" result indicates that the popular assumption that 
anisotropics and multi-refingence are likely to occur in "quantum gravity" leads 
to significant difficulties for the Einstein equivalence principle — since even the 
loosest interpretation of the Einstein equivalence principle would imply the neces- 
sity of a coherent formalism for dealing with all signal comes, and the spacetime 
geometry, in some unified manner. We conclude that, despite the fact that space- 
time anisotropics and multi-refringence are very popularly assumed to be natural 
features of "quantum gravity", and while these features have a straightforward 
"on-shell" implementation in terms of a suitably defined Fresnel equation, there 
is no natural way of extending them "off-shell" and embedding them into a single 
over-arching spacetime geometry 

But to remind the reader: if one steps outside of the usual manifold picture, 
either by adopting non-commutative coordinates, or even more abstract choices 
such as spin foams, causal dynamical triangulations, or string-inspired models, 
then the issues addressed in this chapter are moot — our considerations are rele- 
vant only insofar as one is interested in the first nontrivial deviations from exact 
low-energy Lorentz invariance, and only relevant insofar as these first nontrivial 
effects can be placed in a Finsler-like setting. 

The structure of this chapter This chapter begins with introducing the concept 
of Finsler geometry. This is followed by proving our "no-go" result for the partic- 
ular example of the bi-refringent crystal analogue model. After this we show that 
the result holds in arbitrary bi-metric situation. At the end of this chapter we add 
a section where we explore the general conditions, which any pseudo-Finslerian 
geometry must fulfil in order to give bi-refringence. But as a consequence of our 
"no-go" result, such constructs represent only a "complicated" and non-intuitive 
route for how to recover bi-refringence by pseudo-Finslerian geometry. 
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1.2 Basics of (pseudo-)Finsler geometry 

Mathematically, we define a Finsler function (Finsler norm, Finsler distance func- 
tion) ||2~T[|58ll to be a C-valued function F(x, v) on the tangent bundle to a manifold, 
such that it is homogeneous of degree 1: 

F(x,kv) = kF(x,v), k>0, x e M, veT x M. (1.1) 

This then allows one to define a notion of distance on the manifold, as the minimal 
value of the functional 



S(x{ti),x{tf)) 



*u d -t) 



dt, (1.2) 



which is now guaranteed to be independent of the specific parameterization t. 

By "pseudo-Finsler geometry" we mean Finsler geometry with Lorentzian sig- 
nature. Now by Lorenzian signature of the general Finsler metric (see the equation 
11.4)) we mean, that for any arbitrary vector taken as an argument of the metric we 



obtain matrix with ( — h ++) signature. A pseudo-Riemannian norm is only a 
special case of a pseudo-Finsler norm. For a pseudo-Riemannian manifold with 
metric g a b{x) one would take 

F(x,v) = ^g ab (x)v a v\ (1.3) 

but for a general pseudo-Finslerian manifold the function F(x, v) is arbitrary ex- 
cept for the 1-homogeneity constraint in v and the metric signature constraint. 
Note that in Euclidean signature (where g a b(x, v) is taken to be for any arbitrary v 
a positive definite matrix), the general Finsler function F(x, v) is typically smooth 
except at v = 0. In Lorentzian signature however, F(x, v) is typically non-smooth 
for all null vectors — so that non-smoothness issues have grown to affect (and in- 
fect) the entire null cone (signal cone). As we shall subsequently see below, some- 
times a suitable higher algebraic power, F 2n (x,v), of the pseudo-Finsler norm is 
smooth. 

To ensure smoothness of the (pseudo-)Finsler metric, defined below, it is enou- 
gh to weaken the condition that F(x, v) shall be smooth and to demand only that 
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the square F 2 (x, v) be smooth, except possibly at v = 0. It is standard to define the 
(pseudo-)Finsler metric as 

, , 1 d 2 [F 2 {x,v)} 

which then satisfies the constraint that it is homogeneous of order zero 

g ab (x, kv) = g ab (x,v), k>0. (1.5) 

This can be viewed as a "direction-dependent metric", and is clearly a significant 
generalization of the usual (pseudo-)Riemannian case. One can immediately see 
that pseudo-Riemannian metric fulfills this definition with respect to the pseudo- 
Riemannian norm (|1.3|) . 

Almost all of the relevant mathematical literature has been developed for the 
Euclidean signature case. Because of this assumption, any mathematical result 
that depends critically on the assumed positive definite nature of the matrix of 
metric coefficients cannot be carried over into the physically interesting pseudo- 
Finsler regime, at least not without an independent proof that avoids the posi- 
tive definite assumption. (Unfortunately it is not uncommon to find significant 
mathematical errors in the pseudo-Finsler physics literature due to neglect of this 
elementary point.) Basic references within the mathematical literature include 
EQ|58]|. 

The Legendre transformation between a vector tangent space at the point x 
and its dual: V x — > V*, is defined as 

h(v) = gab(x,v)v a (1.6) 

Then the dual (pseudo-)Finsler norm F* can be defined by the condition: 

F»(l(v)) = F(v). (1.7) 

The dual metric is again naturally obtained as: 

^,, )= i*&a. ( ,8) 

2 dv a dv b 
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All this is a natural generalization from the (pseudo-)Riemannian case. The 
construction of a full (pseudo-)Finsler geometry is in general significantly more 
complicated than in the (pseudo-)Riemannian subcase. But since the definition 
of objects like non-linear connection, Finsler connection, (etc.), is not needed for 
the purpose of this chapter, it will be omitted here and left for the specialized 
literature (see for example [|133 | ] ). 

1.3 Analogue model: Birefringent crystal 

1.3.1 Outline 

Purely for the purposes of developing a useful analogy, which we shall use as a 
guide to the mathematics we wish to develop, we will focus on the optical physics 
of bi-axial bi-refringent crystals. After the basic definitions are presented, we will 
show how various purely spatial 3-space Finsler structures arise. (Many purely 
technical details, when not directly involved in the logic flow, will be relegated to 
the appendix [A]) We again emphasize that uni-axial bi-refringent crystals, which 
are what much of the technical literature and textbook presentations typically fo- 
cus on, are for our purposes rather uninteresting — uni-axial bi-refringent crystals 
"merely" lead to bi-metric Riemannian space-times and are from a Finslerian per- 
spective "trivial". Such crystals are only one particular example demonstrating 
general difficulties with Finsler representation of bi-metric theories. The gener- 
alization from this example to any bi-metric case is presented in the following 
section. 

We shall soon see that even in three-dimensional space there are at least four 
logically distinct Finsler structures of interest: On the tangent space each of the 
two photon polarizations leads, via study of the group velocity, to two quite dis- 
tinct Finsler spacetimes. On the co-tangent space each of the two photon polariza- 
tions leads, via study of the phase velocity, to two quite distinct co-Finsler space- 
times. The inter-relations between these four structures is considerably more sub- 
tle than one might naively expect. 
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Additionally, (apart from some purely technical difficulties along the optical 
axes in bi-axial crystals), each of these four 3-dimensional spatial Finsler struc- 
tures has a natural 4-dimensional extension to a spacetime pseudo-Finsler struc- 
ture. Beyond that, there are reasonably natural ways of merging the two photon 
polarizations into "unified" Finsler and co-Finsler norms, closely related to the 
appropriate Fresnel equation, though the associated Finsler metrics are consider- 
ably more problematic — all these mathematical constructions do come with a 
price — and we shall be careful to point out exactly where the technical difficul- 
ties lie. Finally, using this well-understood physical system as a template, we shall 
(in the spirit of analogue spacetime programme) then ask what this might tell us 
about possible Finslerian extensions to general relativity, and in particular to the 
subtle relationship between bi-refringence and bi-metricity, (or more generally, 
multi-refringence and multi-metricity). 

Specifically, we have investigated the possibility of whether one can usefully 
and cleanly deal with both Finsler structure (anisotropy) and multi-refringence 
simultaneously. That is, given two (or more) "signal cones": Is it possible to natu- 
rally and intuitively construct a "unified" pseudo-Finsler spacetime such that the 
pseudo-Finsler metric specifies null vectors on these "signal cones", but has no 
other zeros or singularities? Our results are much less encouraging than we had 
originally hoped, and lead to a "no-go" result. 

1.3.2 Space versus space-time: Interpretations of the Finsler and 
co-Finsler structures 

The key physics point in bi-axial bi-refringent crystal optics is that the group ve- 
locities, and the phase velocities, are both anisotropic and depend on direction in 
a rather complicated way [29 j. Technical details that would detract from the flow 
of the text are relegated to appendix lAl 
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From group velocity to pseudo-Finsler norms 

We can summarize the situation by pointing out that the group velocity is given 

by 

v 2^ = 92 (n, n) ± y/g 2 (n, n) 2 - g (n, n) (n ■ n) 
9 n 9b(n,n) 

where ^(n, n) and qo(n, n) are known quadratic functions of the direction n and 
are given as 

~ ( \ 2-2-2,2-2-2,2-2-2 /1 1 r\\ 

q Q {n,n) = n x v y v z +n y v x v z + n z v x v y ; (1.10) 

q 2 {n, n) = \{nl{v y 2 + vj 2 ) + nj(vj 2 + v~ 2 ) + n\{v~ 2 + v~ 2 )). (1.11) 

The coefficients in these quadratic forms are explicit functions of the compo- 
nents of the 3x3 permittivity tensor. (See appendix lAl) The function f g (n) so 
defined is homogeneous of degree zero in the components of n: 

Vg{KXl) = v g (n) = v g {h). (1.12) 

The homogeneous degree zero property should remind one of the relevant feature 
exhibited by the Finsler metric. There is a natural connection between the concept 
of group velocity and the geometric objects on a tangent space (rather than a co- 
tangent space). This is given by the fact that group velocity describes how energy 
propagates. 

Let us now first define the quantities 



v g (n) 



^3±(n) = -— — = A/ g 2 (n, n) =f \fth(n, n) 2 - q (n, n) (n ■ n), (1.13) 



or adopt the perhaps more transparent notation 

Idxl 



^3±(dx) = = \Jq 2 (dx, dx) =f y/q 2 (dx, dx) 2 - g (dx, dx) (dx • dx) . (1.14) 



The quantity l/v fl (n) appearing in (|1.14[) is in the literature often called "slow 



ness". (|1.14[) is by inspection a 3-dimensional (Riemannian) Finsler distance de- 
fined on space, having all the correct homogeneity properties, F 3 ±(ndx) = \k\ F 3 ±(dx) 
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Physically, the Finsler distance is in this situation the time taken for the wavepacket 
to travel a distance dx. 

To now extend the construction given above to full (3+1) dimensional space- 
time, we first define a generic 4-vector 



dX = (di;dx), 
and then formally construct 

F 4± (dX) = v/-(dt)2 + F 3± (dx)2. 
That is 

F 4± (dx) = J-(dty 

Even more explicitly, one may write 

F 4± (dX) 



dx ■ dx 

v g (dx) 2 ' 



(1.15) 



(1.16) 



(1.17) 



- (d£) 2 + g 2 (dx, dx) =p A/g 2 (dx, dx) 2 - g (dx, dx) (dx • dx) 
The null cones (signal cones) of F 4 ±(dX) are defined by 

F 4± (dX) = <£> ||dx|| = v g {dx) dt. 



1/2 



(1.18) 



(1.19) 



So far this has given us a very natural pair of (3+l)-dimensional pseudo-Finsler 
structures in terms of the ray velocities corresponding to the two photon polariza- 
tions. 



For future use, let us now formally define the quantity 

ds 4 = {F 4 (dX)} 4 

= {F 4+ (dX)F 4 _(dX)} 2 

= (dt) 4 - 2(dt) 2 g 2 (dx, dx) + g (dx, dx) (dx • dx). 



(1.20) 



This certainly provides an example of a specific and simple 4 <h -root Finsler norm 
that can naturally and symmetrically be constructed from the two polarization 



CHAPTER 1. PSEUDO-FINSLER EXTENSIONS TO GRAVITY 15 

modes, and its properties (and defects) are certainly worth investigating. Phys- 
ically the condition ds = defines a double-sheeted conoid (a double-sheeted 
topological cone) that is the union of the propagation cone of the individual pho- 
ton polarizations. This Finsler norm defines the Finsler geometry naturally uni- 
fying the two original geometries. It is also very close to quartic extension of the 
notion of distance that Bernhard Riemann speculated about in his inaugural lec- 
ture (see El). 

However we shall soon see that when it comes to defining a Finsler space- 
time metric this construction nevertheless leads to a number of severe technical 
difficulties; difficulties that can be tracked back to the fact that we are working in 
non-Euclidean signature. 

From phase velocity to pseudo-co-Finsler norms 

In counterpoint, as a function of wave-vector the phase velocity is 



*-, _ g 2 (k,k)± v /g 2 (k,k)2-g (k,k)(k-k) 
where the quadratics (72 (k, k) and go(k, k) are now given by equations 

g„(k, k) = k 2 xV y z + k 2 yV y z + ky x v 2 y , (1.22) 

and 

q 2 (k, k) = \ (klivl + vl) + k 2 y (vl + vl) + k 2 M + v 2 )) . (1.23) 

This expression is homogeneous of order zero in k, so that 

v p (Kk)=v p (k) = v p (k). (1.24) 

Again, we begin to see a hint of Finsler structure emerging. Because k is a wave- 
vector it transforms in the same way as the gradient of the phase; thus k is most 
naturally thought of as living in the 3-dimensional space of co-tangents to physical 
3-space. Let us now define a co-Finsler structure on that co-tangent space by 



G 3± (k) = v p (k)||k|| = Vg2(k,k)±v/g2(k,k) 2 -go(k,k)(k-k). (1.25) 
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We use the symbol G rather than F to emphasize that this is a co-Finsler structure, 
and note that this object satisfies the required homogeneity property 



G 3 ±(Kk) = |«|G 3 ±(k). 



(1.26) 



Now let us go for a (3+1) dimensional spacetime interpretation: Consider the 4 
co-vector 

k = (w; k) , 

and define 



That is 

More explicitly 

G 4± (k) 



G 4± (k) = ^-^ + G 3± (ky. 

G 4± (k) = 



-oo" 



Up (k)2(k.k). 



u 2 + q 2 (k, k) ± v /g 2 (k,k) 2 -g (k,k)(k-k) 



1/2 



(1.27) 
(1.28) 
(1.29) 

(1.30) 



We again see that this object satisfies the required homogeneity property 

G 4 ±(Kk) = \K\G 4 ±(k), (1.31) 

so that this object is indeed suitable for interpretation as a co-Finsler structure. 
Furthermore the null co-vectors of G 4 are defined by 



G 4 ±{k) = <^> u — v p (k) ||k||, 



(1.32) 



which is exactly the notion of dispersion relation for allowed "on mass shell" wave- 
4-vectors that we are trying to capture. Thus G 4 ± lives naturally on the co-tangent 
space to physical spacetime, and we can interpret it as a pseudo-co-Finsler struc- 
ture. 

We can again define a "unified" quantity 



G 4 (kf = {G 4+ {k) G 4 -{k)} 2 

= c/-2u; 2 g 2 (k,k) + g (k,k)(k-k). 



(1.33) 
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Physically, the condition G±(k) = simultaneously encodes both dispersion re- 
lations for the two photon polarizations. It defines a double-sheeted conoid (a 
double-sheeted topological cone) that is the union of the dispersion relations of 
the individual photon polarizations. The vanishing of G±(k) can be viewed as a 
Fresnel equation, and can indeed be directly related to Fresnel's condition for the 
propagation of a mode of 4-wavenumber k = (to; k). As is the case for F 4 (dX), 
we shall soon see that this construction (once one tries to extract a spacetime co- 
Finsler metric) nevertheless leads to a number of severe technical difficulties; dif- 
ficulties that can again be tracked back to the fact that we are now working in 
non-Euclidean signature. 

1.3.3 Technical issues and problems 

The situation as presented so far looks very pleasant and completely under control 
— and if what we had seen so far were all there was to the matter, then the study 
of pseudo-Finsler space-times would be very straightforward indeed — but now 
let us indicate where potential problems are hiding. 

• Note that up to this stage we have not established any direct connection 
between the Finsler functions F 3 ±(n) and the co-Finsler functions Cr3±(k). 
Physically it is clear that they must be very closely related, but (as we shall 
soon see) establishing the precise connection is tricky. 

• Furthermore, the transition from Finsler distance to Finsler metric requires at 
least two derivatives. Even in Euclidean signature this places some smooth- 
ness constraints on the Finsler distance, smoothness constraints that are non- 
trivial and not always satisfied. 

• Especially, there are problematic technical issues involving the 4-dimension- 
al spacetime Finsler and co-Finsler metrics — certain components of the met- 
ric are infinite, and this time the potential pathology is widespread. (In 
Lorentzian-like signature situations potential problems tend to infect the en- 
tire null cone.) 
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The Finsler and co-Finsler 3-metrics 

The standard definition used to generate a Finsler metric from a Finsler distance is 

to set: 

l d 2 [F 3± (n) 2 ] 

^' (n) = 2 dnW ' (134) 

which in this particular case implies 

1 <9 2 [g 2 (n, n) =f \fth(n, n) 2 - g (n, n) (n • n)] 

^' (n) = 2 d^ • (L35) 

It is convenient to rewrite the quadratics as 

g 2 (n, n) = [q 2 ]ij n l n 3 ] (1.36) 

gb(n, n) = [g ]y n* n J ; (1.37) 
since then we see 

<7y(n) = [g 2 ]ij =F (discriminant contributions). (1.38) 

Unfortunately we shall soon see that the contributions coming from the discrimi- 
nant are both messy, and in certain directions, ill-defined. This is obvious from the 
fact that squares of both Finsler functions are not even everywhere differentiable. 
Similarly we can construct a Finsler co-metric: 

which specializes to 



hij(k) = 1 d"[q 2 (k, k) T ygg k) 2 - g (k, kHl^k)] 
2 dkidkj 

It is again convenient to rewrite the quadratics as 

q 2 (k,k) = [q 2 ] ij k i k j ; (1.41) 

g (k,k) = [g ]^A;,; (1.42) 

since then we see 

/i y (k) = [g 2 ] tJ =F (discriminant contributions). (1-43) 

Again we shall soon see that the contributions coming from the discriminant are, 
in certain directions, problematic. 
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Technical problems with the Finsler 3-metric 

Consider the (ray) discriminant 

D = q 2 {n,n) 2 -q {n,n) (n-n). (1.44) 

There are three cases of immediate (mathematical) interest: 

Isotropic: If v x = v y = v z then D = 0; in this case the two Finsler functions 
F± are equal to ech other. F 3 (dx) then describes an ordinary Riemannian geome- 
try, and F 4 (dX) an ordinary pseudo-Riemannian geometry. This is the standard 
situation, and is for our current purposes physically uninteresting. 

Uni-axial: If one of the principal velocities is distinct from the other two, then 
we can without loss of generality set v x — v y — v and v z — v e . The discriminant 
then factorizes into a perfect square 

(v 2 - vpjn l + n 2 v ) 
2v 2 v; 



J o e 



D = { V ° :; 2 J yJ \ . (1-45) 



In this case it is immediately clear that both F 3± (dx) reduce to simple quadratics, 
and so describe two ordinary Riemannian geometries. Indeed 



F 3+ (n) = ^; (1.46) 



and 



v'l 



,2j_ r ,2 ^2 



^(n) = ^±^ + ^f. (1.47) 

vi vi 



In the language of crystal optics v and v e are the "ordinary" and "extraordinary" 
ray velocities of a uni-axial birefringent crystal. In geometrical language the two 
photon polarizations "see" distinct Riemannian 3-geometries F 3± (dx) and distinct 
pseudo-Riemannian 4-geometries F 4± (dX) — this situation is referred to as "bi- 
metric". This situation is for our current purposes physically uninteresting. 

Bi-axial: The full power of the Finsler approach is only needed for the bi-axial 
situation where the three principal velocities are distinct. This is the only situation 
of real physical interest for us, as it is the only situation that leads to a non-trivial 
Finsler metric. In this case we can without loss of generality orient the axes so that 
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v x > v y > v z . There are now two distinct directions in the x-z plane where the 
discriminant vanishes — these are the called the (ray) optical axes. After some ma- 
nipulations that we relegate to Appendix I A.51 the discriminant can be factorized 
as 



D 



ivivi 



X 



[(n ■ n) - (e x ■ n) 2 ] [(n ■ n) - (e 2 ■ n) 2 ] 



(1.48) 



where the two distinct (ray) optical axes are 



ei, 2 



V r \ V 1 — V 1 ' 



:-* 



y y 



v 2 — v 2 



r \l v 2 x - v 2 y 



(1.49) 



Note that ei 2 are unit vectors (in the ordinary Euclidean norm) so that the dis- 
criminant D vanishes for any n oc e lj2 , and does not vanish anywhere else. We 
can thus introduce projection operators Pi and P 2 and write 



Pi(n, n) = (n ■ n) - (e 1 • n) 2 ; 
P 9 .(n, n) = (n ■ n) - (e 2 • n) 2 . 



21",- 



Combining this with our previous results: 



{^3 ± (n)r 



2v 2 x v 2 z 



g 2 (n, n) T ^ .^ JP^n, n) P 2 (n, n) 



(1.50) 
(1.51) 

(1.52) 



If we now calculate the Finsler metric [<7 3 ±(n)]y we shall rapidly encounter tech- 
nical difficulties due to the discriminant term. To make this a little clearer, let us 
define 

d 2 Jp 1 (n,n)P 2 (n,n) 



[P*(*)]i 



dn l dn? 



since then 



[g 3 (n)] tJ = [fc(n)fo =F %^[A(n)]y. 

Zi ~. u „ 



(1.53) 
(1.54) 



Temporarily suppressing the argument n, we have 



[P 



3\ij 



2 



Pi - /Pi 

■' i/ Pi J V Pi 



(1.55) 
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A brief computation now yields the rather formidable result 



[P 



3\ij 



n V "2 



1 

+ 4 



9iA &p 2 + aP 2 ap 



p 



1/2 



'PlP 2 



J ' diPi djP! ^ - diP 2 djP 2 



pi 



1/2 



^ 



3/2 



P 



3/2 



1 



2 V / ^P 2 
1 



= [ didjPi P 2 + 9^P 2 P 1 ] 



+ 



4VPiP2 



aA «9,p 2 + aA a,-A - aA aPi § - aA ^p 2 § 

PL P2 



. (1.56) 



From this expression it is clear that along either optical axis, (as long as the op- 
tical axes are distinct, which is automatic in any bi-axial situation), some of the 
components of [^3]^, and therefore some of the components of the Finsler metric 
[93±]ij = [q2\ij ± (constant) x [P 3 ] ijr wi\lbe infinite. 

To see this in an invariant way, let u and w be two 3-vectors and consider 



[P 3 (n)](u,w) = [P 3 ] ij u i w i . 



(1.57) 



After a brief computation: 

[P 3 (n)](u,w) 



2 v /P 1 (n,n)P 2 (n ! n^ 

1 
+ 



[ Pi(ll, w) P 2 (n, n) + P 2 (u, w) P 1 (n, n) ] 



4 v /P 1 (n,n)P 2 (n,n 

-A(n,u)A(n,w) P?( "'" 



Pi(n, u) P 2 (n, w) + P 2 (n, u) Pi(n, w 

Pi(n, n 



P 2 (n, u) P 2 (n, w) 



. (1.58) 



Pi(n, n) "" "' P 2 (n, n 

This quantity will tend to infinity as n tends to either optical axis provided: 

• The optical axes are distinct. 

(If the optical axes are coincident then P 1 = P 2 and so P 3 degenerates to 

[P 3 (n)](u, w) -> P 1 (u, w) = P 2 (u, w). (1.59) 
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One recovers the [for our purposes physically uninteresting] result for a uni- 
axial crystal.) 

• One is not considering the special case u = w = n. 
(In this particular special case P 3 degenerates to 



[P 3 (n)](n,n) -* ^Pi(n,n) P 2 (n,n), (1.60) 

which is well-behaved on either optical axis.) 
In summary: 

• The spatial Finsler 3-metric is [g^±\ij genetically ill-behaved on either optical axis. 

• This feature will also afflict the spacetime pseudo-Finsler 4-metric [g4±] a & de- 
fined by suitable derivatives of the Finsler 4-norm F±±. 

• This particular feature is annoying, but seems only to be a technical problem 
to do with the specifics of crystal optics, it does not seem to us to be a critical 
obstruction the developing a space-time version of Finsler geometry. Ulti- 
mately it arises from the fact that the "null conoid" is given by a quartic; this 
leads to two topological cones that for topological reasons always intersect, 
this intersection defining the optical axes. 

• It is the technical problems associated with the (3+1) spacetime Finsler metric, 
to be discussed below, which much more deeply concern us. 

Technical problems with the co-Finsler 3-metric 

The phase discriminant 

J D = g 2 (k,k) 2 -g (k,k)(k-k), (1.61) 

arising from the Fresnel equation (and considerations of the phase velocity) ex- 
hibits features similar to those arising for the ray discriminant. There are three 
cases: 
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Isotropic: If the crystal is isotropic, then D — 0. (This again is for our purposes 
physically uninteresting.) 

Uni-axial: If the crystal is uni-axial, then D is a perfect square 

..{««}' (1 .62) 

and so the co-Finsler structures G%± are both Riemannian: 

G 3+ (k)=v 2 k-n; (1.63) 

G 3 ^)=v 2 e (kl + kl)+v 2 k 2 z . (1.64) 

(This situation again is for our purposes physically uninteresting.) 

Bi-axial: Only in the bi-axial case are the co-Finsler structures G 3 ± "truly" Fins- 
lerian. There are now two distinct (phase) optical axes (wave-normal optical axes) 
along which the discriminant is zero, these optical axes being given by 




ei, 2 = ±i 75— £ ;i M-i , (1.65) 



in terms of which the phase discriminant also factorizes 

D = (t, «~ t, * )2 [(k • k) - (k ■ ei ) 2 ] [(k • k) - (k • e 2 ) 2 ] . (1.66) 

The co-Finsler norm is then (now using projection operators P\ and P 2 based on 
the phase optical axes ei i2 ) 

{G 3± (k)} 2 = g 2 (k, k) :p {VI ~ V " z) y/Pjk, k) P 2 (k, k). (1.67) 

The co-Finsler metric is defined in the usual way 

[ll3±m = i^|2#!l. (1 . 68) 

L J v ; 2 dk l dk3 v ' 

This now has the interesting "feature" that some of its components are infinite 
when evaluated on the (phase) optical axes. That is: The co-Finsler 3-metric is 
[^3±] u generically ill-behaved on either optical axis. This feature will also afflict 
the pseudo-co-Finsler 4-metric [/i4±] a & defined by suitable derivatives of the Finsler 
4-norm G±±. 
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Technical problems with the (3+1) spacetime interpretation 

The (3+l)-dimensional spacetime objects that give the best way how to merge two 
(3+1) pseudo-Finsler norms and pseudo-co-Finsler norms into one geometry are 
the quantities 

F 4 (dX) 4 = {F 4+ (dX) F 4 _(dX)} 2 ; (1.69) 

and 

G 4 {ky = {G 4+ (k) G 4 _(k)} 2 . (1.70) 

as defined in equation (|1.21[) and (|1.33|) . This is tantamount to taking 



F 4 (dX) = v / i 71 4 + (dX)F 4 _(dX); (1.71) 

and 



G A {k) = ^G i+ {k)G A .{k). (1.72) 

Now F 4 (dX) and G±(k) are by construction perfectly well behaved Finsler and 
co-Finsler norms, with the correct homogeneity properties — and with the nice 
and concise physical interpretation that the vanishing of F A (dX) defines a double- 
sheeted "signal cone" that includes both polarizations, while the vanishing of 
G±{k) defines a double-sheeted "dispersion relation" ("mass shell") that includes 
both polarizations. (Thus F 4 (dX) and G±{k) successfully unify the "on-shell" be- 
haviour of the signal cones in a Fresnel-like manner.) 

While this is not directly a "problem" as such, the norms F 4 (dX) and G±(k) do 
have the interesting "feature" that they pick up non-trivial complex phases: Since 
F 4 ±(dX) 2 is always real, (positive inside the propagation cone, negative outside), 
it follows that F±± (dX) is either pure real or pure imaginary. But then, thanks to 
the additional square root in defining F 4 (dX), one has: 

• F 4 (dX) is pure real inside both propagation cones. 

• F 4 (dX) is proportional to \/i = ^4^- between the two propagation cones. 

• F 4 (dX) is pure imaginary outside both propagation cones. 
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Similar comments apply to the co-Finsler norm G 4 (k). 

A considerably more problematic point is this: In the usual Euclidean signa- 
ture situation the Finsler norm is taken to be smooth everywhere except for the 
zero vector — this is usually phrased mathematically as "smooth on the slit tan- 
gent bundle". What we see here is that in a Lorentzian-like signature situation 
the Finsler norm cannot be smooth as one crosses the propagation cones — what 
was in Euclidean signature a feature that only arose at the zero vector of each tan- 
gent space has in Lorentzian-like signature situation grown to affect (and infect) 
all null vectors. The Finsler norm is here at best "smooth on the tangent bundle 
excluding the null cones". (In a mono-refringent case the squared norm, F 4 (cLY) 2 , 
is smooth across the propagation cones, but in the bi-refringent case one has to go 
to the fourth power of the norm, F 4 (cLY) 4 , to get a smooth function.) 

A "no go" result: Unfortunately, when attempting to bootstrap these two rea- 
sonably well-behaved norms to Finsler and co-Finsler metrics one encounters ad- 
ditional and more significant complications. We have already seen that there are 
problems with the spatial 3-metrics [<7 3 ±]y(n) and [/i 3± ] JJ (k) on the optical axes, 
problems which are inherited by the single-polarization spacetime (3+l)-metrics 
[fl'4±]ab(n) and [h 4 ±] ab (k), again on the optical axes. 

But now the spacetime (3+l)-metrics 

[9i]M = 2^^' (L73) 

and 

both have (at least some) infinite components — [gA\ a b{n) has infinities on the en- 
tire signal cone, and \h 4 ] ah {k) has infinities on the entire mass shell. Since the ar- 
gument is essentially the same for both cases, let us perform a single calculation: 

g ab = \ d a d h sJ[Fl F 2 ] (1.75) 



i* 



d b [F 2 + ] ?=- + d b [Fl] |t 



(1.76) 
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so that 

9ab 



1 

4 

+- 



d a d b [Fl] ?=■ + d b d b [Fl] |± 



d a F + d b F_ + d a F„d b F + - d a F + d b F+ 



d a F.d b F_- 



(1.77) 



That is, tidying up 

9ab = 



Pi P 2 



+ - 



d a F + d b F„ + d a F_d b F + - d a F + d b F+ 



d a F_d b F_ 



(1.78) 



' F + F_ 

The problem is that this "unified" metric g ab {n) has singularities on both of the 
signal cones. The (relatively) good news is that the quantity g ab (n) n a n b = F 2 {n), 
and so on either propagation cone F — > 0, so F(n) itself has a well defined limit. 
But now let n a be the vector the Finsler metric depends on, and let w a be some 
other vector. Then 

1 



a,„..b 



g ab (n) n a w 



n a w b d a d b [F 2 

\™ b d b [F 2 ] 

1 

2 



1 



-w 



d b [Fi 



d b [F 2 



+ 



J(b + ]^V)- + ([g.] ai r l V])-l. 



(1.79) 
(1.80) 
(1.81) 
(1.82) 

(1.83) 



2 [^- rju " ' F + w " J J/ F_ 

The problem now is this: g + and #_ have been carefully constructed to be individ- 
ually well defined and finite (except at worst on the optical axes). But now as we 
go to propagation cone "+" we have 



1 F_ 

g ab {n) n a w b ->■ -( [g+] ab n a w b ) — = oo, 

and as we go to the other propagation cone "— " we have 



g ab (n)n a w b ^^([g.] ab n a w b )^ 



oo. 



(1.84) 



(1.85) 
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So at least some components of this "unified" Finsler metric g a b( n ) are unavoid- 
ably singular on the propagation cones. Related (singular) phenomena have pre- 
viously been encountered in multi-component BECs, where multiple phonon mo- 
des can interact to produce Finslerian propagation cones [183J. 

Things are just as bad if we pick u and w to be two vectors distinct from "the 
direction we are looking in", n. In that situation 



a„.,b 



g ah {n) u a w 



g+{u,w) — + g-{u,w) — 

g+(u, n) g_(w, n) + g + (w, n) g_(u, n) 
F + F_ 

F- F + 

-9+ (u, n) g+(w,n)j^- g_ (u, n) g_ (w,n)^ 



(1.86) 



Again, despite the fact that both g + and g have been very carefully set up to be 
regular on the propagation cones (except for the known, isolated, and tractable 
problems on the optical axes), the "unified" metric g a b(n) is unavoidably singular 
there — unless, that is, you only choose to look in the nn direction. 

If we give up the condition of mathematical simplicity in the construction of 
the "unified" space-time (co-)Finsler norms, then F 4 and G4 do not have to fac- 
torize into a product of powers of the Finsler functions for individual polariza- 
tions. We could then look for more complicated ways of building "unified" Finsler 
and co-Finsler structures (constrained mainly by giving the correct light propaga- 
tion cones in the birefringent crystal), and we might be able to find an appropri- 
ate pseudo-Finsler geometry (which might also fulfill some additional reasonable 
physical conditions). The necessary conditions for such a geometry are relatively 
easy to formulate (and this is done in the last section), but because of the lack of 
any intuitive interpretation, and the corresponding lack of direct physical motiva- 
tion, the physical meaning of such an approach is highly doubtful 111451 . 
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1.4 General bi-metric situations 

In the previous part one could observe that the arguments (given the way they 
were constructed) might apply for a large set of different situations than the bi- 
refringent crystal. In fact they generally hold within the class of bi-metric theories. 
Take arbitrary bi-metric theory. Such theories contain two distinct pseudo- 
Riemannian metrics g^ h , so we can define two distinct "elementary" pseudo-Rie- 
mannian norms 

F±(x,v) = Jg±(x)v"v b . (1.87) 

Suppose one now wants a combined Finsler norm that simultaneously encodes 
both signal cones — then the natural thing is always to take 



F(x,v) = v / F + (x,v)F_(x,v); g abcd = g+ ah g cd) . (1.88) 

This construction for F(x, v) is automatically 1-homogeneous in v. Generally 
the vanishing of F(x, v) correctly encodes the two signal cones. So this definition 
of F(x, v) provides a perfectly good Finsler norm. 

However only from the fact that the individual F± are not positive definite 
automatically follows that: 

• F(x,v) is proportional to \fi — ^ between the two propagation cones, 
hence picks up a non-trivial phase. 

• As a result this Finsler norm is at best "smooth on the tangent bundle ex- 
cluding the null cones". 

For the "unified" metric, from the fundamental definition we see 

g a b(x,v) = 



9ab p "T" 9 a b p 



d a F + d h F + - d a F + d b F^ - d a F_d b F_?± 

r _|_ t — 



(1.89) 



This "unified" and "natural" Finsler metric g a b(x,v) has necessarily singularities 
on both of the signal cones. 
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The quantity g ab (x, v)v a v b = F 2 (x, v) is still being well defined, but if v a is the 
vector the Finsler metric depends on, and w a some other (arbitrary) vector, then 

g ab (x,v)v a w b 

= l -v a w b d a d h sjFlFl = l -w b d h ^f~K 



h> 



d b [Fl]^ + d b [Fl]^_ 



(9> a w b )^ + (g~ ab v a w b )^ 



(1.90) 



As we go to propagation cone "+", necessarily 



g ab (x,v)v a w b -> l -( g + b v a w b )^ = oo, (1.91) 

and as we go to the other propagation cone "-", necessarily 

g ab (x,v)v a w b -v ±(g- b v a w b )^- = oo. (1.92) 

So at least some components of this "unified" Finsler metric g ab (x,v) are (com- 
pletely generally) unavoidably singular on the propagation cones. And again, 
things are just as bad if we pick u and w to be two vectors distinct from v. 

In other words we see that all the problems which arose in the particular case 
of bi-refringent bi-axial crystal are in fact also features of arbitrary bi-metric situ- 
ation. 



1.5 General construction 



1.5.1 General constraints 



In the previous sections we observed that an "intuitive" way of encoding bi- 
metricity into pseudo-Finsler geometry leads to significant problems which seem 
inevitable and unavoidable. The other problematic part which was omitted is 
whether such intuitive construction leads to appropriate pseudo-Riemannian limit 
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for, (in some frame), slowly moving objects. Here we discuss all the general con- 
strains on a pseudo-Finsler geometry recovering arbitrary bi-refringence. 

What are the specific physical constraints given by bi-refringent theories on 
our geometry? We consider it meaningful to impose the following constraints: 

a) Locally there must exist a coordinate frame in which holds the following: 
At any arbitrary point from the domain of these coordinates take within this 
frame the purely time-oriented, (v, 0) vector. Then on some neighborhood 
of this vector the pseudo-Finsler norm approaches the pseudo-Riemannian 
norm. Also each "constant proper time" hypersurface, specified by F 2 (v ) = 
— Tq, is connected, and contains vectors (v, 0). 

b) The pseudo-norm must break Lorentz invariance as encoded in the "Fres- 
nel" equation giving the bi-refringence. 

Let us discuss these conditions a little bit: The first is just requirement that we 
want to recover, in some frame, low-energy physics as we know it (in a pseudo- 
Riemann form). It also means that if we have a massive particle, it should not 
be able to have two distinct sets of four-velocities between which it is not able 
to undergo a smooth transition. The last constraint it gives is that any massive 
particle should be able to move arbitrarily slowly with respect to the used "pre- 
ferred" frame. The second condition is trivial: it is the specific contribution of the 
given bi-refringent model (unlike the first condition, which can be considered as 
generic). 

1.5.2 How one proceeds for general bi-metric situations 

Co-metric structure construction 

Take the general bi-metric situation and define 

G(x,k) 4 = G + (x,k) 2 G_(x,k) 2 (1.93) 

where 

G±(x,k) 2 = gf{x)k a k b . (1.94) 
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This suggests a generic candidate for the pseudo-co-Finsler norm, which can be 
written as G(x,k) 2 oc G(x, k) 4 or more precisely G(x,k) 2 = M ± (x, k) G(x, k) 4 , 
where M\(x,k) is an otherwise arbitrary R-valued function fulfilling following 
constraints: 

x d[M 1 (x,k)G(x,k) 4 ] . , .. . , Tn Jr 

a ) ' is a bijection from V* -»• V; 

b) Mi(x, fc) is in A; a homogeneous map of degree —2; 

c) Mi (a;, k) is smooth on F*/{0}; 

d) Mi (x, k) is inside both signal cones negative, outside both positive, between 
them nonzero; 

e) Mi (x, k) G(x, k) 4 must approximate the pseudo-Riemann norm (squared) lo- 
cally in some coordinates for k close to (v, 0) and such vectors must lie on 
some hypersurface (mass-shell) given by: M ± (x, k) G(x, k) 4 = — m 2 , which 
must be always connected. 

These conditions follow trivially from the conditions we imposed on any physi- 
cally meaningful Finslerian geometry 

Metric structure construction 

Take the given bi-metric situation and define 

F(x,v) 4 = F + (x,v) 2 F_(x,v) 2 , (1.95) 

with 

F±(x,v) 2 = g±(x)v a v b . (1.96) 

The generic Finsler pseudo-norm is in this case expressed exactly in the same way 
as was the co-Finsler pseudo-norm in the previous case, by the function F(x, v) 2 = 
M 2 (x,v) F(x,v) 4 , where M 2 (x,v) is again an arbitrary function fulfilling exactly 
the same conditions as the Mi(x, k) function, we just have to exchange V* for V. 
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Interconnecting the metric structure with the co-metric structure 

Now the last step in putting constraints on geometric construction is to relate the 
Finsler and co-Finsler structures. The items which are as yet undetermined are the 
functions M±(x, v), and M 2 (x, k). So they have to fulfill the last condition, which is 
the condition of forming the full united geometry This condition is given by the 
Legendre transform relation: 

n 2 

= F(x,v) 2 . (1.97) 



G i.w*,«)» 



In our language, the function M 2 (x, v) must be connected to M\{x, k) by: 
^ / d[M 2 (x,v)F(x,v) 4 ] \ G ( d[M 2 (x,v)F(x,v)T * 



= M 2 (x,v) F(x,v)\ (1.98) 

Thus we see that if we have somehow found an appropriate M 2 (x,v), then Mi (x, k) 
will be uniquely determinecQ. 

To find solutions of such an abstract (although precisely formulated) mathe- 
matical exercise is not an easy task. At this early stage it is quite doubtful whether 
general solutions exist. 

1.6 Conclusions 

In conclusion: Our "no-go" theorem suggests that while pseudo-Finsler space- 
times are certainly useful constructs, in bi-refringent situations it does not appear 
possible to naturally and intuitively construct a "unified" pseudo-Finsler space- 
time such that the pseudo-Finsler metric is null on both "signal cones", but has 
no other zeros or singularities — it seems physically more appropriate to think of 
physics as taking place in a single topological manifold that carries two distinct 
pseudo-Finsler metrics, one for each polarization mode. 



4 It holds both ways, also Mq,(x, k) is uniquely determined by Mi (x, k). 
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This means that in the case of high-energy Lorentz violations it is highly doubt- 
ful whether one has to follow the idea to try to formulate the theory geometrically 
There is no analogy between this case and the very successful Minkowski geomet- 
ric formulation of special theory relativity This leads to a considerable concern 
about how the equivalence principle might be sustainable if we admit the high- 
energy Lorentz violations. 

Physically this might suggest that high energy Lorentz violations - if they ac- 
tually occur in nature - should be "universal". That is, all particles should see 
the same Lorentz violation. In this case one could at least have a single "signal 
cone", avoid all the problems that were demonstrated in this chapter and so have 
a reasonable chance of satisfying the Einstein equivalence principle. 



Chapter 2 

Analytic results for highly damped 
quasi-normal modes 

2.1 Introduction 

Brief introduction into the topic Black hole quasi-normal modes are physically 
intuitive gravitational perturbations of various types of black hole spacetimes. 
Take black hole spacetimes with spherical symmetry, where the metric is of the 
form: 

9^ = ~f(r)dt 2 + f(r)- l dr 2 + r 2 dfi 2 . (2.1) 

One decomposes the general perturbation into a tensorial generalization of spher- 
ical harmonics (for a detailed introduction into the topic see [123]): 



oo I 10 



v = EE E^j^){(^u(m)} . (2.2) 



Z=0 m=—l n=l 



They split into scalar, vector and tensor perturbations. Moreover we can split 
them with respect to parity into two sets: axial and polar perturbations. Both fol- 
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low the one-dimensional equation!^ 

^|M - ^|^> + V(^ lm{ t, X ) = 0. (2.3) 

Here we transformed the coordinate r to the tortoise coordinate x = f -rf\- The 
V(a;) function is called the Regge-Wheeler potential (axial perturbations) or Zerilli 
potential (polar perturbations). Both Regge-Wheeler and Zerilli potentials (in gen- 
eral different) depend on the background geometry, the spin of the perturbation, 
and the wave-mode number of the perturbation. 

If we are interested only in perturbations with the "harmonic" time depen- 
dence ^i m (t,x) = e mt ifji m (x), we obtain: 

^^ - (V(x) - u 2 )^ lm ^) = 0. (2.4) 

Despite the fact that this equation formally reminds us of the Schrodinger equa- 
tion for L 2 Hilbert space self-adjoint operator eigenf unctions, it should be kept in 
mind that here the whole situation is very different. As a result of this fact the cu 2 
"eigenvalues" are in general non-real numbers. 



The quasi-normal modes (QNMs) are solutions of (|2.4|) with the boundary con- 
ditions giving purely outgoing radiation: 

ipim{x) -> C ± e ±iu)X x -> T oo. (2.5) 

These particular perturbations are of a general interest because: 

• The boundary conditions are the physically intuitive ones. 

• It can be proven H123I1 that after some time scale these perturbations become 
dominant within arbitrary black hole perturbation. 

• In the field of quantum gravity, there exists Hod's conjecture ||77|, and more 
recently Maggiore's 111 1011 conjecture, concerning the connection between the 
highly damped QNMs and the black hole area spectrum. 



1 This equation reminds us of the Klein-Gordon equation with a potential, describing a relativis- 
tic scalar particle scattering. 
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• The QNMs related to asymptotically anti-de Sitter black holes are interesting 
for people working with AdS/CFT correspondence. 

The second point suggests QNMs describe the characteristic "sound" of black 
holes. They are in principle observable in black hole oscillations and ring-down 
phenomena. As we see the modes are characterized by the QNM frequencies u 
(QNFs - quasinormal frequencies). It can be generally proven that for the QNM 
frequencies we have: 

• the quasinormal frequencies to form an infinite but countable set, 

• the ui-s are complex with Im(u) > 0, and hence describe stable perturba- 
tions, 

• they have real parts symmetrically spaced with respect to the imaginary 
axis. 

The literature lists many techniques that have been used to calculate the QNM 
frequencies. The basic ones are: 



WKB inspired approximations H57ll74ll82ll90lH85l; 



phase-amplitude methods EHH; 

continued fraction approximations 111051 11061 11071 ; 



monodromy techniques |J5Tll65llll8llll9 | l ; 



Born approximations B21 ITOITI51IT25I ; 



• approximation by analytically solvable potentials. 

Some of these techniques were used only to calculate the fundamental (least damp- 
ed) QNM frequencies (like the approximation by the analytically solvable Poschl- 
Teller potential H571D , other techniques were used to estimate the asymptotic be- 
havior of highly damped frequencies as well. 
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The basic focus and results of our work The main focus of this work is the an- 
alytic results for highly damped QNM frequencies. These can be obtained by two 
different methods: approximation by analytically solvable potentials, and mon- 
odromy techniques. Unfortunately the first method was previously used only 
to estimate the fundamental frequencies, and one of the main contributions of 
this work is using this method to also explore the highly damped QNMs. Unlike 
the method of approximation by analytically solvable potentials, the monodromy 
technique was used many times to understand the asymptotic QNM behavior. 
There are striking similarities between our results obtained by the analytically 
solvable potentials and the known monodromy results. These allow us to analyze 
the monodromy results in a new way bringing much deeper understanding about 
the behaviour of the asymptotic QNM frequencies encoded in the monodromy 
formulae. 

A problem of special interest is the following: for the highly damped frequen- 
cies the asymptotic behavior 

u n = uo + in ■ gap + 0(n~ m ), m > 0, (2.6) 



(where "gap" denotes some real constant), was often observed |118illl9| . On the 



other hand situations have been observed where this behavior seems to fail. The 
analysis of the equations derived by both the analytically solvable potentials and 
monodromy techniques gives us an indication to when such behavior is to be 
expected. It also tells us how the gap spacing is given. 

The structure of this chapter In the first part of this chapter we analyze the 
highly damped QNMs by using the idea of approximation by analytically solv- 
able potentials. We verify our method by applying it to the Schwarzschild black 
hole, where the results are known and widely accepted. After that we explore the 
much less known Schwarzschild-de Sitter (S-dS) case. As a result of our approach 
we will prove interesting theorems about the highly damped QNM behavior for 
the S-dS black hole. We will also discuss our results in the context of black hole 
thermodynamics. In the second part we explore the complementary set of ana- 



CHAPTER 2. ANALYTIC RESULTS FOR HIGHLY DAMPED QNMS 38 

lytic results obtained by monodromy calculations (related to many different types 
of black holes) and prove that all of the results follow the patterns discovered by 
our approximation. This means we are able to generalize our theorems for almost 
every analytic result presently known. This suggests the behavior discovered is 
very generic also between different black hole spacetimes. 

2.2 Approximation by analytically solvable potentials 

2.2.1 Introduction 

As previously mentioned, one of the ways to derive analytic approximate expres- 
sions for the highly damped QNMs is to approximate the real Regge-Wheeler- 
/Zerilli potential by analytically solvable potentials. In the past this method has 
been used to give a formula for the fundamental QNM frequencies. Ferrari and 
Mashoon ||57|| used the Poschl-Teller (Eckart) potential 

V(x) = V /cosh 2 (ax) (2.7) 

to approximate the Regge- Wheeler potential at the peak (by fitting the V and a 
parameters by the peak height and peak curvature). The QNM frequencies of the 
Poschl-Teller potential are given by the formula 



-'„ ±\/V -^--ia(n + ^). (2.8) 



The widely accepted result (see for example 111191 ) for highly damped QNM 
frequencies of Schwarzschild black hole is 

u n = ±^k + %k (n + ~) + O (n- 1 / 2 ) , (2.9) 

where k is the surface gravity at the black hole horizon. The interesting observa- 
tion is the following: 

• the asymptotic formula (|2.9|) can be obtained by the appropriate fitting of 
Poschl-Teller potential, by setting a = k, 
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• such Poschl-Teller potential qualitatively recovers the behavior of one taio 
of the Regge- Wheeler potential. 

This fact seems not to be a coincidence as one can also turn the logic around. The 
wavepacket formed of highly damped QNMs close to the peak will quickly spread 
out from the peak region to the region of the tails of the potential. That suggests 
the tails are the most important factor in determining the highly damped QNMs. 
There is also a general observation that the wavelength (given by the Re(u)) is 
higher for highly damped modes, so this supports the view that these modes must 
be more sensitive to the asymptotic behavior of the potential. So the general ex- 
pectation is that a good approximation to the tails of the Regge-Wheeler/Zerilli potentials 
should give a good qualitative estimate for the behavior of the highly damped quasinormal 
modes. 

To give this statement an exact meaning it is appropriate to exactly define what 
we mean by a good qualitative vs. quantitative asymptotic estimate. We say that 
sequence of QNFs u ln quantitatively matches with the sequence u 2n (for the asymp- 
totic QNFs), if 

lim (u ln - u 2n ) = 0. (2.10) 

n— >oo 

We say that sequence uj\ n matches qualitatively with the sequence u 2n if 

,. \Uln — U)2n\ -.. l^ln — U)2n\ „ /<-iii\ 

lim : : = lam ; ; = 0. (2.11) 

n— >oo \Uli n \ n— >oo \0J2n\ 

So in the first case the "error" goes to zero, in the second case only the "relative 
error" (error compared to the result) goes to zero. One can observe that if ui n and 
u)2n are both equispaced, then if they match quantitatively they must be equal. On 
the other hand for qualitative matching only the same gap spacing is required (they 
can have different coq modes). The second statement can be proven immediately 



2 By the "tail of the potential" we mean the potential in one of the asymptotic regions, given 
as |ar| >> 0. By the "peak of the potential" we mean the potential in the region near its global 
maximum, typically near x = 0. 
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from the definition 



lim 



I Win - U2n\ |gap(wi n ) - gap(w 2 r. 



>oo \UJi2n\ 



|gap(wi, 



(2.12) 



2n 



where gap(tui n ), gap(o;2 n ) denote the "gap" constants from (|2.6|) related to sequences 
Wi n/ w 2n . Then the condition of qualitative matching implies 



gap ( w in ) - gap(w 2 „) = 0, 



(2.13) 



and nothing else. 

All this means the question whether the tails of the real potentials can be ap- 
proximated by piecewise smooth analytically tractable potentials is of high inter- 
est. As we will see this is the case of 1 and 2 horizon situations. 



One horizon situations 

This is the situation in the Schwarzschild geometry. Since it is asymptotically flat, 
only one side of the potential is approaching the (black hole) horizon. 

For the specific case of a Schwarzschild black hole the tortoise coordinate is 
given by 



dr 
dx 



2m 
r 



x(r) = r + 2m In 



r — 2m 
2m 



and the Regge- Wheeler potential is 



VU(rY) = | l- — 



£(£+!) 2m(l-s 2 ) 



(2.14) 



(2.15) 



Here s is the spin of the particle and £ is the angular momentum of the specific 
wave mode under consideration, with £ > s. Asx— > — oo we have r — > 2m and 



V(x) — > exp 



x — 2m 
2m 



+ !) + (! -s' 
(2m) 2 



V exp(2nx) 



(2.16) 



where k is the black hole surface gravity. The Zerilli potential is more complicated, 
but leads to the same asymptotic behavior. This specific behaviour in terms of 
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the surface gravity generalizes beyond the Schwarzschild black hole and for an 
arbitrary black hole in an asymptotically flat spacetime one has 

{Vn_ exp(— 2/dxl), x — y — oo: 
(2.17) 
V 0+ (2m) 2 /x 2 , x ->■ +oo. 

For the highly damped modes this behavior suggests we shall fit the real poten- 
tials by a proper combination of Poschl-Teller and 1/x 2 (inverted harmonic oscila- 
tor) potential. 

Two horizon situations 

If one turns to asymptotically de Sitter black holes (or more generally any two 
horizon system) the situation is different, but following the same patterns — the 
Regge-Wheeler potential and the Zerilli potential have the asymptotic behaviour 

{Vn_ exp(— 2k-\x\), x —y — oo: 

(2.18) 
Vo+ exp(—2K + \x\), x —¥ +oo; 

where the two surface gravities are now (in general) distinct. It is this situation 
that we will model using a piecewise Eckart potential (Poschl-Teller potential). A 
considerable amount of analytic information can be extracted from this model, 
information which in the concluding discussion we shall attempt to relate back to 
"realistic" black hole physics. 

3 (and more) horizon situations 

There is no simple or practicable way of dealing with three-horizon situations 
using semi-analytic techniques. 

2.2.2 Schwarzschild black hole 
Potential 

We work with the Regge-Wheeler/Zerilli equation 
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- ip"(x) + (V(x) - u 2 )^(x) = 0. (2.19) 

The potential recovering the tails behavior of the given Regge-Wheeler/Zerilli 
potential is 

{V _ sech 2 (x/6_) for x < 0; 
(2.20) 
V 0+ /(x + a) 2 for x > 0. 

Here 6 is defined by b- = k~ x . The potential can be discontinuous at the origin 
and the infinite peak can be placed anywhere in the negative part of the real line, 
so a can be taken to be an arbitrary positive real number. 

Wavefunction 



We are interested in the solutions of the Regge-Wheeler/Zerilli equation (|2.4|) with 
the boundary constraints 

ip + (x -> +oo) -> C + e~ iu>x - V-(z-^-oo) ^C_e + ^ x , (2.21) 

defining the quasi-normal modes. The solutions are (for Vq+ > —1/4) 



V+(x) = C+\/ — [J a+ (u(x + a)) - e- iQ +V_ Q+ (w(x + a))] , (2.22) 

V oo 

where J a+ (x), J- a+ (x) are Bessel functions with 

a + = ^1 + 4^0+ (2.23) 

and 

V>_ = CLe*"* 2 Fi Q + a_, ^ - a_, 1 + ib-tu, l + ^_ 2x/b ) , (2-24) 

where 2-^1 (—) is the hypergeometric function with 
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\ - V b 2 _ for V b 2 _ < 1/4; 
'»- = { (2.25) 

iJv bl - \ for V b 2 _ > 1/4. 



Junction condition 

We need to match the functions ip±(x) and their derivatives at the origin. We can 
combine the two matching conditions in such a way, that one of them will give 
the equation for quasinormal modes, while the other will only relate the normal- 
ization constants C±, hence will not be of pressing interest. The first equation 
we obtain by equating the il>'±(0)/il>±(0) ratios. (In special situations where ^±(0) 
might accidentally equal zero one might need to perform a special case analysis. 
The generic situation is ip± (0) ^ 0, and will prove sufficient for almost everything 
we need to calculate.) 

For the V+(°)AM ) ratio ** holds: 

g+(g) _ 1 
i/>+(Q) 2a 

u J( Q+ -i)(atj) - J( a++ i)(au) - e- ia+7T [j ( _ a+ _i)(oa;) - J ( _ a++ i ) (aa;)] 

2' J a+ (aco) - e- ia + n J- a+ (acu) 

This equation is obtained after using known differential identities for Bessel 
functions (see identity [HT9l in appendix IB .4[) . 

The key step in obtaining the ratio ip'_(0)/ip-(0) is to calculate the logarithmic 
derivative. By choosing the variable z = 1/(1 + e~ 2x ^ b -), note that x = maps into 
z = 1/2. Then using the Leibnitz rule and the chain rule one has: 



ip'jO) 1 d ln{ 2 Fi (i + a_ ± - a_, 1 + *_w, z) } 



V>-(0) 2fe_ dz 

Invoking the differential identity (|B.8|) in appendix IB.31 we see 



(2.27) 



2=1/2. 



^-(0) 2F1 (h + «_, \ — a-,ib-u,z) 

— iuj 



^_(0) '" 2-Fi (£ + «-, 5 -a_,l + i&_w,z) 



(2.28) 



2=1/2. 



CHAPTER 2. ANALYTIC RESULTS FOR HIGHLY DAMPED QNMS 44 

Now using Bailey's theorem (|B.7|) to evaluate the hypergeometric functions at 
z — > i we have the exact result 



V>_(o) 


r, -p^a_+io;b_ 


- + |)r(- Q -+^- 


■ + !) 


V_(o) 


"" l V 2 


- + i)rr Q -+^- 


■ + i) 



(2.29) 



The asymptotic approximations 

The exact junction condition we wish to apply at x = is 



(2.30) 



^+(0) ^-(0)' 

We see that equating the right sides of (|2.29|) and (|2.26|) gives a messy equation and 
the presence of the Gamma and Bessel functions above makes this exact junction 
condition intractable. Fortunately, from the beginning we are interested only in 
the asymptotic quasi-normal modes. This means Im{uj) 3> and Im(co) ^> Re(cu). 
In such a case we can use asymptotic formulas for both the Bessel and Gamma 
functions. 



For Bessel functions one takes the asymptotic expansion (|B.10|) from the ap- 
pendix EH In the dominant (zero-th) order, for Re(x) bounded and Im(x) — > 
oo, the functions P(a, x) and Q(a, x), (see IB.11I and IB.12I in appendix IB.4[) , be- 
have as P(a,x) — > 1, while Q(a,x) — > 0. Hence for x ~ Im(x) ^> the 



expansion (|B.10[) gives: 



j a(l)w J „»(,_). (2.31) 



By substituting this into (|2.26|) , after some additional calculations using trivial 



trigonometric relations we obtain a very simple asymptotic formula: 

g| K __L +/m(w)K/m(w) . (2 .32) 

The case (|2.29|) is slightly more difficult. In this case consider the following: if 
to has a large positive imaginary part, then the Gamma function arguments above 
tend towards the negative real axis, a region where the Gamma function has many 
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poles. This is computationally inconvenient, and to obtain a more tractable result 
it is extremely useful to use the reflection formula (|B.4[) of appendix IB.2l to derive 



- ot—+iu)b- 



+ 



r 



■ a—+iu)b- 



+ *: 



r(i- 


a-+iu)b— \\ 
2 4> 


sm(7r[ 2 


+il) 


r(i- 

r(- 


a-+iuib- 3\ 
2 4/ 

a_+iwfc_ , 3\ 
2 4/ 


sm(7r[ 2 
sin(7r[ 2 


h-ID 

-+JD 



r( 
r( 



ct—+iu)b- 



Slll 7T 



-iu)b- 



+ iJ + i 



a_ +iuib- 



+ ¥. 



x tan 7r 



a_ + io;6_ 



+ 



(2.33) 



This leads to the exact result 

— iu)b- 



V>-(0) 

V_(o) 



2 r(- 

x tan 



+ 



-iujb- 



-icub- , 1 



-iuib- 



71 



a. 



4' 

+ io;6_ 



+ 



+ 



tan 7r 



-a_ + iub- 1 
2 + 4 



(2.34) 



If u has a large positive imaginary part, then the Gamma function arguments 
above now tend towards the positive real axis, a region where the Gamma func- 
tion is smoothly behaved — all potential poles in the logarithmic derivative have 
been isolated in the trigonometric functions. We can also use one of the trigono- 
metric identities (|B.1|) of appendix IB. II to rewrite this as 



C(0) 



2 r 



—ct-—iu)b- 



+ 



-iu)b- 



+ 



COS 7T<X 



^_(o) b. r(- 



-iu)b- 



+ h)ri 



-iujb- 



X 



.) - cos(7i[iujb^ + 1/2]) 



+ 



COS 7TQ!_ 



cos(7r[ioob_ + 1/2]) 



which we can rewrite (still an exact result) as 

-p/ a-—iojb 



C(Q) 

^_(0) 






-iub- 



+ 



+ 4) cos(7ra_) + sin(i7ro;6_^ 



-iuib- 



+ *)H 



-iuib- 



+ 



COS 7TQ!_ 



sinu7ru;fe_ 



(2.35) 



We have already seen how to eliminate the Bessel functions in the approxima- 
tion that Im(<jj) is very large and Re(u) bounded. Now we have to do the same 
with Gamma functions. Fortunately this is possible as well. As long as we are 
primarily focussed on the highly damped QNFs (Im(u) — ► 00, Re(u) bounded) 
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we can employ the Stirling approximation in the form (|B.6|) indicated in appendix 
IB.2l to deduce 



Y\( —ct——iuib— 


+ i) 


1 y 2 


-p/ — a—— iuib— 


+ 1) 


1 \ 2 



-cv_ — iub^ 1 



Im{oS)b- 



x 

4 

1 + 



1 + 



Im(ub- 



Im(tub- 



This allows us to approximate the junction condition by 

cos(7ra_) + sin(i7ru;6_) 



leading directly to 



and hence 



cos (7ra_) — sin(z7ra;6_ 



sin(iTTub-) 



in 
u) = - — = inn. 
b- 



(2.36) 



(2.37) 



(2.38) 



(2.39) 



Discussion 

This result means that the QNMs are purely imaginary and equispaced with the 
gap being given by the surface gravity at the black hole horizon. This result can be 
considered to qualitatively (but not fully quantitatively) match with the well known 
result (1Z91 - 



2.2.3 Schwarzschild - de Sitter black hole 
Potential 

The model we are interested in investigating is 

- ip"(x) + (V(x) - io 2 ) ip(x) = 0, 



(2.40) 
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with 

{Vo- sech 2 (x/b-) for x < 0; 
(2.41) 
V 0+ sech 2 (a;/6 + ) for x > 0. 

Here we again take b± = k^. 1 . We will allow a discontinuity in the potential at 
x = 0. The standard case that is usually dealt with is for 



Vo- = V 0+ = V ; b_=b + = b; V(x)= ,__ l2 ° • (2.42) 



Vo_ 

cosh 2 (x/b) 



A related model where V - = V 0+ = V but b + ^ &_ has been explored by 



Suneeta [1591, but our current model is more general, and we will take the analysis 



much further. 



Wavefunction 



We start again by imposing quasi-normal boundary conditions (outgoing radia- 
tion boundary conditions) 

V>+(x -)■ +oo) -► C+e-™*; ^_(x -> -oo) -> C_e +iwx . (2.43) 

On each half line (x < 0, and x > 0) the exact wavefunction (see especially page 
405 of the article by Beyer [23]) is now: 

Mx) = C ± eTiUJX 2F1 Q + a±, \ - a±, 1 + ib±u, x + ^ ±2x/b± ) , (2-44) 

where 



<l-V bi for F 6|<l/4; 
n± = { (2.45) 



/i/y 6|-i for V Q b 2 ± >l/A. 



The same as in previous section, "all" we need to do is to appropriately match 
these wavefunctions at the origin. 
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Junction condition 

The junction condition is again 

<Mo) i>-((>) ( 

and the same remarks related to zero values of ^±(0) as in the previous section 
hold. The ratio ijj±(0)/ip±(0) can be rewritten through Gamma functions as: 



V4(o) 


2 r( « ± + 2 ^ ±+ | )r( -« ± +^ ±+ |) 


^±(0) 


T b ± r( Q ± + f 6 ± + 1) r(- Q± +^± + 1) 



(2.47) 



The asymptotic approximations 



Here both sides of the junction condition contain Gamma functions in the same 
way as the right side of the junction condition from the previous section. Hence 
all the analysis can be repeated (rewriting some of the Gamma functions through 
the basic Gamma functions relations and using the Stirling approximation for the 
asymptotic modes). Following exactly the same arguments one obtains the fol- 
lowing approximation to the junction condition: 

cos(7r« + ) + sm(i7iub + ) cos(7ra_) + sin^incub-) 

cos(7ra + ) — sm(iiiub + ) cos(7ia) — sm(^7^UJb) , 

which is accurate up to fractional corrections of order O (l/Im(ub±)). This is now 
an approximate "quantization condition" for calculating the QNFs being asymp- 
totically increasingly accurate for the highly-damped modes. 

QNF condition 

The asymptotic QNF condition above can, by cross multiplication and the use of 
trigonometric identities, be rewritten in any one of the equivalent forms: 

sin(— i-Kub + ) sm{—moob) = cos(7r« + ) cos(7r«_); (2.49) 

smh.(7iub + ) smh(iiub) = — cos(7ia + ) cos(7ia); (2.50) 
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cos(— inuj\b + — b-]) — cos(— i7iu[b + + b-]) = 2 cos(7ra + ) cos(na); (2.51) 

cosh(7TUj[b + — b-]) — cosh(7ru[b + + &_]) = 2 cos(7ro; + ) cos(7r«_). (2.52) 

Which particular form one chooses to use is a matter of taste that depends on ex- 
actly what one is trying to establish. It is sometimes useful to split the asymptotic 
QNF condition into real and imaginary parts. To do so note 

cos(A + iB) = cos(A) cos(iB) — sin(A) sm(iB) 

= cos(A) 00811(5) - i sin(A) smh(B), (2.53) 

so that 

cos(— i7iu[b + — bJ\) = cos( y Im(u)7i[b + — bJ\) cosh(_Re(o;)7r[6 + — &_]) 

-ism(Im(u)ii\b + - 6_|) sinh( J Re(a;)7r|6 + -6_|). (2.54) 

Therefore the asymptotic QNF condition implies both 

cos(Im(cj)7r[b + — &_]) cosh.(Re(u)ir[b + — b-]) 

= cos(Im(u)7r[b + + 6_]) cosh.(Re(u)7i[b + + 6_]) + 2 cos(7r« + ) cos(7ra_), (2.55) 

and 

sin(/m(a;)7r|6 + — 6_|) sinh(i?e(w)7r|6 + — 6_|) 

= sm(Im(uj)7i[b + + fe_]) sinh( J Re(w)7r[6 + + &_]). (2.56) 

We shall now seek to apply this QNF condition, in its many equivalent forms, to 
extract as much information as possible regarding the distribution of the QNFs. 

Some general observations 

We shall start with some general observations regarding the QNFs. 

1. Note that the a± are either pure real or pure imaginary. 

2. Consequently cos(7ra + ) cos(7ra_) is always pure real G [— 1, +oo). 
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3. If cos(7ra + ) cos(7ra_) > 0, then there are no pure real QNFs. 

Proof: Consider equation (|2.50[) and note that under this condition the LHS 



is positive while the RHS is negative. 

4. If cos(7r« + ) cos(7ra_) < 0, then there is a pure real QNF. 

Proof: Consider equation (|2.50[) and note that under this condition the RHS 
is positive. The LHS is positive and by continuity there will be a a real root 

uj G (0, oo). 

5. If cos(7ra + ) cos(7ra_) > 1, then there are no pure imaginary QNFs. 



Proof: Consider equation (|2.49|) and note that under this condition the LHS 
< 1 while the RHS > 1. 

6. There are infinitely many pure imaginary solutions to these asymptotic QNF 
conditions provided cos(7tq; + ) cos(7ra_) < Q(b + ,b-) < 1; that is, whenever 
cos(na + ) cos(7ra_) is "sufficiently far" below 1. 

Proof: Define 

Q(6 + , 6.) = max ( cos (l"H & + Z M) - cos(|u;|7r[6 + + 6-]) 1 ^ L ^ 



Then by inspection equation (|2.51[) will have an infinite number of pure 
imaginary solutions as long as 

cos(7ra + ) cos(7r«_) < Q(b+,b-). (2.58) 

7. For any purely imaginary u there will be some choice of b±, a± that makes 
this a solution of the asymptotic QNF condition. 

Proof: Consider the specific case 

u, = %±H = ' (a ; + 5 ) , (2.59) 

b + o_ 

and note this satisfies the QNF condition but enforces only two constraints 
among the four unknowns b±, V ±. 
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Rational ratios for the f alloff 

Suppose b + /b_ is rational, that is 

r = -eQ, (2-60) 

and suppose we now define 6* by 

b + = p+b*; b- =pJb*\ 5* = hcf(6 + ,6_), (2.61) 

where "hcf" means "highest common factor". (p + , p_ are relatively primely) Then 
the asymptotic QNF condition is given by 

sin(— iu)7rp + b*) sin(— iunrp-b*) = cos(7ra; + ) cos(7ro;_). (2.62) 

If w* is any specific solution of this equation, then 

i2n ( 1 1 \ 

uj n = w* + -. — = w* + i2n 1cm — — (2.63) 

o* \b+ b_J 

will also be a solution. (Here "lcm" stands for "least common multiplier".) But are 
these the only solutions? Most definitely not. For instance, consider (for rational 
b + /b-) the set of all QNFs for which 

Imiu) < ^, (2.64) 

and label them as 

uv a e {1,2,3... N}. (2.65) 

Then the set of all QNFs decomposes into a set of families 

S« = "o,a + ^; ae{l,2,3...iV}; n G {0,1,2,3...}; (2.66) 

where N is yet to be determined. But for rational &+/&_ we can rewrite the QNF 
condition as 

cos(— iojTTb it \p + — p_|) — cos(— iumb i „\p + + p~]) = 2 cos (vra + ) cos (na_) . (2.67) 



3 Later in this chapter will the symbols Pi/pj automatically mean a rational number given by 
the relatively prime integers pi, pj. 
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Now define z = exp(uirb*), then the QNF condition can be rewritten as 

z \p+-p-\ + z -\p+-p-\ _ z \p++p-} _ z -\p++p-] = 4 cos (,r a+ ) cos (^^ ^ (2.68) 

or equivalently 

^2[p + +p_] _ ^|p + -p_| + [p + +p_] 

-4 cos (7ra + ) cos (vra_) ^-^-l - ^-|p+-p-l+b++P-] + 1 = 0, (2.69) 

that is 

2 2[p + + P _] _ Z 2 P+ _ Z 2 P _ _ 4 cog ^ a+) cog ( 7rQ ,_^ ) Z [ P++P -] + x = o. (2.70) 

This is a polynomial of degree iV = 2(p + + p_), so it has exactly N roots z a . In 
terms of equi-spaced families of QNM the situation is a little bit more tricky. If 
p + ■ p_ is odd then the gap spacing is only -£- and the generic number of different 
families is only p + + p_, while if p + ■ p_ is even, then the gap is given by f- and the 
generic number of families is 2{p + + p_). This can be summarized as: 

Theorem 2.3. Take b+/b- = p+/p~ G Q and 6* = hcf (b+, 6_). Tfoe QNFs are, zyff/z ffoe 
imaginary part of the logarithm lying in [0, 2ir) given by 

\n(z a ) i2n J 2(p+ + p _)\ 

1 — — a G < 1,2,3. . . JV < > n G {0, 1,2,3, ... }, 



a;. 



/i. a 



7T&* gh { g 

(2.71) 

w/f/z g — 2 /or p + • p_ orfrf anrf <? = 1 /or p + ■ p_ euen. 

So for rational b + /b_ all modes form equi-spaced families, all families have 
the same gap spacing and are characterized by distinct offsets ln(z )/(7r&*). That 

is: Arbitrary rational ratios of b + /b_ automatically imply the u n = offset + in ■ gap 
behaviour. 

Irrational ratios for the falloff 

Now suppose b + /b- is irrational, that is 

K = hcf (&+,&_) = 0. (2.72) 
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Then all of the "families" considered above only have one element 

^o,a a e {1,2, 3 . . . oo}. (2.73) 

That is, there will be no "pattern" in the QNFs, and they will not be regularly 
spaced. (Conversely, if there is a "pattern" then b + /b- is rational.) Stated more 
formally, it is possible to derive a theorem as below. 

Theorem 2.4. Suppose we have at least one family of equi-spaced QNFs such that 

uj n = u + in- gap , (2.74) 

then b+ lb- is rational. 



Proof. If we have a family of QNFs of the form given in equation (|2.74|) then we 
know that Vn > 

cos(— iuoir\b + — 6_| + nKir\b + — fe_|) — cos(—iuo7r[b + + 6_] + nK7i\b + + b-\) 

= cos(— ioj 7i\b + — b-\) — cos(— ioj n[b + + &_]). (2.75) 

Let us write this in the form Vn > 

cos(A + nJ) - cos(B + nL) = cos(A) - cos(B), (2.76) 

and realize that this also implies 

cos{A + [n + 1] J) - cos(B + [n + 1]L) = cos{A) - cos{B), (2.77) 

and 

cos(A + [n + 2]J)- cos(B + [n + 2]L) = cos(A) - cos(B). (2.78) 

Now appeal to the trigonometric identity (based on equation (|B.3|) ) 



cos(A + [n + 2] J) + cos{A + nJ) = 2cos(J) cos{A +[n + l]J), (2.79) 

to deduce 

cos(J) cos(A +[n + l]J)- cos(L) cos(B + [n + \}L) = cos(A) - cos(B). (2.80) 
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That is, Vn > we have both 

cos(A +[n + 1] J) - cos(B + [n + l]L) = cos(A) - cos(B), (2.81) 

and 

cos(J) cos(A + [n + 1] J) - cos(L) cos(B + [n + 1]L) = cos(A) - cos(-B). (2.82) 

The first of these equations asserts that all the points 

(cos(A+[n + l]J), cos(B + [n + 1]L)) (2.83) 

lie on the straight line of slope 1 that passes through the point (0, cos B — cos A). 
The second of these equations asserts that all the points 

(cos(A+[n + l]J), cos(B + [n + 1]L)) (2.84) 

also lie on the straight line of slope cos(J)/cos(L) that passes through the point 

(0, [cos B — cos A]/ cos L). We then argue as follows: 

• If cos J 7^ cos L then these two lines are not parallel and so meet only at a 
single point, let's call it (cos A*, cos B*), whence we deduce 

cos(A+ [n + l]J) =cosA; cos(B + [n+ l]L) = cos5*. (2.85) 

But then both J and L must be multiples of 2n, and so cos J = 1 = cos L 
contrary to hypothesis. 

• If cos J = cos L/ 1 then we have both 

cos(A +[n + 1] J) - cos(B + [n + 1]L) = cos(A) - cos( J B), (2.86) 

and 

cos(J) [cos(A + [n + 1] J) - cos(5 + [n + 1]L)] = cos(A) - cos(S). (2.87) 

but these are two parallel lines, both of slope 1, that never intersect unless 
cos( J) = 1. Thus cos J = 1 = cos L contrary to hypothesis. 
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• We therefore conclude that both J and L must be multiples of 2% so that 
cos J = 1 = cos L (in which case the QNF condition is certainly satisfied). 

But now 

b + + o_ L 
and therefore 

b f e Q. (2.89) 

That is: Rational ratios ofb + /b- are implied by the oo n = offset + in gap behaviour. 

□ 

Explicit examples 

• If b + = b_ = &*, but we do not necessarily demand a + = a_, then the 
asymptotic QNFs are exactly calculable and are given by 

icos _1 {l — 2 cos(7r« i ) cos(7ra_)} in 
UJ n = 1 I ^ ^ + -. (2.90) 

Proof: The asymptotic QNF condition reduces to 

1 — cos(—i2 / K0jb*) = 2 cos(7ro; + ) cos(7ro;_), (2.91) 

whence 

cos(— i2n(jb*) = 1 — 2 cos(7ra + ) cos(7ro;_). (2.92) 

This is easily solved to yield 

— i2nuj n b if = cos -1 {1—2 cos(7ra + ) cos(7ra_)} + 2mr, (2.93) 

whence 

icos -1 {1 — 2 cos(7ra + ) cos(7r«_)} in ,*„.. 

UJ n = 1 ) + t-. (2.94) 

Comment: These QNFs are pure imaginary for cos(7ra + ) cos(7ra_) < 1, and 
off-axis complex for cos(7ra + ) cos(7ra_) > 1. We can always, for convenience, 
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choose to define Re(cos~ 1 (x)) G [0, 27r); then because cos _1 (-) is double val- 
ued we see 

9 G [0, 7r] and cos# = a; implies cos _1 (:r) = {0, 7r — 0}, (2.95) 

G [ff, 27r) and cos0 = x implies cos -1 (x) = {9, 3n — 9}. (2.96) 

With this notation 

in i cos -1 {1 — 2 cos(7rai+) cos(7ra_)j 
w„ = wo + 7-; w = i I +J ^ J -L; (2.97) 

with < Im(u>o) < 1/6* and n G {0, 1, 2, 3, . . . }. Note that because of the 
double valued nature of cos -1 (-) there are actually two branches of QNFs 
hiding in this notation — which we will need if we wish to regain the known 
standard result when we specialize to a_ = a + . (We shall subsequently 
generalize this specific result, but it is explicit enough and compact enough 
to make it worthwhile presenting it in full. Furthermore we shall need this 
as input to our perturbative analysis.) 

• If b + — 35_, that is b + = |6* and 5_ = |5*, but we do not necessarily demand 
a + = a_, then the asymptotic QNFs are calculable and are given by 

i , / 1 ± \/9 — 16 cos(7ra + ) cos(-n-a)) \ 2in ,„„„.> 

., = - cos- (_* ~ j + V P ' 98) 

Proof: To see this note that in this situation \b+ — 6_| = fo* = (6+ + 6-)/2. 
Therefore the QNF condition reduces to 

cos(— i7rujK) — cos(—i27iujb*) = 2 cos(7ro; + ) cos(7ra_), (2.99) 

implying 

cos(— iiruub*) — 2cos 2 (— iiruub*) + 1 = 2 cos(7ro; + ) cos(7ra_). (2.100) 

That is 

2cos 2 (— iitub*) — cos(-iiruiK) — 1 + 2 cos(na + ) cos(7ra_) = 0, (2.101) 
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whence 



1 ± y/1 + 8(1 - 2 cos(7ra+) cos(7ra_)) 
cos(— mub*) = , (2.102) 



so that 



1± ^9-16 cos(7r« + )cos(7ra_)) 
cos(— mub*) = , (2.103) 



implying 



i / 1 ± \/9 — 16 cos(7ro; i ) cos(7ra_)) , 
ira n 6 t = cos" 1 ¥ i — ±Z 1 ^Z + n2 7r. (2.104) 



Finally 



% -, / 1 ± \/9 — 16 COs(7T« i ) COs(7T«_)) \ 2l7l 

oj n = -r cos" 1 ^ j — + — . (2.105) 

Comment: This gives us another specific example of asymptotic off-axis com- 
plex QNF's — now with b + ^ b_ . Note that because of the ± and the double- 
valued nature of cos _1 (-) there are actually 4 branches of QNFs hiding in this 
notation. 

This particular trick can certainly be extended to the cubic and quartic poly- 
nomials, for which general solutions exist. 



- The quadratic corresponds to 

&+ + &_ 

- The cubic corresponds to 

b+-b_ 



2; 6 + = 3 6_. (2.106) 



b + + b. 
- The quartic corresponds to 



3; b + = 2b_. (2.107) 



^-^ = 4; b + = \b_. (2.108) 

b + + o_ 3 
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Some approximate results 

A number of approximate results can be extracted by looking at special regions of 
parameter space. 

Case o_ ps b + Suppose o_ ps b+. Then the quantity — i(b+ — b_)uo is slowly vary- 
ing over the range where —i(b + + b_)tu changes by 2n. Let w* be any solution of 
the approximate QNF condition, and define o* = (b + + 6_)/2. Then for nearby 
frequencies we are trying to (approximately) solve 

cos(— iu*7r[b + — b-]) — cos(— iujn[b + + b-]) = 2 cos(7ra + ) cos(7ra_), (2.109) 

that is 

cos(— iou*Tr[b + — b-]) — cos(—iui2'Kb*) = 2 cos(7ra + ) cos(7ra_), (2.110) 

and the solutions of this are approximately 

in i • i r ii °+ + 6_ _ . . . . 

w n ps w* + — valid for \n\ <C — r . (2.111) 

o* \b + — 0-| 

Thus approximate result will subsequently be incorporated into a more general 
perturbative result to be discussed below. 

Case 6_ <C 6 + Now suppose o_ <C b+. Then the quantity —ib_oo is slowly varying 
over the range where —ib + co changes by 2n. Let cu* be any solution of the approxi- 
mate QNF condition, then for nearby frequencies we are trying to (approximately) 
solve 

sin(— iunb + ) sin(— iu*7rb-) = cos(7ro; + ) cos(7ra:_), (2.112) 

and the solutions of this are approximately 

uo n ps a>* + -. — valid for \n\ <C -^-. (2.113) 

o+ o_ 
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Case a_ « 1/2 This corresponds to 

\/ -6 2 « 0, (2.114) 

in which case the QNF condition becomes 

sin(— iujnb + ) sin(— iuiirb-) «0. (2.115) 

Therefore one obtains either (the physically relevant condition) 

— iub+ = n =>• cj = - — , (2.116) 
or (the physically uninteresting situation) 

— iub- = n =>- lo — — . (2.117) 

Note: If you go to the limit a_ = 1/2 by setting Vo_ = 0, one sees on physical 
grounds that b_ is irrelevant, so it cannot contribute to the physical QNF. Alterna- 
tively if you hold Vo_ 7^ but drive 6_ — >■ 0, then these QNF's are driven to infinity 
— and so decouple from the physics. Either way, the only physically interesting 
QNFs are uj = in/b + . We explore these limits more fully below. 

Some special cases 

A number of special cases can now be analyzed in detail to give us an overall feel 
for the general situation. 

Case a_ = = a + This corresponds to 

Vo-b 2 _ = l - = V 0+ b 2 + , (2.118) 

in which case the QNF condition becomes 

sm(iunb + ) sm(iuirb^) = 1. (2.119) 
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Let us look for pure imaginary QNFs. (We do not claim that these are the only 
QNFs.) This implies that we must simultaneously satisfy both 



or both 



That is both 



sin(—iuj7ib + ) = sm(—iu7Tb-) = 1, (2.120) 

sin(— iojitb + ) = sm^—iunb-) = — 1. (2.121) 

(2.122) 
(2.123) 





— iub + = 2n + 


+ 


2 ; 


—iub- = 2n_ H — , 

2' 




or both 


— iub + = 2n + 


— 


1 
2 : 


—ibjb- = 2n_ . 

2 




Therefore either 












b + 2n + + \ 


An + + 1 
4n_ + 1 ' 




or 


b + 2n + - \ 4n+ 


- 1 


6_ " 2n_ + \ 


6_ 2n_ - \ 4n_ 


-1 



(2.124) 

In either case we need b+/b- to be rational, so that b + = p + b* and 6_ = p_6*. 
This special case is thus evidence that there is something very special about the 
situation where b + /b_ is rational, as we know from the theorem (|2.3|) . Then either 

- w = — :L r^-; -*w = r- 2 -! (2.125) 



or 



Now write 



2n+ - i 2n_ - ± , 

*w = + , 2 ; -iw = =— ^. (2.126) 



n + = m + + np+; n_ = m_ + np_; (2.127) 

with m + < p + and m_ < p_. (While n G {0, 1, 2, 3, . . . }.) Then in the first case 

(2m + + l n] ,f2m_ + i n] in , n .,„. 

u = M — r^ + r r = M — r^ + r r = w * + r> (Z128) 

[ p + b* o* J y p-b* o* J o* 

while in the second case 

277U — t: n) f 2m_ -^ n] in 



u=i Y-fr JL + T =i T^+f f= w * + r (2 - 129) 

As we already know all this is only a consequence of our general result (|2.3[) . 
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Case Vo_ = We can best analyze this situation by working directly with the 
exact wavef unction. If V Q - = then a_ = 1/2 and 

ib' CO) 
V-(0) = 1; #_(0) = +iw; ^^ = +z W . (2.130) 

The exact QNF boundary condition is then 

iu = -f- l ' f l 1 h r - (2.131) 

But this we can rewrite as 

_ ^ r(^^ + f ) r(^± + j) cos(7ra + ) - sin(-t7ra;6 + ) 
fe + r (- Q +-^ b + + 1) Y( a +- iub + + 1) cos(7ra+) + sin(-nru;6 + ) ' 

This certainly has pure imaginary roots. If we write to = i\ui\ then asymptotically 
(| a; | — » oo) this becomes 

_ cos(7ra + )-sin(7r|a;|6 + ) 
cos(7ra; + J + sin(7r|a;|o + ) 

implying 

IVY 

sin(7r|(jj|6 + ) = 0; =>■ ir\cu\b + = mr; =>- u> = — . (2.134) 

K 

This agrees with our previous calculation for a_ « 1/2 and as expected gives the 
Schwarzschild result (|2.39|) . 



Case 6_ — > This is best dealt with by using a Taylor expansion to show that 

#-(0) 



That is 



+iu + Vb-fo- +0{bi). (2.135) 



^-(0) 

lim - — — — = +iu. (2.136) 

b-^o ^-(0) v ; 

The analysis then follows that for the case Vo- = above, and furthermore agrees 

with our previous calculation for cn_ « 1/2. 
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Case b _ — > oo This is best dealt with by using the Stirling approximation to- 
gether with a Taylor expansion to show that 

^| = zv^^Tv+0(l/6 2 ). (2.137) 

That is 

lim ^=^ = i y / u 2_ v (2.138) 

The exact QNF boundary condition is then 



„ y/ a + +iu)b + _|_ 3\ j-i/ — a + +iai6 + . 



i^uj 2 - Vol = -— v J 2 . . f v * . f . (2.139) 

But this we can rewrite as 

_ « T(^±^ + |) r(^± + f ) cos(7ra + ) - smHTra^j 
6+ r(~ a+ ~^ b+ + 1) r( a+ ~ iujb+ + i) cos(vra+) + sin(-27ra;6 + ) 

If we write a; = z|a;| then asymptotically (\lu\ — > oo, with Vb- held fixed, implying 
that Vq- effectively decouples from the calculation), this becomes 



cos(7ra + ) — sin(7r|o;|6 + ] 
cos(7ra + ) + sin (tt \ ui \ b + ) 

implying 



(2.141) 



sin(7r |o;|£» + ) = 0; =>■ 7r|u;|6 + = ni[\ =>- u = —. (2.142) 

The importance of this observation is that it indicates that for "one sided" poten- 
tials it is only the side for which the potential has exponential falloff that con- 
tributes to the "gap". This again confirms (|2.39|) , since this is the case when cos- 
mological horizon surface gravity goes to and we get the result appropriate to a 
Schwarzschild black hole. 
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Systematic first-order perturbation theory 

Sometimes it is worthwhile to adopt a perturbative approach and to estimate 
shifts in the QNFs from some idealized pattern. Define 

6= 6+ + 6__ A=|6 + -6_|; (2.143) 

and rewrite the asymptotic QNF condition as 

cos(— inuA) — cos(— i27rub) = 2 cos(7ra + ) cos(7ra_), (2.144) 

where we are implicitly holding a± fixed. When A = we have previously seen 
that the QNF are explicitly calculable with 

icos -1 {1 — 2 cos(7ra; + ) cos(7T«_)} in 

Un = ^ + y (Z145) 

Can we now obtain an approximate formula for the the QNF's when A ^ 0? It is 
a good strategy to define the dimensionless parameter e by 

A = 2 e b, (2.146) 

and to set 

u = u + 5u; 5u = 0(e); (2.147) 

so that the asymptotic QNF condition becomes 

cos(— i2rr[uj + Su>]eb) — cos(— i2ir[uj + Su>]b) = 2 cos(7ra + ) cos(7ra_). (2.148) 

Then to first order in e 

cos(— i27icjeb) — cos(— i2ir[u + Su]b) = 2 cos(7ro; + ) cos(7tq;_), (2.149) 

where implicitly this approximation requires e|5o;|6 <^C 1. Subject to this condition 
we have 

cos(— i27r[u + Su]b) = cos(—i27rcjeb) — 2 cos(7ra + ) cos(ira^), (2.150) 



CHAPTER 2. ANALYTIC RESULTS FOR HIGHLY DAMPED QNMS 64 

whence 

— i27r[co n + 5co n ]b = cos -1 {cos(— i2ixCj n eb) — 2 cos(7T« + ) cos(7ra:_)} + 2irn, (2.151) 

so that 

cos -1 {cos(— i2nu n eb) — 2 cos(7ra + ) cos(7ra_)} in 

w„ + du n = % — — h — . (2.152) 

2-ko b 

But we know that the unperturbed QNFs satisfy w n = w + in/b, so we can also 
write this as 

.cos -1 {cos(— i2nCj n eb) — 2 cos(7ra + ) cos(7ra_)} 

du n = i — w - (2.153) 

2-Kb 

Using the definition of co n this can now be cast in the form 

cos -1 {cos(— i2iru n eb) + cos(— i27ru n b) — 1} 

buj n = % — w , (2.154) 

2-Kb 

or the slightly more suggestive 

. .cos -1 |cos(— i2nco n b) + cos(— i2i[Cj ri tb) — 1} 

<5w„ = i — w , (2.155) 

27TO 

which can even be simplified to 

cos -1 {cos(— i27tu b) + cos(— i27iu n eb) — 1} ^_^. 

dw n = z — w - (2.156) 

2nb 

Note that this manifestly has the correct limit as e — > 0. Note that we have not 
asserted or required that Co n tb <ti 1, in fact when n » 1/e this is typically not true. 
(Consequently cos(— i2iiCj n eb) is relatively unconstrained.) Note furthermore that 
Im(5uj n ) < 1/b. 

Discussion 

The key lesson to be learned from our semi-analytic model for the QNFs is that the 
commonly occurring w n = offset + i gap • n behaviour is common but not universal. 
Specifically, in our semi-analytic model the key point is whether or not the ratio 
6 + /6_ is a rational number. 
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Relation to the real physical Regge-Wheeler / Zerilli potentials. The very in- 
teresting question is the following: If the physically relevant results are by this approx- 
imate potentials recovered only qualitatively (so only the correct gap structure is recov- 
ered), does also the multi-family splitting transfer to the physically relevant case? This is 
because all the different families have the same gap structure, so the multi-family 
splitting might be only uninteresting artefact of our approximate model. The an- 
swer to this question comes with the following theorem: 

Theorem 2.4.1. Take for b+/b- = p+/p~ G Q the gap structure as gap = 1/K, where 
b* = b±/p± and assume that the number of equi-spaced families is bounded with respect 
to b + ,b— Then the function u(b+, bJ) is discontinuous at every b+, b~; b+/b- G Q for 
infinite number of QNM frequencies. 

Proof. Take any b+/b- = p+/p~ and arbitrary monotonically growing sequence of 
primes P t . Take a sequence 



It is obvious that 



(bi + , b t .) = ( -^U + , ^-|ft_ ) . (2.157) 



lim(&,+ ,&,_) = (&+,&_). (2.158) 

I— >oo 



But then 



6,. = "* Scd(P, - 2, P+ ) p 159) 

-n — t 
Here "gcd" means "greatest common divisor" and that means gcd(P/ — 2,p + ) is 
from the definition integer. Also it clearly holds that 1 < gcd(P; — 2,p + ) < p + , 
hence it is upper and lower bounded and 6* = b±/p±. But then 

p i 

lim gap(6 i+ , b^) = lim = —J- = oo. (2.160) 

l ^°° l ^°° h gcd(Pi - 2, p + ) 

But if the number of families is bounded with respect to b + , b^ then only finite 

number of modes can be obtained at b + , 6 as a limit 
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lim u N m(bi+,bi-) = u N (b+,bJ). (2.161) 

This proves the theorem. 

D 

But this just means that if we do not want to end up with extremely strongly 
discontinuous function u(b+, b_), we have to accept the fact that the given 1/6* gap 
dependence automatically implies unbounded number of QNM families (with respect to 
b + , bJ), as it was in the case of our analytically solvable potential. 

We also suspect that it might be possible to generalize the model potential 
even further — the "art" would lie in picking a piecewise potential that is still 
analytically solvable (at least for the highly damped modes) but which might be 
closer in spirit to the Regge-Wheeler (Zerilli) potential that is the key physical 
motivation for the current work. (Of course if we temporarily forget the black hole 
motivation, it may already be of some mathematical and physical interest that we 
have a nontrivial extension of the Eckart potential that is asymptotically exactly 
solvable — one could in principle loop back to Eckart's original article and start 
asking questions about tunnelling probabilities for electrons encountering such 
piecewise Eckart barriers.) 

Relation to the black hole thermodynamics. The last point we would like to 
briefly discuss is the relation of our results to the conjectured connection between 
the highly damped QNMs and the black hole thermodynamics. There are two ba- 
sic conjectures: The first conjecture, due to Hod |77|, gives some strong arguments 
supporting the idea that there is a connection between the real part of asymptotic 
QNMs and the quantum black hole area spacing. The basic principle underlying 
this conjecture is Bohr's correspondence principle. The conjecture was also used 
in the context of Loop Quantum Gravity by Dreyer ||53l]. The second conjecture, 
such that it modifies Hod's original proposal, is due to Maggiore 111101 . It solves 
some controversies of Hod's conjecture, (see BUI), and gives the relation between 



CHAPTER 2. ANALYTIC RESULTS FOR HIGHLY DAMPED QNMS 67 

black hole mass spectra and highly damped QNMs as: 



AM « H ■ A^Re{uj n ) 2 + Im(u n ) 2 « h ■ AIm(u n ), n » 0. (2.162) 

Here we naturally assume that Re(co n ) is bounded and Im(uj n ) is growing to in- 
finity But it means, that if we are interested in mass quantum A min M: 

A min M « ft • (Jm(w„) - Jm(w„_i)) , n » 0. (2.163) 

The area quantum for the Schwarzschild black hole one obtains from the formula: 

A min A = 32vrM A min M = 8ttZJ = const. , (2.164) 

(where l p denotes the Planck length). In our case, if we naively extrapolate Mag- 

giore's conjecture to S-dS spacetime and the ratio of surface gravities is rational, 

we obtain black hole mass quanta asQ 

2 
A min M fa h ■ gap(a;„) = h ■ - lcm (k + , «_), n » 0. (2.165) 

9 

Here all the symbols are defined as in the proof of the theorem [23] and in the equa- 
tion (|2.12|) . Now there is a fascinating result by Choudhury and Padmanabhan 



[43J, that clearly "fits" very well into all these ideas. It says that in the S-dS space- 
time there exists a coordinate system with globally defined temperature, (hence 
some kind of thermodynamic equilibrium), if and only if the ratio of surface grav- 
ities is rational. (For a very nice discussion of what might be the physical meaning 
of rational ratios of surface gravities and the role of highly damped QNMs as a 
potential source of information see again |43|.) The global temperature is given as 

T sds = Mi^z). (2.166) 

Then the gap spacing is, in our case: 

l \ 2 1 t \ 2 U tf \ ^P+P- T sds , n - srjs 

gap{oo n ) = - icm(« + , k_j = - p+p- hct{K + , k_j = . (2.167) 

9 9 9 



4 As the reader might have noted: we claim that the expression for the gap in the spacing of the 
QNMs in the asymptotic formula, (the one derived in this section), is an exact result, "unharmed" 
by the fact that we used only approximate potential. 
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This is a different result from the Schwarzschild spacetime, where it is 

gap(u; n ) = 2vrT s . (2.168) 

But the fact that something special happens with both the S-dS spacetime ther- 
modynamics and highly damped QNMs when the surface gravities have rational 
ratio is very interesting. Especially because it allows us to relate the constant gap 
in the QNM spacing and global spacetime temperature, and strikingly they both 
exist under the same condition. Moreover, in terms of Maggiore's conjecture, the 
existence of equispaced families is very interesting, because the gap in the QNM 
spacing multiplied by Planck constant is simply the quantum of black hole mass. 
These considerations suggest that something very interesting is happening, but 
the topic needs clearly further exploration. What is written here has much more a 
character of ambiguous indications than of a well founded ideas, but it gives very 
exciting suggestions for future work. Even more because, as we will show in the 
next section, the link between rational ratio of surface gravities and equispaced 
families of highly damped QNMs is very generic for the multi-horizon black hole 
spacetimes. 

2.4.1 Conclusions 

The arguments that highly damped modes qualitatively depend on the tails of the 
potential were confirmed for Schwarzschild black hole by deriving the behavior 
of (|2.9|) . That was the only straight computational test that could have been done. 
After this the method was used to obtain new results for Schwarzschild-de Sitter 
black holes. One of the nice features of the semi-analytic model related to S-dS 
black holes is that a quite surprising amount of semi-analytic information can be 
extracted, in terms of general qualitative results, approximate results, perturba- 
tive results, and reasonably explicit computations. Our results also might have 
important implications for the black hole thermodynamics. 
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2.5 Monodromy results 

2.5.1 Introduction 

In monodromy approaches one works with an analytic continuation into the com- 
plex radial plane. They are part of a wider approach, the phase integral method 
(for details see [5]). One has to choose branch cuts (as a part of analytic continua- 
tion), identify the singular points of the solutions in the complex plane, locate the 
Stokes and Anti-Stokes lines and calculate the monodromy around the singulari- 
ties Biei^l^lll^lMllM ITMITmiTMITMITITllM^ 

While many technical details differ, both between the semi-analytic and mon- 
odromy approaches, and often among various authors seeking to apply the mon- 
odromy technique, there is widespread agreement that not only the semi-analytic 
approximation, but also the monodromy approaches lead to QNF master equa- 
tions of the general form: 

l>exp(£^)=0. (2.169) 

A=l W Ki J 

Here Ki is the surface gravity of the i-th horizon, H is the number of horizons, 
the matrix Zm always has rational entries (and quite often is integer-valued). The 
physics contained in the master equation is invariant under substitutions of the 
form Zm — > %Ai + (1, ■■•, 1)a^«/ where the hi are arbitrary rational numbers. Either 
Y2i %Ai = 0, or it can without loss of generality be made zero. Furthermore N is 
some reasonably small positive integer. (In fact iV < 2H + 1 in all situations we 
have encountered, and typically N > H.) The C A are a collection of coefficients 
that are often but not always integers, though in all known cases they are at least 
real. Finally in almost all known cases the rectangular N x H matrix Za% has 
rank H, and the QNF master equation is almost always irreducible (that is, non- 
factorizable). We shall first demonstrate that all known master equations (whether 
based on semi-analytic or monodromy techniques) can be cast into this form. Then 
we will generalize the results for rational /irrational surface gravities ratios from 
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the special cases obtained within the approximations by the analytically solvable 
potentials to this more general case (|2.169|) . 



2.5.2 Particular results 

Survey of the monodromy results 

One horizon: In the one-horizon situation there is general agreement that the 
relevant master equation is 



cxp 



( 



7TU) 



V K 



+ 1 + 2cos(ttj) = 0. 



(2.170) 



Unfortunately there is distressingly little agreement over the precise status of the 
parameter j. References H |22l |TT5l EM E2Q1 11211 assert that this is the spin of 
the perturbation under consideration, but with some disagreement as to whether 
this applies to all spins and all dimensions. In contrast in reference 11511 a partic- 
ular model for the spacetime metric is adopted, and in terms of the parameters 
describing this model these authors take 

3 = y " 1- (2-171) 

Here, (and also later in this section), the symbol d denotes spacetime dimension, 
and q is a parameter related to the power with which the general d-dimensional 
black hole metric coefficients 



g{r) 



(2.172) 



behave close to the singularity, (see [51]). Reference [109] asserts that for spin 1 
perturbations 

Be this as it may, there is universal agreement on the form of the QNF master 
condition, and it is automatically of the the form of equation (|2.169|) , with H = 1 
and N = 2 terms. The vector Ca and the matrix Zai are 



Gj. 



+1 

1 + 2 COs(7T?') 



Ul 



-1 





(2.174) 



CHAPTER 2. ANALYTIC RESULTS FOR HIGHLY DAMPED QNMS 



71 



By multiplying through by exp(— ttu>/(2k)) we can re-cast the QNF condition as 

expg) +{l + 2cos(7tf)}exp(~) = 0. (2.175) 

This now corresponds to 



C, 



+1 

1 + 2 COs(7Tj) 



<A\ 



+ 1/2 
-1/2 



(2.176) 



and in this form we have J2a ^ai = 0. (This is one of rather few cases where it is 
convenient to take the Z Ai to be rational-valued rather than integer-valued.) 



Two horizons: For two-horizon situations the analysis is slightly different for 
Schwarzschild-de Sitter spacetimes (Kottler spacetimes) versus Reissner-Nord- 
strom spacetimes. 

• For Schwarzschild-de Sitter spacetimes there is general agreement that the 
relevant master equation for the QNFs is 



( ItU 7TuA ( TXUJ 7YU) 

\ 1 + 2cos(7r?) \ cosh 1 + cosh 



0. 



(2.177) 



We shall again adopt conventions such that k± are both positive. Again, 
there is unfortunately distressingly little agreement over the precise status 
of the parameter j. References DH [221 [1181 EH [1201, HUB assert that this is the 
spin of the perturbation under consideration, but with some disagreement 
as to whether this applies to all spins and all dimensions. In contrast in 
reference [65] a particular model for the spacetime metric is again adopted, 
and in terms of the parameters describing this model they take 



^_1 

2 



(2.178) 



One still has to perform a number of trigonometric transformations to turn 
the quoted result of reference ||65H for d ^ 5 



tanh ( — ) tanh I — ) = ^- . 

k + ) \k_) tan 2 (7rj/2) - 1 



(2.179) 
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into the equivalent form (|2.177|) above. For d = 5 the authors of 111411 assert 
the equivalent of 

{1 + 2 cos(nj)} sinh ( ™ + ™ )+ sinh (™ - ™) = 0. (2.180) 



Reference 111091 again asserts that for spin 1 perturbations 

. _ 2(d-3) 
J d-2 ' 



(2.181) 



Be this as it may, there is again universal agreement on the form of the QNF 
master condition, and converting hyperbolic functions into exponentials, it 
can be transformed into the form of equation (|2.169|) , with H = 2 and N = A 
terms. The vector Ca and matrix Za% are 



cr 



1+2 COs(7T?') 

+ 1 

+ 1 
1 + 2 COs(7Tj) 



J Ai 



+ 1 +1 

+ 1 -1 

-1 +1 

-1 -1 



(2.182) 



Note that we explicitly have J2 A Z Ai = 0. There are two exceptional cases: 
- If j — 2m + 1 with m € Z then 



Ca 



-1 
+ 1 
+ 1 
-1 



JAi 



+1 


+1 " 


+1 


-1 


-1 


+1 


-1 


-1 



(2.183) 



In this situation the QNF master equation factorizes 

0. (2.184) 

This appears to be the physically relevant case for spin 1 particles. The 



sinh ( — ) sinh I — I 
K+J V K - J 



relevant QNF spectrum is that of equation (|2.115[) 
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If cos(7r/') = — ~, which does not appear to be a physically relevant situ- 
ation but serves to illustrate potential mathematical pathologies, then 



C: 





+1 
+1 





Z Ai 



+1 +1 

+1 -1 

-1 +1 

-1 -1 



(2.185) 



But in this situation the top row and bottom row do not contribute to 
the QNF master equation and one might as well delete them. That is, 
one might as well write 



Gj. 



+1 
+1 



J Ai 



(2.186) 



This is a situation (albeit unphysical) where the matrix Zm does not 
have maximal rank. The QNF master equation degenerates to 



. ( T\UJ 1XUj\ 

smn J 

In this situation the QNF spectrum is 



0. 



(2.187) 



Wr, 



K. 



(2.188) 



with no restriction on the relative values of k±. This situation is how- 
ever clearly non-generic (and outright unphysical). 



For Reissner-Nordstrom spacetime one has 111191 



exp 



2ttlo 



2{1 + cos(7r/')} exp 



2ttlu 



{1 + 2cos(ttj)} = 0, (2.189) 



where n+ is the surface gravity of the outer horizon and k_ is the surface 
gravity of the inner horizon. There is again some disagreement on the status 
of the parameter j. Reference | |119 | | now takes j — \ for spin 0, and j — | 
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for spins 1 and 2 (in any dimension). Reference [121] asserts that for general 
dimension 

d — 3 r , „ _ , . 3d — 7 



J 



for spin 0, 2, and 



J 



for spin 1. (2.190) 



2d -5 r ""' " 2d -5 

Be this as it may, there is universal agreement on the form of the QNF master 
condition, and it is automatically of the form of equation (|2.169|) , with H = 2 
and N = 3 terms. The vector Ca and matrix Z Ai are 

+i i r +2 o " 

C A = 2{l + cos(7rj)} ; Z Ai = -2 . (2.191) 

1 + 2cos(ttj) 0_ 

If we multiply through by a suitable factor then we can write the QNF con- 
dition in the equivalent form 



cxp 



2-KLU TTLd 

+ 



+{1 + 2cos(7Tj)} exp 

This corresponds to 

+ 1 
C A =-- 2{l + cos(vrj)} 



+ 2{1 + cos(7tj)} exp 



ITLU 



ITLU 



3k + 3/t- 



J Ai 



-2/3 
-1/3 
-1/3 



TIU 

3k, + 
= 0, 



1/3 

-2/3 

1/3 



2ttu 
3k- 



(2.192) 



(2.193) 



1 + 2C0S(7TJ) 

In this form we now explicitly have J2 A Z Ai = 0. (This is one of rather few 
cases where it is convenient to take the Z Ai to be rational-valued rather than 
integer-valued.) Returning to the original form in equation (|2.189[) , there are 
two exceptional cases: 

- If j = 2m + 1 with m E Z, (this does not appear to be a physically 
relevant situation but again this serves to illustrate the possible mathe- 
matical pathologies one might encounter), then 



Cz 



+1 



-1 



2 








-2 









(2.194) 
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But then (without loss of information) one might as well eliminate the 
second row, to obtain 



C t 



' +1 " 


j %Ai — 


' +2 


-1 








(2.195) 



Furthermore, since K- now decouples, we might as well eliminate the 
second column, to obtain 



C A 



+1 
-1 



Zm 



2 




(2.196) 



The QNF master equation then specializes to 

27ru/ 



cxp 



1 = 0. 



(2.197) 



If cos(ttj) — —\, which does not appear to be a physically relevant situ- 
ation but serves to illustrate potential mathematical pathologies, then 



a 



+i 

+i 

o 



J Ai 



+2 

-2 




(2.198) 



But in this situation the bottom row does not contribute to the QNF 
master equation and one might as well delete it. That is, one might as 
well write 



Ci 



+1 
+1 



J Ai 



+2 
-2 



(2.199) 



We can rearrange the terms in the master equation to have the QNF 
master equation specialize to 



'2nu 2nu\ 
exp | + 1+1 = 0. 



- 



(2.200) 
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This corresponds to 

C A = 



+1 
+1 



J Ai 



+2 +2 




(2.201) 



Note that in this exceptional case Zai is not of maximal rank. In this 
situation the QNF spectrum is 

(2n + l)i k + k_ 



uj r , 



(2.202) 



K + + K_ 

with no restriction on the relative values of n±. This situation is how- 
ever clearly non-generic (and outright unphysical). 



Three horizons: For three horizons the natural example to consider is that of 
Reissner-Nordstrom-de Sitter spacetime. References [121111411 agree that (for d ^ 
5) 

, [txuj ttlu\ r .... f-KOJ nu\ 
cosh + {1 + cos^j)} cosh 1 

\K + K C J \K + K C J 

+2{1 + cos(Trj)} cosh (— + — + —) = 0. (2.203) 

V K - K + K c ) 

Here k± refer to the inner and outer horizons of the central Riessner-Nordstrom 
black hole, while kq is now the surface gravity of the cosmological horizon. All 
these surface gravities are taken positive. In contrast for d — 5 one has 



sinn + {1 + cos(7rjJ| smh — 






+2{1 + cos(7rj)} sinh 



(2-koj 



7TU) 7TUJ 

K + K C 



0, (2.204) 



Again 



J 



d-3 



for spin 0, 2, and 



J 



3d- 7 



for spin 1. 



(2.205) 



2d -5 r 2d -5 

There is universal agreement on the form of the QNF master condition, and con- 
verting hyperbolic functions into exponentials, it can be transformed into the form 
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of equation (|2.169[) , with H = 3 and N = 6 terms. The vector Ca and matrix Z Ai 

are 

+1 

1 + COs(7Tj) 
2{1 + C0S(7T/)} 

±2{1 + cos(ttj)} 

±{1 + COs(7T/')} 
±1 



c, 



J M 



+1 


-1 


" 


+1 


+1 





+1 


+1 


+2 


-1 


-1 


-2 


-1 


-1 





-1 


+1 






(2.206) 



Generically, Za% has maximal rank H = 3. Note that we explicitly have J2a ^Ai 
0. 

The only exceptional case is cos(nj) = — 1 in which case 



C, 



+1 





±1 



J Ai 



1 


-1 


" 


1 


+1 





1 


+1 


+2 


1 


-1 


-2 


1 


-1 





1 


+1 






(2.207) 



But then the 2 to 5 rows decouple and may as well be removed, yielding 



C, 



+1 
±1 



J Ai 



+ 1 -1 
-1 +1 



(2.208) 



The 3 rd column, corresponding to k_, now decouples and may as well be removed, 
yielding 



C, 



+1 
±1 



J Ai 



+ 1 -1 
-1 +1 



(2.209) 



Note that in this exceptional case Za% is not of maximal rank, and the QNF master 
equation degenerates to 



. ' TIUJ TlUj\ 

cosh | = 0, 



' -nu 7TU>\ 
or sinn ( = 0, 

K+ K C J 



(2.210) 
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respectively. In this situation the QNF spectrum is 

(2n + l)z k + k c 



UJ,, 



or 



Wr, 



ITIK+Kc 



2|K + — Kc| |/€_|_ — Kc\ 

respectively, with no restriction on the relative values of k± . This situation is how- 
ever clearly non-generic (and outright unphysical). 



Rewriting the analytically solvable potentials results into our general form 

One horizon For highly damped QNFs the master equation is derived in previ- 
ous section and also in references I1501 I153H , in a form equivalent to 



This means 



f7TUJ\ 

smn I — J = 0. 



C A 



(2.212) 



and 



(2.213) 



>Ai 



1 
-1 



(2.214) 



Two horizons For highly damped QNFs the master equation is derived in pre- 
vious section and also in references 111501 11511 11531 in a form equivalent to 

flTUJ 7TUj\ /tCU 7TuA 

cosh 1 — cosh + 2 cosf 7ra + ) cos(7ra_) = 0. (2.215) 

\K + K- J \ K + K - J 

This means 



C, 



1 

1 

-1 
-1 
2 cos(7ra + ) cos(-7ra_ 



(2.216) 
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and 



J Ai 



1 


1 " 


1 


-1 


1 


-1 


1 


1 





_ 



(2.217) 



Note that a± = . A — -Jr and hence a± = \ is a physically degenerate case corre- 
sponding either to V ± = (in which case the corresponding k± is physically and 
mathematically meaningless), or k± = oo, (in which case the QNF master equa- 
tion is vacuous). Either of these situations is unphysical so one must have a± ^ |. 
This QNF condition above is irreducible (non-factorizable) unless a± = m+ ~ with 
m E Z. This occurs when 

V ± = -m(m + 1) 4, (2.218) 



and in this exceptional situation the QNF master equation factorizes to 

sinh I — ) sinh I — I = 0, 



(2.219) 



and hence takes the same shape as (|2.184[) . 



2.5.3 Rational ratios of surface gravities 

From master equation to polynomial 

First let us suppose that the ratios Rij = k,{/kj are all rational numbers. This is 
not as significant a constraint as one might initially think. In particular, since the 
rationals are dense in the reals one can always with arbitrarily high accuracy make 
an approximation to this effect. Furthermore since floating point numbers are 
essentially a subset of the rationals, all numerical investigations implicitly make 
such an assumption, and all numerical experiments should be interpreted with 
this point kept firmly in mind. 
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Provided that the ratios i2y = K{/kj are all rational numbers, it follows that 
there is a constant k* and a collection of relatively prime integers m; such that 

Kt = — • (2.220) 

mi 

The QNF master equation then becomes 

N / H 



A=l \i=l 



J2 Ca exp [J2 Z M m — = 0, (2.221) 



Now define z = exp(7rw/K*), and define a new set of integers m A = X^=i -^At m i- 
(There is no guarantee or requirement that the fh A be relatively prime, and some of 



the special cases we had to consider in the previous section and in reference H151B 
ultimately depend on this observation.) Then 



N 



J2c A z lflA = 0. (2.222) 



A=l 

This is (at present) a Laurent polynomial, as some exponents may be (and typically 
are) negative. Multiplying through by z^ minin converts this to a regular polynomial 
with a nonzero constant z° term and with degree 

D = m max - m min . (2.223) 

If we write ffiA = friA ~ "Vin then the relevant regular polynomial is 

N 



J2C A z™ A =0- (2.224) 



A=l 



Note that the polynomial is typically "sparse" — the number of terms iV is small 
(typically iV < 2H + 1) but the degree D can easily be arbitrarily large. There are 
at most D distinct roots for the polynomial z a , and the general solution of the QNF 
condition is 

^ a »= K * +2inK*; ae{l,...,D}; n e {0, 1,2, 3, . . . }. (2.225) 

7T 
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If the rfiA are not relatively prime, define a degeneracy factor g = hcf {m^}. Then 
the roots will fall into D/g classes where the g degenerate members of each class 
differ only by the various g-th roots of unity In this situation we can somewhat 
simplify the above QNF spectrum to yield 

u a , n = ^^ + — ; ae{l,...,D/g}; ne {0,1,2,3,...}. (2.226) 

n g 

We again emphasize that behaviour of this sort certainly does occur in practice. 
There is no guarantee or requirement that the fn A be relatively prime, and some of 



the special cases we had to consider in the previous section and in reference ||151|1 
ultimately depend on this observation. 

So let us summarize what we just proved in a theorem: 

Theorem 2.6. Take quasi-normal frequencies to be given by the equation (|2.221|) and, 
take the ratio of arbitrary pairs of surface gravities to be rational. Then the quasi-normal 
frequencies are given by the formula (|2.226[) , where z a are, (in general), D/g solutions of 
the equation (|2.224|) . 

There is a (slightly) weaker condition that also leads to polynomial master 
equations and the associated families of QNFs. Suppose that we know that the 
ratios 

R AB = g|i ZAi/Ki eQ (2.227) 

are always rational numbers. Then it follows that there is a set of integers rh A such 
that 

J2^ = ^, (2.228) 

where the rh A are all relatively prime. (Note k* does not have to equal «;*). This 
is actually a (slightly) weaker condition than R^ = Ki/kj being rational, since it 
is only if Z&i is of rank H that one can derive Rij E Q from R AB E Q. Assuming 
Rab € Q the QNF master equation becomes 

TV , x 

Y,C A exp (m A —) =0. (2.229) 

A=\ \ K * ' 
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This can now be converted into a polynomial in exactly the same manner as pre- 
viously leading to families of QNFs as above. Provided both R AB and R^ are 
rational we can identify R* = K*/g. 

Factorizability 

Now it is mathematically conceivable that in certain circumstances the master 
equation might factorize into a product over two disjoint sets of horizons 



A=l \i=l " 



' N 2 

^C 2A exp 



A=l 




0. (2.230) 



Physically one might in fact expect this if the horizons indexed by i G {1, . . . , Hi} 
are very remote (in physical distance) from the other horizons indexed by i E 
{1, . . . , H 2 }. If such a factorization were to occur then the QNFs would fall into 
two completely disjoint classes, being independently and disjointly determined 
by these two classes of horizon. 

2.6.1 Irrational ratios of surface gravities 

We now wish to work "backwards" to see if the existence of a family of equi- 
spaced QNFs can lead to constraints on the ratios Rij = Ki/Kj. Such an analysis 
has already been performed for the specific class of QNF master equations arising 
from semi-analytic techniques, and we now intend to generalize the argument to 



the generic class of QNF master equations presented in equation (|2.169|) . Let us 



therefore assume the existence of at least one "family" of QNFs of the form: 

u n = ujq + in gap; n G {0, 1, 2, 3, . . . }. (2.231) 



Then we are asserting 

N / H 



A=l \i=\ 



Z Ai tt(w _+ in gap) 



Y.C*>s»[Y. )=* »e{0,1.2,8,...}. (2.232) 
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That is 

^2 S c a exp f ^2 A \, 1 \ ex P ( inn S a P ^2 ~^~ ) = 0; (2.233) 

ne{0,l,2,3,...}. 
We can rewrite this as 

N 

Y^ D A exp (2vrm J A ) = 0; n 6 {0, 1,2,3, . . .}. (2.234) 

A priori, there is no particular reason to expect either the D A or the J a to be real. 

Case 1 

One specific solution to the above collection of constraints is 

N 

^2D A = 0; exp(2vri J A ) = r. (2.235) 

A=l 

Furthermore, as long as no proper subset of the D A 's sums to zero, we assert that this 
is the only solution. To see this let us define 

\ A = exp(2mJ A ); M AB = (A^) 6 " 1 ; A,B e {1,2,3, ... ,N}. (2.236) 



Then M AB is a square N x N Vandermonde matrix, and then equation (|2.234[) 
implies 



N 



J2D A M AB = 0, (2.237) 

A=l 

whence det(M AB ) = 0. But from the known form of the Vandermonde determi- 
nant we have 

det(M AB ) = H (X A - \ B ) = 0, (2.238) 



A>B 
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implying that at least two of the A^ are equal. Without loss of generality we can 
shuffle the A^'s so that the two which are guaranteed to be equal are Ai and A 2 . 
Then equation (|2.234[) implies 

N 

(D 1 + D 2 )Af- 1 + J2 Da{\a) B - 1 = 0; B e {1, 2, 3, . . . , N - 1}. (2.239) 

A=3 

But by hypothesis D x + D 2 ^ 0, so this equation can be rewritten in terms of a 
non- trivial reduced (N — 1) x (N — 1) Vandermonde matrix, whose determinant 
must again be zero, so that two more of the Aa's must be equal. Proceeding in this 
way one reduces the size of the Vandermonde matrix by unity at each step and 
finally has 

X A = r, (2.240) 

as asserted. We then see 

ln(V) 
J A = -i -^ + m A ; m A G Z. (2.241) 

Expressed directly in terms of the surface gravities this yields 

Yz Al ^ = -i l ^ + 2m A] m A EZ. (2.242) 

By assumption, we have asserted the existence of at least one solution to these 
constraint equations. (Otherwise the family we used to start this discussion would 
not exist.) We have seen that we can choose to present the master equation in such 
a manner that Y^,a=i ^ m = ®' ^ u ^ then 

= -i 311 N + 2 V m A ; m A e Z. (2.243) 

A=l 



This implies that 



l l l^ = 2^3A = q eQ. (2.244) 

71 J\ 



That is, there is a rational number q such that 



H 



8=1 



Y^Z M — = q + 2m A ; qeQ; m A eZ. (2.245) 
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This is already enough to imply that the ratios Rab are rational. If in addition 
Zai is of rank H then, (either using standard row-echelon reduction of the aug- 
mented matrix, or invoking the Moore-Penrose pseudo-inverse and noting that 
the Moore-Penrose pseudo-inverse of an integer valued matrix has rational ele- 
ments), we see that for each horizon the ratio (gap)/«$ must be a rational number, 
and consequently the ratios R^ = Ki/Kj must all be rational numbers. 
That is: 

Theorem 2.7. Take quasi-normal frequencies to be given by the equation (|2.221[) , and 
suppose that we have a family of QNFs as described by equation (|2.231|) . If no proper 
subset of the Da's from (|2.221|) sums to 0, and the matrix Z At from the same equation is 
of rank H, then the ratio of arbitrary pair of surface gravities must be rational. 

Case 2 

More generally, if some proper subset of the Da's sums to zero, subdivide the N terms 
A e {1, 2, 3, ... , N} into a cover of disjoint irreducible proper subsets B a such that 

N 

Y^D A = 0. (2.246) 

AGBa 



Then the solutions of equation (|2.234[) are uniquely of the from 



exp(2mJ A eB a ) = -Wsa = r a . (2.247) 



It is trivial to see that under the stated conditions this is a solution of equation (|2.234[) , 
the only technically difficult step is to verify that these are the only solutions. 
One again proceeds by iteratively using the Vandermonde matrix M A b = (^a) b ~ 1 
and considering its determinant. Instead of showing that all of the Xa's equal 
each other, we now at various stages of the reduction process use the condition 
YIiA&b Da = to completely decouple the corresponding \AeB a = r a from the 



remaining \A<£B a - Proceeding in this way we finally obtain equation (|2.247|) as 
claimed. 
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We then see 

JAet3 a = -i -7T- + m AeBa ; m AeBa G Z. (2.248) 

Expressed directly in terms of the surface gravities this yields 



H 



gap . ln(r a 



V Z AeBa AJ± = -i -^L + 2m AeBa ■ m AeBa G Z. (2.249) 

^— ' Ki 7T 

i=l 

With the obvious notation of a(A) denoting the index of the particular disjoint set 
B a that A belongs to, we can write this as 

v^ gap . ln{r a(A) } «, m . 

> z^i = — « — — h 2myi; m^ G Z. (2.250) 

*-^ Ki n 

i=i 

This result now is somewhat more subtle to analyze. Let A and i? both belong to 
a particular set B a . Then 

y^{Z Ai - Z Bi ) = +2{m A -m B }; m A ,m B e Z; A,BeB a . (2.251) 

., Ki 

i=l 

That is 

gap= *{™A-m B } _ gapei? _ A ,BEB a ; (2.252) 

K ■ 
1=1 

so we see that the gap is real. (Furthermore, the gap is seen to be a sort of "integer- 
weighted harmonic average" of the Ki.) But reality then implies that r a = e l( ^ a so 
that 

} Z Ai = — — + 2m a] rn A G Z. (2.253) 



By using J2 A %Ai = we see that 



V^ 0a|#a| ~ ,,, OC/1 x 

- = }, -a? — e ^' (2.254) 

7T ^^ A/7T 
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so that 

E gap <j>a(A)-(f> of . „ rt „,„ 

Z^j = — — h 2 mi - m ; rriA G Z. (2.255) 

Kj 7T 

i=l 

Unfortunately in the general case there is little more than can be said and one has 
to resort to special case-by-case analyses. One last point we can make is that even 
though in this situation the Rij = Ki/kj are sometimes irrational we can make the 
weaker statement that 



V"^ \Zai ~ Zbi} 



K. 



E{Zci — Z Di } 
K ■ 

1=1 



e Q; A, B,C,De B a ; (2.256) 



That is, certain weighted averages of the surface gravities are guaranteed to be 
rational. If we wish to analyze whether rational ratios of R^ = ki/kj are implied 
in each of the particular cases of interest, we need to: 

a) Check if there exists some cu giving non-trivial subsets B a , leading to H2.246H . 



b) Analyze the sets of equations (|2.251[) implied by such an u . 



By proceeding in this way we are able to prove that periodicity of the QNFs im- 
plies rational ratios for the surface gravities in the following physically interesting 
cases: 



a) For j = 2m in equation (|2.177[) . 



b) For equation (|2.189|) when j ^ 2m + 1 and cos(7r/') ^ — ~. 



c) For equation (|2.203[) when j = 2m. 



2.7.1 Analysis of particular cases 

Now explore the familiar cases, which can serve also as particular examples 
described by this theorem. In the first part of each particular example we ex- 
plore whether periodicity implies rational ratios of the surface gravities (within 
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this particular case). If the rational ratios are implied, we can use the results of 
section 12.5.31 to determine the families and the gap structure. Despite having the 
analysis from section l2.5.3t we will derive (for each case) the gap structure also by 
analysing the equations (|2.251|) . This serves as: 

• a consistency check, 

• to bring more understanding in how the ideas used in the section |2.6.1| work. 

This derivation is made in the second part of each case analysis (to be exact, it is 
made for illustrative reasons for one set splitting only). 

Some common notation 

We will use the following notation: 

/* = £ — > ( Z257 ) 

i=i 

but if there is some index A, such that f^ = —Ja, we rename it to A = —A. So 
always f-A = — !a- We will also use 171a,b = mA — rriB and for the equation (|2.251[) 
of the form 

gap(w n )(/A - Ib) = 2m A ,B (2.258) 

we use the symbol Ea,b- Furthermore let us add one conceptual explanation: We 
say that equations E a ,b and E^ ^ are linearly independent if there do not exist 
such integers itia,b, m-A b> that t ne equations Ea,b and E^ ^ will be linearly depen- 
dent in the usual sense of the word. If Ea.b and E^ ^ are not linearly independent, 
we say they are linearly dependent. 

2 Horizons case, S-dS black hole by monodromy calculations: 
(1) spin 1 perturbations, (2) spin and 2 perturbations 

Take the first case, which is formula derived by the use of monodromy techniques 
for the S-dS black hole: 

cosh (™-™)+[ 1 + 2 oosfa)] cosh (™ + ™)= 0. (2.259) 
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Here j = for spin 0, 2 and j — 1 for spin 1. 
This becomes for spin 1 perturbation: 

/ TTLJ 7TU \ / TYUJ __ 7TUJ \ / TYUJ i 7TUJ \ / 7TUJ , 7TOJ \ 

e u - <V+e U - -+ j - e ^- + «+ j - e ^- + <V = 0, (2.260) 

and for spin and 2 perturbations: 

/ 1T(jJ 7VLJ \ / 7TUJ 7TUJ \ { 7TUJ i 7TUJ \ / 7Ttx; i TVLO \ 

e U_ K+ J +e U_ K+ ^ + 3e U_ + K+ ) + 3e U_+ K+ ) = . (2.261) 

In the case of spin 1 perturbation we already know that ratios of surface gravities 
might be completely arbitrary and we still get periodic solutions. 

The question of surface gravities rational ratios The matrix Za% has rank 2 (= 
H) hence if some subset of Da does not sum to 0, the rational ratio of surface 
gravities is implied. 

Let us have a look what happens here: The functions fi are given as 

and we have also /_i and /_ 2 . If the coefficients .D^ split into two sets each 
having two elements (which is the only way how they can be non-trivially split in 
this case), there are two equations and three possibilities how one can split them: 

• {£>i,£Li} and {D 2 ,D- 2 }, 

• {D U D 2 } and {D_ u D_ 2 }, 

• {£>!,£>_ 2 } and {D 2 , £>_!>. 

In terms of equations this leads to the following combinations: 

• Ex_x, E 2 ~2 gives 2 linearly independent equations, 
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• -E12, -E-1-2 gives only 1 linearly independent equation, 

• -Ei _2, #2,-i gives only 1 linearly independent equation. 

Because we have two surface gravities and two linearly independent equations 
the first combination of equations leads to the condition that the ratio of surface 
gravities must be rational. That means one needs to explore only the second and 
the third combination of equations. Here we have to check the step a) from the end 
of the previous section. If there exists u giving us B a sets leading to the second, 
or the third combination of equations, it must fulfill the following conditions: 

• In the case of the second combination (#1,2, #1, 2) it must fulfil the equa- 
tions 



TTWo 
K + 

TUJO 



i2Tcmi 2 + In 



C'i 



Co 



i27rm_i _2 + In 



C_i 



C_, 



(2.262) 
(2.263) 



In the case of the third combination (Ex _ 2 , #-1,2) it must fulfil the equations 



nu 



TCUo 



K. 



il-nvrix -2 + In 



Cx 



i27rm_x 2 + In 



C_ 2 



C, 



But for the second combination this means: 



In 



Cx 



C, 



-In 



C-x 



and for the third combination this means: 



In 



Cx 



C_ 



-In 



C_ 



C-x 



Co 



From (|2.260|) we see that for spin 1 perturbation it holds 

\n I \n I _ \n 1 _ \n 1 _ 1 

w — U2 — o_ x — o_ 2 — 1. 



(2.264) 
(2.265) 

(2.266) 

(2.267) 

(2.268) 
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As a result of (|2.268[) the conditions (|2.266[) and (|2.267[) are trivially fulfilled, since 
all the logarithms are 0. But for spin and 2 \C\\ = |C_i| = 1 and |C 2 | = |C_ 2 | = 
3, which means that in each case we get ln(|). This means there is no way how to 
fulfil (|2.266[) , or (|2.267|) , hence split the coefficients in such way that we do not get 



the surface gravities rational ratio condition. This means that for spin and spin 2 
periodicity implies the rational ratio of surface gravities. 

For the spin 1 perturbation we already showed we can find explicit uiq-s lead- 
ing to two different families of QNMs, each family related to different horizon 
surface gravity (inn±). This means that in the case of spin 1 perturbation rational 
ratios are not implied. Note also that each of the two families is just a result of 
different splitting of coefficients Da- 

"Gap" derivation by using our approach (spin 1 perturbation) As previously 
noted, the general solutions split into two families, one related to one surface grav- 
ity (to = innJ), the other to another (cu = inn + ). Now take one solution from the 
set {inn-} (for example u = IkJ) and substitute it in (|2.260[) to get the coefficients 
D A . We see that the two set splitting of D A is: 



D x = -e K + , D_ 2 = e K + V, 1 D_! = -e »+ , D 2 = e «+ i- . (2.269) 

(We can easily observe that the sum of those couples of coefficients within each 
set is 0). Now, since E\ _ 2 and ELi i2 are linearly dependent, the second equation 
can be taken only as a definition of m_i j2 . The only independent equation is then: 

gap(w„) = mi_ 2 «— 

To obtain from this equation the basic gap structure, we have to choose the integer 
rxi\ _ 2 to be such, that we obtain the "narrowest" gap. This is obtained by mi _ 2 = 
1, confirming the known result. 

If we start with the second set of solutions related to the cosmological horizon 
surface gravity, we can proceed completely analogically and verify also the second 
result. The gap will be in such case given by the cosmological horizon surface 
gravity. 



CHAPTER 2. ANALYTIC RESULTS FOR HIGHLY DAMPED QNMS 92 

"Gap" derivation by using our approach (spin and 2 perturbation) Here we 
proved we have 2 linearly independent equations (the rest of the equations are 
only definitions of itia,b terms): 

• gap(w n ) = -mi >2 K+, 

• gap(w n ) = -mi _ 2 «— 

Now, they immediately lead to the condition — = mi ' 2 , hence surface gravities 
ratio being rational. 

To explore the gap structure one has to explore all the other independent con- 
ditions as well (here there are three independent conditions obtained by equating 
4 terms). This is because definitions of itia,b might put some additional constrains 
on the gap function. (The constraints come from the fact that the vtia,b have to 
always be integers. This will be seen in the analysis of the following case.) But in 
this case the third independent condition can be given as 

gap(w n ) ( J = mi_i (2.270) 

\k_ k + J 

and is consistent with taking the "narrowest" gap as 

/ x P- P+ 1 

gap(w„) = k* = K_mi,_ 2 = K-P- = n + m 12 = K + p + = — = — = —. 

(Here 6* is the highest common divisor of 6_ , b + with respect to integers.) This 
means: —7711,2 = P+ arid — 7771,-2 = P— 

2 Horizons, S-dS black hole by analytically solvable potentials 

The equation we obtained is the following: 

( 1XUJ lXL>j\ ( 1XUJ lXL>j\ . . . . /~~r-M\ 

cosh — cosh 1 = 2cos(7ro; + ) cos(7T«_). (2.271) 

\K- K + J \K- K + J 
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The question of surface gravities rational ratios In the previous section we 
already proved that in this case the rational ratios are implied by the periodicity 
We can prove the same by analyzing all the possibilities how to split coefficients 
into different sets {D A }. It is an alternative proof to the one in the previous 
section. It is definitely a more complicated proof, but its advantage is that it is 
part of more general approach, which generates proofs for all the cases where the 
implication holds. 

Let us show the "ugly" way of proving it by the new general method: The 
matrix Za% is of rank 2 (= H) and hence if the Da coefficients do not split into 
nontrivial subsets, the rational ratios of surface gravities is proven. Now let us 
explore what ways of splitting the coefficients one can obtain and what they gen- 



erally mean. We can rewrite the equation (|2.271[) to get the following functions: 



A - i: + i' / 2 - it + i' h- 3 (;~ + i:)> 



ft ~ K_ + K + ' f 5 ~ 2 l K_ K+ j • 



Let us explore what happens when we split the D A coefficients into two sets, 
one having three and the other two elements (these are the only possible nontriv- 
ial ways of splitting the coefficients). This gives for each splitting 3 equations. 
One can analyze all the combinations of equations in the following way: Pick f ^ 
belonging to a coefficient in the three element set and take all combinations of 
such equations Ea,b (related to all the possible ways of splitting the coefficients 
in which D ^ is in the set of three elements), that they do not involve the given f A . 
If these combinations give 2 independent equations we are finished, and do not 
need to explore the 3-rd equation involving fj. If there is only one independent 
equation from the given couple, we need to explore if there is a way how to add an 
equation involving fa (E^ B ), by keeping the number of independent equations 
to be still 1. The I/D letters in the table mean linearly independent/ dependent: 
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h 


combination 1 


I/D 


combination 2 


I/D 


combination 3 


I/D 


h 


-^2,3) -^4,5 




^2,4j -^3,5 


I 


^2,5) ^3,4 


I 


h 


-El,3, -^4,5 




-El,4, -^3,5 


I 


-^1, 5, -^3,4 


I 


h 


-El, 2, -^4,5 




-El,4, £-2,5 


I 


-El, 5, -£"2,4 


I 


h 


-El, 2, -£"3,5 




-El,3, -^2,5 


I 


-El, 5, -£"2,3 


I 


h 


-^1, 2) -^3^4 




-El,3> -^2,4 


D 


-El,4, -^2,3 


D 



We see that only the last two cases give linearly dependent equations, but it 
is very easy to check that in both of these cases it holds, that if we add arbitrary 
third equation of the type E 5tB (hence involving / 5 ), the number of linearly inde- 
pendent equations grows to 2. This is a result of the fact that in _E 13 , E 2 ^ there is 
only k + present and in £^4, £2,3 there is only k present. This proves that here 
the rational ratios are implied by the periodicity. 



The "gap" derivation by using our approach If there are no non-trivial ways 
of splitting the coefficients, then we get from our result (|2.251|) 4 independent 
conditions (although, in the sense defined, only 2 linearly independent equations) 
and they are: 

1) gap(w n ) (±- - ±) = mi, 2 , 

2) gap(w„) = -mi )3 «+ , 

3) gap(w n ) = mi,4«- , 

4) gap(w n ) (± - i) = 2mi, 5 • 

Now any two from the first three equations^ give surface gravity rational ratio 
condition. Choose for this purpose, for example, equations 2) and 3). (Notice here 
that in order to get the surface gravities ratio positive, my, 777,1,4 must have oppo- 
site signs.) Equation 1) is then uninteresting since it is only defining rri\ $, without 
constraining the possible gap function. Now the whole problem is encoded in 



5 The same holds for any two of the last three equations. 
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the equation 4) and particularly in the fact that there is 2 on the right side of the 
equation. In fact it leads to m^ 4 + m^ 3 = 2m lj5 , and that means one can fulfil the 
equation 4) only if both 7711,4, 7711,3 are odd or both are even. This means that the 
definition of m lj5 constrains the gap function. Now realize that only the ratio of 
m i,3) m i,4 is determined by the ratio of the surface gravities p_/p + , so there is still 
freedom to multiply both, p + and p-, by the same arbitrary scaling integer. We 
have to determine the integer to fulfill the equation 4) and simultaneously to give 
the "narrowest" gap. But then, if the ratio 

«- = K = P± 

is in its most reduced form given by one odd and one even number (p + ■ p is 
even), we have to take both mi j4 , 7711,3 even by: mi j4 = 2p_ and —7711,3 = 2p+. Then 
the gap will be gap(w n ) = 2/t* = £- (k* = n_p^ = K + p + — £- ) . In the case both 
P-,p+ are odd numbers (hence p_ • p + is odd), the 4)-th equation is automatically 
fulfilled and to get the "narrowest" gap one chooses the scaling integer to be ±1, 
hence 777.1,4 = p_ and 7771,3 = ~P+- Then the gap is obtained as gap(u n ) = k* 






2 Horizons, R-N black hole 

Take as another example the equation (|2.189|) with j = 0. We can rewrite it as: 

e™ ( ^ + - ) + 3e-™ ( ^ + - ) + 3e^ ( - ~^ ] = 0. (2.272) 

There cannot be any non- trivial subset of coefficients Da summing up to 0. (This is 
because there are three non-zero coefficients, so it cannot happen that two of them 
will give by the summation). That means we have always two independent 
conditions: 

1) gap(w n ) = m lj3 K + , 

2) gap(w n ) = -777 2 ,3K- • 

But 1) and 2) are also linearly independent equations, so this immediately implies 
the same results, as in the previous case. These are: 
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• the rational ratio of surface gravities given as -m li3 /m 2i 3 (note again that to 
get the ratio positive, m 1)3 m 2j3 have to be chosen with the opposite signs), 

• the "narrowest" gap given by m 13 = p + and m 2 , 3 = —p~, hence 

gap(c; n ) = K* = n_p_ = K + p + = — . 

So in this case the rational ratios are implied by the periodicity. 
3 Horizons, Monodromy results: R-N-dS black hole 



Take the equation (|2.203[) with any j, for which the coefficients are non-zero (trivial 



example is standard j = ). 

The question of surface gravities rational ratios Now this is a case with six 
functions Ja'- 

h = — -—, h = — + — , h = — + — ~— and /_!, /_ 2 , /_ 3 . 

Here we have 3 surface gravities and the rank of Zai matrix is 3, so if there is no 
nontrivial subset of Da coefficients summing to 0, the surface gravities rational 
ratios are implied. 

Let us explore the nontrivial ways of splitting the coefficients. The splitting of 
coefficients giving the lowest number of independent conditions is the splitting 
into three sets, each having 2 elements. It gives 3 equations (and we need three 
linearly dependent). There are three basic ways how to split the coefficients: 

• The first splitting is giving E lt _i, E 2 -2, E 3 _ 3 , but these are necessarily 3 
linearly independent equations as /i, / 2 , fz are linearly independent set of 
functions. 

• The other type of splitting is E AlA2 , E_ A2 ,a^ e -m-a x , A = ±1,±2,±3, 
\Ai\ ^ \Aj\ if i 7^ j. But all these ways of splitting the coefficients give 3 
linearly independent equations as well. 
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• The third basic way how to split coefficients into three two element sets gives 
E Al _ Al , E A2i a 3 i E_ A2 _ A[i , Ai being defined as before. This splitting gives 
only 2 linearly independent equations. 

That means in the last case we have to take the second step. Fortunately here we 
can use for different C A coefficients multiplying different cosh(-) terms the same 
argument as in the case of S-dS and spin and 2 perturbations. This argument 
excludes the possibility of such ways of splitting the coefficients. 

Any splitting into two sets, one having 2 and the other 4 elements is just more 
constrained version of some "three set each having two elements" splitting. This 
means the arguments provided in the previous case transfer automatically to this 
case. 

The last, deeply nontrivial splitting is when D A coefficients split into two sets, 
each having three elements. In such case we have 4 independent conditions. There 
are two possible basic ways of splitting the equations: 

• E Al _ Al , E AuA2 , E_ A2tA3 , E A3 _ A3 

• E AuA2 , E A2tA . A , E- Al ,-A 2 , E- A2 -a 3 

Now the first splitting leads to three independent equations, but the second split- 
ting only to two independent equations. So in the second splitting case one has 
to proceed to the next step: One can observe that if there exists an cu giving the 
second splitting, then the following two equations have to be fulfilled. Take: 



K + KC' 




^±(A + £> 




5 = ± ( -2- + -i- + ±) 

\K- K + KC J 


and 


K = 1 + cos(ttj). 
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Then the following holds: 

z a + Kz p + 2Kz 5 = 0, (2.273) 

z~ a + Kz-P + 2Kz~ s = 0. (2.274) 

K 

z a +KzP 



But then by substituting to (|2.274[) z 5 = — za + Kz p we will obtain 



K 
and then we get the result: 



1 3h ' 2 + z-f + z~("-» 



2K 2]j \ K J 

For K > 4= we have 2; Q_/3 > 0. We see that if we take j = \, the only nontrivial 
cases ( K ^ ) are i^ = 1, 2 and fulfil this condition. But the equation (|2.273[) is 
in fact: 

Z *-P + K + 2Kz {a ~ p) ^ = 0. (2.275) 



Define by u = z a P. Then by solving (|2.275|) we obtain 

5-/3 , / u + K 



a - p ou \ 2K 
But since the logarithm argument is for K > 4= negative (and w is in such case 
positive), it means that ^5§ rnust be nonreal. But considering how a, (3,5 are de- 
fined, this cannot happen for surface gravities being real numbers. Hence it is 
proven that there is no ojq giving the splitting considered. This means we proved 
that for (|2.203[) the periodicity implies rational ratios of surface gravities. 



The "gap" analysis by using our approach If we consider again only such u -s, 
for which there is no nontrivial splitting of Da coefficients, we obtain the follow- 
ing 5 independent conditions: 



gap(w n )-J- = 7712,3/ 
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• gap(w n )^ =1711-2, 

• gap(o; n )^ = -rai,2, 

• g ap(u; n )(^ : -i) = mi,_i # 

• gap(w n )^- = m_ 2 ,-3- 

Now note: the last two equations are uninteresting, since they are only definitions 
of mi_i and m_ 2 ,-3 without putting any constraint on the gap function. On the 
other hand the first three equations tell us that the ratio between arbitrary two 
surface gravities is a rational number. Then to get the "narrowest" gap, consider 

the following: 

K- _ p+_ __ fe+ k_ _ pc_ _ bc_ 

k + pi_ &_' k c Vi- b_' 

Further k* = ■£- — K_pi_ = n + p + and «£ — y- — K-P2- = ^cVc- Then the gap 
must be given by gap(w n ) = 1chi{k*, k'^} = I x k* = Z 2 «C (h, h G ^)- That means the 
m integers must be taken as: m 2>3 = h ■ pi_ , mi -2 — h • P+ and — my = Z2 • Pc- 

2.7.2 Conclusions 

It seems to be very hard (if not impossible) to find a more constrained form of the 
monodromy results than the formulation (|2.169|) . Since in (|2.169|) it is not a general 



result (but seems to be generic enough) that the rational ratios of surface gravities 
are implied by periodicity (see S-dS spin 1 perturbations as counter-example), we 
developed a general method how to prove the implication for every case in which 
it holds. The necessary step was to prove the theorem at the beginning of the sec- 
ond section. By using our method we proved that the implication (periodicity — > 
rational ratios) holds in every monodromy case (from the cases described before), 
apart of the one given counter-example. 



Chapter 3 

Multiplication of tensorial 
distributions 



3.1 Introduction 

This chapter is devoted to a topic from the field of mathematical physics which is 
closely related to the general theory of relativity It offers possible significant con- 
ceptual extensions of general relativity at short distances /high energies, and gives 
another arena in which one can conceptually /physically modify the classical the- 
ory. But the possible meaning of these ideas is much wider than just the general 
theory of relativity It is related to the general questions of how one defines the 
theory of distributions for any geometrically formulated physics describing interac- 
tions. 

The main reasons why we "bother" Let us start by giving the basic reasons 
why one should work with the language of distributions rather than with the old 
language of functions: 

• First there are deep physical reasons for working with distributions rather 
than with smooth tensor fields. We think distributions are more than just a 
convenient tool for doing computations in those cases in which one cannot 

100 
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use standard differential geometry. We consider them to be mathematical 
objects which much more accurately express what one actually measures in 
physics experiments, more so than when we compare them to the old lan- 
guage of smooth functions. The reason is that the question: "What is the 
'amount' of physical quantity contained in an open set?" is in our view 
a much more reasonable physical question, (reasonable from the point of 
view of what we can ask the experimentalists to measure), than the ques- 
tion: "What is the value of quantity at a given point?". But "point values" 
as "recovered" by delta distributions do seem to give a precise and reason- 
able meaning to the last question. There is also a strong intuition that the 
"amount" of a physical quantity in the open set Cli U f2 2 , where fii, f2 2 are 
disjoint is the sum of the "amounts" of that quantity in fii and fi 2 . That 
means it is more appropriate to speak about distributions rather than gen- 
eral smooth mappings from functions to the real numbers (the mappings 
should also be linear). 

• The second reason is that many physical applications suggest the need for 
a much richer language than the language of smooth tensor fields. Actually 
when we look for physically interesting solutions it might be always a matter 
of importance to have a much larger class of objects available than the class 
of smooth tensor fields. 

• The third reason (which is a bit more speculative) is the relation of the lan- 
guage defining the multiplication of distributions to quantum field theories 
(but specifically to quantum gravity). Note that the problems requireing 
distributions, that means problems going beyond the language of classical 
differential geometry, might be related to physics on small scales. At the 
same time understanding some operations with distributions, specifically 
their product has a large formal impact on the foundations of quantum field 
theory, particularly on the problem with interacting fields. (See, for example, 
||73| .) As a result of this it can have significant consequences for quantum 
gravity as well. 
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The intuition behind our ideas Considerations about language intuitiveness 
lead us to an interesting conclusion: the language of distributions (being con- 
nected with our intuition) strongly suggests that the properties of classical tensor 
fields should not depend on the sets having (in every chart) Lebesgue measure 
0. If we follow the idea that these measure zero sets do not have any impact on 
physics, we should naturally expect that we will be able to generalize our lan- 
guage from a smooth manifold into a piecewise smooth manifold (which will 
bring higher symmetry to our conceptual network). The first traces of piecewise 
smooth coordinate transformations are already seen, for example, in [100]. 

The current situation is "strange" Now it is worth noting how strange the cur- 
rent situation of the theory of distributions is: we have a useful and meaningful 
language of distributions, which can be geometrically generalized, but this lan- 
guage works only for linear physics. But linear physics is only a starting point (or 
at best a rough approximation) for describing real physical interactions, and hence 
nonlinear physics. So one naturally expects that the "physical" language of dis- 
tributions will be a result of some mathematical language defining their product. 
Moreover, at the same time we want this language to contain the old language of 
differential geometry (as a special case), as in the case of linear theories. It is quite 
obvious that the nonlinear generalization of the geometric distributional theory 
and the construction of generalized differential geometry are just two routes to 
the same mathematical theory. The natural feeling that such theory might exist is 
the main motivation for this work. The practical need of this language is obvious 
as well, as we see in the numerous applications ||TTll6Tl[73l[76l [T04l[T56l[T58l[T73l . 
(Although this is not the main motivation of our work.) But is it necessary that 
such more general mathematical language exists as a full well defined theory? 
No, not at all. The potentially successful uses of distributions that go behind the 
Schwarz original theory might be only "ad hoc" from the fundamental reasons. 
Take the classical physics. The success of such uses of distributions here might 
not be a result of some more general language than smooth tensor calculus being 
a classical limit of the more fundamental physics. There might be only hidden 
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specific reasons in the more fundamental physics why such "ad hoc" calculations 
in some of the particular cases work. But it is certainly very interesting and im- 
portant to explore and answer the question whether: (i) such a full mathematical 
langauge exists, and (ii) to which extent it is useful to the physics community 
For the first question this work suggests a positive answer, to answer the second 
question much more work has to be done. 

So what did we achieve? The main motivation for this work is the develop- 
ment of a language of distributional tensors, strongly connected with physical in- 
tuition. (This also means it has to be based on the concept of a piecewise smooth 
manifold.) This language must contain all the basic tensorial operations in a gen- 
eralized way, enabling us to understand the results the community has already 
achieved, and also the problems attached to them. It is worth to stress that some 
of this motivation results from a shift in view regarding the foundations of clas- 
sical physics (so it is given by "deeper" philosophical reasons), but it can have 
also an impact on practical physical questions. The scale of this impact has still 
to be explored. We claim that the goals defined by our motivation (as described 
at the beginning) are achieved in this work. Particularly, we have generalized all 
the basic concepts from smooth tensor field calculus (including the fundamental 
concept of the covariant derivative) in two basic directions: 

• The first generalization goes in the direction of the class of objects that, (in 
every chart on the piecewise smooth manifold), are indirectly related to sets 
of piecewise continuous functions. This class we call D'^ EA {M). For a de- 
tailed understanding see section 13.4.31 

• The second generalization is a generalization to the class of objects naturally 
connected with a smooth manifold belonging to our piecewise smooth man- 
ifold (in the sense that the smooth atlas of the smooth manifold is a subatlas 
of our piecewise smooth atlas). This is a good analogy to the generalization 
known from the classical distribution theory. The class of such objects we 
call D'J^ So) (M). For detailed understanding see again the section 13.4.31 



CHAPTER 3. MULTIPLICATION OF TENSORIAL DISTRIBUTIONS 104 

We view our calculus as the most natural and straightforward construction achiev- 
ing these two particular goals. The fact that such a natural construction seems to 
exist supports our faith in the practical meaning of the mathematical language 
here developed. 

The last goal of this chapter is to suggest much more ambitious, natural gen- 
eralizations, which are unfortunately at present only in the form of conjectures. 
Later in the text we provide the reader with such conjectures. 

The structure of this chapter The structure of this chapter is the following: In 
the first part we want to present the current state of the Colombeau algebra the- 
ory and its geometric formulations. We want to indicate where its weaknesses are. 
This part is followed by several technical sections, in which we define our theory 
and prove the basic theorems. First we define the basic concepts underlying our 
theory. After that we define the concept of generalized tensor fields, their impor- 
tant subclasses and basic operations on the generalized tensor fields (like tensor 
product). This is followed by the definition of the basic concept of our theory, the 
concept of equivalence between two generalized tensor fields. The last technical 
part deals with the definition of the covariant derivative operator and formulation 
of the initial value problem in our theory. All these technical parts are followed 
by explanatory sections, where we discuss our results and show how our theory 
relates to the practical results already achieved (as described in the first part of the 
chapter). 

3.2 Overview of the present state of the theory 

3.2.1 The theory of Schwartz distributions 

Around the middle of the 20-th century Laurent Schwartz found a mathemati- 
cally rigorous way for extending the language of physics from the language of 
smooth functions into the language of distributions. Physicists such as Heavi- 
side and Dirac had already given good physics reasons for believing that such a 
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mathematical structure might exist. 

The classical formulations of the distribution theory were directly connected 
with MJ 1 and were non-geometric. Distributions in such a formulation are typically 
understood as continuous, linear maps from compactly supported smooth func- 
tions (on M n ) to real numbers. (The class of such functions is typically denoted 
by D(R. n ). The dual to such space, which is what distributions are, is typically 
denoted by D'iW 1 ).) Here the word "continuous" refers to the following topol- 
ogy on the given space of compactly supported smooth functions: The sequence 
of smooth compactly supported functions fi(xi) converges for I — > oo to a smooth com- 
pactly supported function f(xi) iff an arbitrary degree derivative with respect to arbitrary 
variables 1/' 7™ k [Yli=i m k — n) converges uniformly on each compact W 1 subset to 

d n f(x t ) 



o rrn q 77r k • 

ox 1 ■■■ox k 

The alternative classical way of formulating the theory of distributions is to 
extend the space of test objects to be such that they still follow appropriate fall-off 
properties. The properties can be summarized as: 

f(xi) belongs to such space if it is smooth and 

d a f{x l 
xu.TxlUoo V~ x "'~ fe dx™\..dxZ 



(V«, V/3) __ lim __ ( x? . -.xT ^ mi J \% ) = 0, (3.1) 



k k 

t=i i=i 

This is a topological space with topology given by a set of semi-norms P a> p de- 
finecu as 



sup 



ni n k & J \ X i 

I ■ ■ ■ J^l 



dx 1 1 ...dx k 

k k 



(3.2) 



rrij = a. 



Naturally, the space of continuous and linear maps on such a space is more re- 
stricted as in the first, more common formulation. Objects belonging to such duals 



1 We admit that the notation P a .p, might be somewhat misleading, since the a, j3 values do not 
specify the semi-norm in a unique way. 
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are called "tempered distributions". The advantage of this more restricted version 
is that Fourier transform is a well defined mapping on the space of tempered dis- 
tributions (see [122]). 



If we refer to the more common, first formulation, the space of distributions 
accommodates the linear space of smooth functions by the mapping: 

/(zi)-> / f( Xi )---d n x. (3.3) 

Here f(xi) is a smooth function on M n and J Rn f(x) ■ ■ ■ dx is a distribution de- 
fined as the mapping: 

*(xi) -> f f{xi)*{xi) d n x, *{xi) e D(R n ). (3.4) 

Moreover, by use of this mapping one can map into the space of distributions 
any Lebesgue integrable function (injectively up to a function the absolute value 
of which has Lebesgue integral 0) . The distributional objects defined by the im- 
ages of the map (|3.3[) are called regular distributions. The map (|3.3|) always pre- 



serves the linear structure, hence the space of Lebesgue integrable functions is a 
linear subspace of the space of distributions. But the space of distributions is a 
larger space than the space of regular distributions. This can be easily demon- 
strated by defining the delta distribution 

6(V) = #(0) (3.5) 

and showing that such mapping cannot be obtained by a regular distribution. 
Note particularly that delta distribution had an immediate use in physics in the 
description of point-like sources. (Its intuitive use in the work of Paul A.M. Dirac 
before the theory of Schwartz distributions was found, was one of the main physics 
reasons why people searched for such language extension.) 

These considerations show that the space of distributions is significantly larger 
than the space of smooth functions. The space of distributions is classically taken 
to be a topological space with the weak (<r— ) topology. In this topology the space 
of regular distributions given by smooth functions is a dense set. 
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Moreover one can continuously extend the derivative operator from the space 
of C 1 (M n ) functions to the space of distributions by using the definition: 

T x i(^) = -T(^ jXi ), (3.6) 

(where T denotes a distribution). This means that C°°(IR n ) functions form not 
only a linear subspace in the space of distributions, but also a differential linear 
subspace. It looks like there stands "almost" nothing in the way of fully and satis- 
factorily extending the language of C 1 (W n ) functions to the language of Schwartz 
distributions. Unfortunately there is still one remaining trivial operation and this 
is the operation of multiplication. That means we need to obtain some distribu- 
tional algebra having as subalgebra the algebra of Lebesgue integrable functions 
factorized by functions the absolute value of which has Lebesgue interal 0. Unfor- 
tunately, shortly after Laurent Schwartz formulated the theory of distributions he 
proved the following "no-go" result [138J: 

The requirement of constructing an algebra that fulfills the following three con- 
ditions is inconsistent: 

a) the space of distributions is linearly embedded into the algebra, 

b) there exists a linear derivative operator, which fulfills the Leibniz rule and 
reduces on the space of distributions to the distributional derivative, 

c) there exists a natural number k such, that our algebra has as subalgebra the 
algebra of C k (W n ) functions. 

This is called the Schwartz impossibility theorem. There is a nice example show- 
ing where the problem is hidden: Take the Heaviside distribution H. Suppose 
that a) and b) hold and we multiply the functions /distributions in the usual way. 
Then since H m = H the following must hold: 

H' = (H m )' = 5 = mH m ~ 1 5 . (3.7) 

But this actually implies that 6 = 0, which is nonsense. 
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The closest one can get to fulfill the conditions a) — c) from the Schwarz impos- 
sibility theorem is the Colombeau algebra (as defined in the next section) where 
conditions a) — c) are fulfilled with the exception that in the c) condition k is taken 
to be infinite. This means that only the smooth functions form a subalgebra of 
the Colombeau algebra. This is obviously not satisfactory, since we know (and 
need to recover) rules for multiplying multiply much larger classes of functions 
than only smooth functions. In the case of Colombeau algebras this problem is 
"resolved" by the equivalence relation, as we will see in the following section. 

3.2.2 The standard W 1 theory of Colombeau algebras 
The special Colombeau algebra and the embedding of distributions 

The so called special Colombeau algebra is on R n defined as: 

G(R n ) = MM n )/A/"(R n ). (3-8) 

Here £a/(R™) (moderate functions) is defined as the algebra of functions: 

W 1 x (0, 1] -> R (3.9) 

that are smooth on R n (this is usually called £(R n )), and for any compact subset 
K of R™ (for which we will henceforth use the notation K CC R") it holds that: 



& 



Va G N™, 3p G N suchthafl sup \D a f e ( Xi )\ < 0(e~ p ) as e ->■ 0. (3.10) 

The A/"(R n ) (negligible functions) are functions from £(R n ) where for any K CC 
R n it holds that: 

Va G N£, VpeNwe have sup \D a f e ( Xi )\ < 0(e p ) as e ->■ 0. (3.11) 

The first definition tells us that moderate functions are those whose partial deriva- 
tives of arbitrary degree (with respect to variables Xi) do not diverge faster then 



2 In this definition, (and also later in the text), the symbol Nq denotes a sequence of n members 
formed of natural numbers (with included). 
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any arbitrary negative power of e, as e — v 0. Negligible functions are those mod- 
erate functions whose partial derivatives of arbitrary degree go to zero faster than 
any positive power of e, as e —> 0. This simple formulation can be straightfor- 
wardly generalized into general manifolds just by substituting the concept of Lie 
derivative for the "naive" derivative used before. 

It can be shown, by using convolution with an arbitrary smoothing kernel (or 
mollifier), that we can embed a distribution into the Colombeau algebra. By a 
smoothing kernel we mean, in the widest sense a compactly supported, smooth 
function p 6 (xj), with e G (0, 1], such that: 

• supp(p e ) -> {0} for (e -> 0), 

• / Rn p e (xi) d n x -»■ 1 for (e -)> 0), 

• V77 > 3tf, Ve G (0, V ) sup Xi \p e ( Xi )\ < C. 

This most generic embedding approach is mentioned for example in J100| (in 
some sense also in ||104ip . 

More "restricted" embeddings to QiW 1 ) are also commonly used. We can 
choose for instance a subclass of mollifiers called A (M. n ), which are smooth func- 
tions from /}(IR n ) (smooth, compactly supported) and (i.e. [72]) such that 

Ve holds I p e (x l )d n x = 1. (3.12) 

Their dependence on e is giverjj as 

Pe ( Xl ) = -p ( X A . (3.13) 



Sometimes the class is even more restricted. To obtain such a formulation, we 
shall define classes A m (M n ) as classes of smooth, compactly supported functions, 
such that 

/ x\-" x[ 4>(xi)d n x = 5 0k for i-\ \- j = k < m. (3.14) 



3 Later in the text will the notation p e (xi),ip e {%i) (etc.) automatically mean dependence on the 
variable e as in l|3.13t . 



CHAPTER 3. MULTIPLICATION OF TENSORIAL DISTRIBUTIONS 110 

Clearly A m+1 (W l ) c A m (W l ). Then the most restricted class of mollifiers is taken 
to be the class 74°°(IR n ). This approach is taken in the references ||76l 11241 [1561 158, 
IT731 . 



Even in the case of the more restricted class of mollifiers the embeddings are 
generally non-canonical [721I173L The exception are smooth distributions, where 



the difference between two embeddings related to two different mollifiers is al- 
ways a negligible function. 

The full Colombeau algebra and the embedding of distributions 

What is usually considered to be the canonical formulation of Colombeau algebras 
in R n is the following: The theory is formulated in terms of functions 

R n x ^0( R n) _^ R ( caU them F j _ (3 ^ 

The Colombeau algebra is defined in such way that it is a factor algebra of mod- 
erate functions over negligible functions, where: 

• Moderate functions are functions from F that satisfy: 

Vm G NJ, VK CC R n 3N G N such that if e A N (R n ) } there are a, p > 0, 

such that sup \D m F(0 e ,Xi)\ < ae~ N if < e < p. (3.16) 

Xi£K 

• Negligible functions are functions from F that satisfy: 

Vm e N, VK CC R n , Vp G N 3g G N such that if e A q (W n ), 3a, p > 0, 

we have sup \D m F{(j) € ,Xi)\ < ae p if < e < p. (3.17) 

XiEK 

Then ordinary distributions automatically define such functions by the convo- 
lution (B6H5DHlQl|TZa etc.)C 



B -»■ B v 



1 /j/j -Xi 



c v 



G A u (R n ). (3.18) 



4 Here the "B x " notation means that x is the variable removed by applying the distribution. 
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Important common feature of both formulations 

All of these formulations have two important consequences. Given that C denotes 
the embedding mapping: 

• Smooth functions (C°°(IR' 1 )) form a subalgebra of the Colombeau algebra 
(C(f)C(g) = C(f ■ g) for /, g being smooth distributions). 

• Distributions form a differential linear subspace of Colombeau algebra (this 
means for instance that C(f') = C'(f) ). 

The relation of equivalence in the special Colombeau algebra 

We can formulate a relation of equivalence between an element of the special 
Colombeau algebra f e (xi) and a distribution T. We call them equivalent, if for 
any e D(W n ), we have 

lim / f e (x i )<P(x i )d n x = T(<P). (3.19) 

Then two elements of Colombeau algebra f e (xi), g e {xi) are equivalent, if for any 
(f>(xi) e D(R n ) 

lim f (f t ( Xl ) - g e ( Xi )) </>(xi) d n x = 0. (3.20) 



For the choice of A 00 (M") mollifiers the following relations are respected by the 
equivalence: 

• It respects multiplication of distribution by a smooth distribution 111731 in 
the sense that: C(f ■ g) « C(f) ■ C(g), where / is a smooth distribution and 

g e D'(W n ). 

• It respects (in the same sense) multiplication of piecewise continuous func- 
tions (we mean here regular distributions given by piecewise continuous 
functions) H46l . 

• If g is a distribution and / ~ g, then for arbitrary natural number n it holds 
D n f^D n g EDI. 
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• If / is equivalent to distribution g, and if h is a smooth distribution, then / • h 
is equivalent to g ■ h ||50| . 



The relation of equivalence in the full Colombeau algebra 

In the canonical formulation the equivalence relation is again formulated either 
between an element of the Colombeau algebra and a distribution, or analogously 
between two elements of the Colombeau algebra: If there 3m, such that for any 

G A m (R n ), and for any *(xj) G D(R n ), it holds that 

lim / (/(&, Xi ) - <?(&, Xi))*(xi) d n x = 0, (3.21) 

then we say that / and g are equivalent (/ ~ g). For the canonical embedding and 
differentiation we have the same commutation relations as in the non-canonical 
case. It can be also proven that for /i . . . /„ being regular distributions given by 
piecewise continuous functions it follows that 

C(A)...C(/ n )«C(A-/n), (3.22) 

and for / being arbitrary distribution and g smooth distribution it holds that 

C(f)-C(g)nC(f-g). (3.23) 

How this relates to some older Colombeau papers 

In older Colombeau papers Ii46l |47| all these concepts are formulated (equiva- 
lently) as the relations between elements of a Colombeau algebra taken as a sub- 
algebra of the C°°(-D(IR n )) algebra. The definitions of moderate and negligible 
elements are almost exactly the same as in the canonical formulation, the only 
difference is that their domain is taken here to be the class D(W l ) x M n (being 
a larger domain than 74°(IR n ) x W 1 ). The canonical formulation is related to the 
elements of the class C°°(D(W 1 )) through their convolution with the objects from 
the class D(R n ). It is easy to see that you can formulate all the previous relations 
as relations between elements of the C oc (D(W n )) subalgebra (with pointwise mul- 
tiplication), containing also distributions. 



CHAPTER 3. MULTIPLICATION OF TENSORIAL DISTRIBUTIONS 113 

3.2.3 Distributions in the geometric approach 

This part is devoted to review the distributional theory in the geometric frame- 
work. How to define arbitrary rank tensorial distribution on arbitrary manifolds 
by avoiding reference to preferred charts? Usually we mean by a distribution rep- 
resenting an (m, n) tensor field an element from the dual to the space of objects 
given by the tensor product of (m, n) tensor fields and smooth compactly sup- 
ported /c-form fields (on k dimensional space). That means for example a regular 
distributional (m, n) tensor field Bfr" is introduced as a map 

T®u^ J B^ a T^u. (3.24) 

This is very much the same as to say that the test space are smooth compactly 
supported tensor densities T„""^ [60. 61J. The topology taken on this space is the 
usual topology of uniform convergence for arbitrary derivatives related to arbi- 
trary charts (so the convergence from W 1 theory should be valid in all charts). The 
derivative operator acting on this space is typically Lie derivative. (Lie derivative 
along a smooth vector field £ we denote Lg.) This does make sense, since: 

• To use derivatives of distributions we automatically need derivatives along 
vector fields. 

• Lie derivative preserves p-forms. 

• In case of Lie derivatives, we do not need to apply any additional geometric 
structure (such as connection in the case of covariant derivative). 



There is an equivalent formulation to |6i~|, given by [173]], which takes the 



space of tensorial distributions to be D'(M) ® T™(M). Here D'(M) is the dual 
to the space of smooth, compactly supported fc-form fields (k dimensional space). 
Or in other words, it is the space of sections with distributional coefficients. In 
111011 the authors generalize the whole construction by taking the tensorial distri- 
butions to be the dual to the space of compactly supported sections of the bundle 
E* ® Vb/ 1-9 . Here Vol 1 ~ q is a space of (1 — g)-densities and E* is a dual to a 
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tensor bundle E (hence to the dual belong objects given as E <g> Vol q ). In all these 
formulations Lie derivative along a smooth vector field represents the differential 
operatoiB 

Let us mention here that there is one unsatisfactory feature of these construc- 
tions, namely that for physical purposes we need much more to incorporate the 
concept of the covariant derivative rather than the Lie derivative. There was some 
work done in this direction H75tlllltll74| . but it is a very basic sketch, rather than 
a full and satisfactory theory It is unclear (in the papers cited) how one can obtain 
for the covariant derivative operator the expected and meaningful results outside 
the class of smooth tensor fields. 

3.2.4 Colombeau algebra in the geometric approach 

Scalar special Colombeau algebra 

For arbitrary general manifold it is easy to find a covariant formulation of the spe- 
cial Colombeau algebra. At the end of the day you obtain a space of e-sequences 
of functions on the general manifold M . For the non-canonical case the defini- 
tions are similar to the non-geometric formulations, the basic difference is that Lie 
derivative plays here the role of the W 1 partial derivative. Thus the definition of 
the Colombeau algebra will be again: 

g(Q) = £ M {n)/M{n). (3.25) 

£m(^) (moderate functions) are defined as algebra of functions Vt x (0, 1] — > M, 
such that are smooth on fi (this is usually called £(0)) and for any K CC Vt we 
insist that 

VA; e Nq, 3p e N such that V£i...£ fc which are smooth vector fields, 

sup \L^...L ih f e (x)\ < 0(e- p ) as e -> 0. (3.26) 



5 We can mention also another classical formulation of distributional form fields, which comes 
from the old book of deRham |52] (it uses the expression "current"). It naturally defines the space 
of distributions to be a dual space to space of all compactly supported form fields. 
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A/"(fi) (negligible functions) we define exactly in the same analogy to the non- 
geometric formulatioro 

\/K CC Q, VA; G Nq, Vp G N, V&...& which are smooth vector fields, 

sup |%...L Sfe / e (x)| < 0(e p ) ase^O. (3.27) 

Tensor special Colombeau algebra and the embedding of distributions 

After one defines the scalar special Colombeau algebra, it is easy to define the 
generalized Colombeau tensor algebra as the tensor product of sections of a tensor 
bundle and Colombeau algebra. This can be formulated more generally |101] in 



terms of maps from M to arbitrary manifold. One can define them by changing 
the absolute value in the definition {3.27} to the expression "any Riemann measure 
on the target space". Then you get the algebra [10111103111731 : 

T C (X,Y) =T M (X,Y)/N(X,Y). (3.28) 

The tensor fields are represented when the target space is taken to be the TM 
manifoldZl It is clear that any embedding of distributions into such algebra will 
be non-canonical from various reasons. First, it is non-canonical even on R n . An- 
other, second reason is that this embedding will necessarily depend on some pre- 
ferred class of charts on M. The embedding one defines as H101I1 : 

We pick an atlas, and take a smooth partition of unity subordinate to V ai , (©.,, 
supp(6) C V a . j G N) and we choose for every j, £, G D(V aj ), such that £,- = 
1 on supp(£j). Then we can choose in fixed charts an y4 00 (R™) element p, and the 
embedding is given by 

oo 

E (((fc( e i ° <K) * *) ° ^) . ( 3 - 29 ) 



6 This version is due to [76 j. There are also different definitions: in 11011 the authors use instead 
of "for every number of Lie derivatives along all the possible smooth vector fields" the expression 
"for every linear differential operator", but they prove that these definitions are equivalent. This 
is also equivalent to the statement that in any chart holds: <f> € £a/(K") (see |101|). 

7 This is equivalent to Q{X) <E) T(X, E) tensor valued Colombeau generalized functions 11031 . 
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where ipa- is a coordinate mapping and * is a convolution. 

The equivalence relation 



Now let us define the equivalence relation in analogy to the W l case. Since in |101| 
the strongest constraint on the mollifier is taken, one would expect that strong re- 
sults will be obtained, but the definition is more complicated. And in fact, stan- 
dard results (such as embedding of smooth function multiplying distribution is 
equivalent to product of their embeddings) are not valid here 111011 . That is why 
the stronger concept of /^-association is formulated. It states that U G T c is k 
associated to function /, if 

\im L^.-.L^Ue - /) -► (3.30) 

uniformly on compact sets for all / < k. The cited paper does not contain a precise 
definition of k equivalence between two generalized functions, but it can be easily 
derived. 

The older formulation of scalar full Colombeau algebra 

If we want to get a canonical formulation, we certainly cannot generalize it straight 
from the W 1 case (the reason is that the definition of the classes A n (R. n ) is not 
diffeomorphism invariant). However, there is an approach providing us with a 
canonical formulation of generalized scalar fields [68] . The authors define the 
space £ (M) as a space of C°°(M x A°(M)), where A°(M) is the space of n-forms 
(n-dimensional space), such that J to = 1. Now the authors define a smoothing 
kernel as C°° map from 

M x I -> A°(M), (3.31) 

such that it satisfies: 

(i) \/K CC M 3e , C > Vp G K Ve < e , supp (j)(e,p) C B eC (p), 
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(ii) VK CC M, VJfe, I 6 N , VXi, • • • X k , Y x ■ ■ ■ Y x smooth vector fields, 

su P P GA', g eM ||^Yi • ■ ■ L Yl (L' Xl +L Xh )~- (L' Xk ■ ■ ■ L Xk )$(e,p)(q) \\ = 

= 0(e-(" +1 )). 

Here U is defined as: 

L' x f(p,q) = L x (p -> f(p,q)) = lf((Fl*)(p),q)\ . (3.32) 

B eC is a ball centered at C having radius e measured relatively to arbitrary Rie- 
mannian metric. Let us call the class of such smoothing kernels A°(M). Then in 
116811 classes A m (M) are defined as the set of all $ G A°(M) such that V/ G C°°(M) 
and \/K ccM (compact subset) it holds: 



sup 



/(p)- / /(g)^(c,p)(?) 



0(e m+i ). (3.33) 



Moderate and the negligible functions are defined in the following way: 

R G £(M) is moderate if \/K ccMVfceN 3iV G N \fX h ....X k (X h ...X k are 

smooth vector fields) and V$ G A°(M) one has: 



sup \\L Xl L Xk (R($(e,p),p))\\ = 0(e- n ). (3.34) 

pdK 

R G f (M) is negligible if VisT CC M, Vfc, / G N 3m G N VX l5 ...X fe (X l5 ...X fe are 
again smooth vector fields) and V$ G A m (M) one has: 

sup \\L Xl ...L Xk (Rme,p),p)\\ = 0(e l ). (3.35) 

Now we can define the Colombeau algebra in the usual way as a factor algebra 
of moderate functions over negligible functions. Scalar distributions, defined as 
dual to n-forms, can be embedded into such algebra in a complete analogy to the 
canonical IR n formulation. Also association is in this case defined in the "usual" 
way (integral with compactly supported smooth n-form field) and has for mul- 
tiplication the usual properties. However any attempt to get a straightforward 
generalization from scalars to tensors brings immediate problems, since the em- 
bedding does not commute with the action of diffeomorphisms. This problem 
was finally resolved in 1169 II . 
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Tensor full Colombeau algebra and the embedding of distributions 

The authors of reference IJ69H realized that diffeomorphism invariance can be achiev- 
ed by adding some background structure defining how tensors transport from 
point to point, hence a transport operator. Colombeau (m, n) rank tensors are 
then taken from the class of smooth maps C°° (to, q, B) having values in (T™) q M, 
where u E A°(M), q E M and B is from the class of compactly supported trans- 
port operators. After defining how Lie derivative acts on such objects and the 
concept of the "core" of a transport operator, the authors of reference 11691 define 
(in a slightly complicated analogy to the previous case) the moderate and the neg- 
ligible tensor fields. Then by usual factorization they obtain the canonical version 
of the generalized tensor fields (for more details see BUI). The canonical embed- 
ding of tensorial distributions is the following: The smooth tensorial objects are 
embedded as 

tip, tu,B) = J t(q)B(p, q)u(q)dq (3.36) 

where as expected u E A°(M), t is the smooth tensor field and B is the transport 
operator. Then the arbitrary tensorial distribution s is embedded (to s) by the 
condition 

s(u, P , B) ■ t(p) = is, B(p, .) • tip) ® w(.)) , (3.37) 

where on the left side we are contracting the embedded object with a smooth 
tensor field t, and on the right side we are applying the given tensorial distribution 
s in the variable assigned by the dot. It is shown that this embedding fulfills all 
the important properties, such as commuting with the Lie derivative operator [69 J. 
All the other results related to equivalence relation (etc.) are obtained in complete 
analogy to the previous cases. 
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The generalized geometry (in special Colombeau algebras) 



In HlOll |102[ 11031 the authors generalized all the basic geometric structures, like 



connection, covariant derivative, curvature, or geodesies into the geometric for- 
mulation of the special Colombeau algebra. That means they defined the whole 
generalized geometry 

Why is this somewhat unsatisfying? 

However in our view, the crucial part is missing. What we would like to see is an 
intuitive and clear definition of the covariant derivative operator acting on the dis- 
tributional objects in the canonical Colombeau algebra formulation, on one hand 
reproducing all the classical results, and on the other hand extending them in the 
same natural way as in the classical distributional theory with the classical deriva- 
tive operator. Whether there is any way to achieve this goal by the concept of gen- 
eralized covariant derivative acting on the generalized tensor fields, as defined in 



[103J, is unclear. Particularly it is not clear whether such generalized geometry 
can be formulated also within the canonical Colombeau algebra approach. There 
exist definitions of the covariant derivative operator within the distributional ten- 
sorial framework [75111 11] . (The reference ||75f gives particularly nice application 
of such distributional tensor theory to signature changing spacetimes.) But these 
approaches are still "classical", in the sense, that they do not fully involve the op- 
eration of the tensor product of distributional tensors. There cannot be any hope 
of finding a more appropriate, generalized formulation of classical physics with- 
out finding such a clear and intuitive definition of both, the covariant derivative 
and the tensor product. All that can be in this situation achieved is to use these 
constructions to solve some specific problems within the area of physics. But as 
we see, the more ambitious goal can be very naturally achieved by our own con- 
struction, which follows after this overview. 
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3.2.5 Practical application of the standard results 

Now we will briefly review various applications of the Coulombeau theory pre- 
sented before. 



Classical shock waves 

The first application we will mention is the non-general-relativistic one. In [48] the 
authors provide us with weak solutions of nonlinear partial differential equations 
(using Colombeau algebra) representing shock waves. They use a special version 
of the Colombeau algebra, and specifically the relation 

H n H' « — !— H', (3.38) 

n+ 1 

(which is related to mollifiers from A°(R") class). The more general analysis re- 
lated to the existence and uniqueness of weak solutions of nonlinear partial dif- 
ferential equations can be found in [124j. 

Black hole "distributional" spacetimes 

In general relativity there are many results obtained by the use of Colombeau al- 
gebras. First, we will focus on the distributional Schwarzschild geometry, which 
is analysed for example in [76] . The authors of j|76l start to work in Schwarzschild 
coordinates using the special Colombeau algebra and A°°(IR n ) classes of mollifiers. 
They obtain the delta-functional results (as expected) for the Einstein tensor, and 
hence also for the stress-energy tensor. But in Schwarzschild coordinates there 
are serious problems with the embedding of the distributional tensors, since these 
coordinates do not contain the point. As a result, if one looks for smooth em- 
beddings, one does not obtain an inverse element in the Colombeau algebra in the 
neighbourhood of for values of e close to 0. (Although there is no problem, if we 
require that the inverse relation should apply only in the sense of equivalence.) 
Progress can be made by turning to Eddington-Finkelstein coordinates |76"H. The 
metric is obtained in Kerr-Schild form, in which one is able to compute R a b, G a b 
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and hence T a b as delta-functional objects (which is expected). (The (1,1) form of 
the field equations is used since the metric dependence has a relatively simple 
form in Eddington-Finkelstein coordinates.) This result does not depend on the 
mollifier (see also |173 | 1 ), but one misses the analysis of the relation between dif- 
ferent embeddings given by the different coordinate systems lj. Even in the case 
of Kerr geometry there is a computation of y/~{jf)R IJ ' v (where g ab = g' ab + fk a k b ) 
given by Balasin, but this is mollifier dependent ffTTl [173 J. Here the coordinate 
dependence of the results is even more unclear. 

Aichelburg metric 

There exists an ultrarelativistic weak limit of the Schwarzschild metric. It is taken 
in Eddington-Finkelstein coordinates u = t + r* , v = t — r* , (where r* is tortoise 
coordinate) by taking the boost in the weak v — > c limit. We obtain the "delta func- 
tional" Aichelburg metric. Reference 111041 provides a computation of geodesies 



in such a geometry. The authors of ||104 | ] take the special Colombeau algebra (and 
take A (]R n ) as their class of mollifiers), and they prove that geodesies are given 
by the refracted lines. The results are mollifier independent. This is again ex- 
pected. Moreover, what seems to be really interesting is that there is a continuous 
metric which is connected with the Aichelburg metric by a generalized coordinate 
transformation fl00lll73l . 



Conical spacetimes 

The other case we want to mention are conical spacetimes. One of the papers 
where the conical spacetimes are analysed is an old paper of Geroch and Traschen 
||6T|. In [|6T| it is shown that conical spacetimes can not be analysed through the 
concept of gt-metric. These are metrics which provide us with a distributional 



8 It seems that the authors use relations between R^ v and components of g^ u obtained in 
Eddington-Finkelstein coordinates by using algebraic tensor computations. Then it is not obvi- 
ous, whether these results can be obtained by computation in the Colombeau algebra using the « 
relation, since in such case some simple tricks (such as substitution) cannot in general be used. 
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Ricci tensor in a very naive sense. The multiplication is given just by a simple 
product of functions defining the regular distributions. A calculation of the stress 
energy tensor was given by Clarke, Vickers and Wilson [44|, but this is mollifier 
dependent (although it is coordinate independent 1117310 . 

3.3 How does our approach relate to current theory? 

General summary How does our own approach (to be described in detail in 
the next section) relate to all what has presently been achieved in the Colombeau 
theory? We can summarize what we will do in three following points: 

• We will define tensorial distributional objects, and the basic related opera- 
tions (especially the covariant derivative). The definition directly follows our 
physical intuition (there is a unique way of constructing it). This generalizes 
the Schwartz IR n distribution theory. 

• We will formulate the Colombeau equivalence relation in our approach and 
obtain all the usual equivalence results from the Colombeau theory. 

• We will prove that the classical results for the covariant derivative opera- 
tor (as known within differential geometry) significantly generalize in our 
approach. 

The most important point is that our approach is fully based on the Colombeau 
equivalence relation translated to our language. That means we take and use 
only this particular feature of the Colombeau theory and completely avoid the 
Colombeau algebra construction. 

The advantages compared to usual Colombeau theory What are the advan- 
tages of this approach comparing to Colombeau theory? 

• First, by avoiding the algebra factorization (which is how Colombeau al- 
gebra is constructed) we fulfill the physical intuitiveness condition of the 
language used. 
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• Second, we can naturally and easily generalize the concept of the covari- 
ant derivative in our formalism, which has not been completely satisfacto- 
rily achieved by the Colombeau algebra approach.]. This must be taken as 
an absolutely necessary condition that any generalization of a fundamental 
physical language must fulfill. 

• It is specifically worth discussing the third advantage: Why is the classi- 
cal approach so focused on Colombeau algebras? The answer is simple: 
We want to get an algebra of C°°(M) functions as a subalgebra of our al- 
gebra. (This is why we need to factorize by negligible functions, and we 
need to get the largest space where they form an ideal, which is the space 
of moderate functions.) But there is one strange thing: all our efforts are 
aimed at reaching the goal of getting a more rich space than is provided by 
the space of smooth functions. But there is no way of getting a larger dif- 
ferential subalgebra than the algebra of smooth functions, as is shown by 
Schwartz impossibility result. That is why we use only the equivalence re- 
lation instead of straight equality. But then the question remains: Why one 
should still prefer smooth functions? Is not the key part of all the theory 
the equivalence relation? So why is it that we are not satisfied with the way 
the equivalence relation recovers multiplication of the smooth objects (we 
require something stronger), but we are satisfied with the way it recovers 
multiplication within the larger class? Unlike the Colombeau algebra based 
approach, we are simply taking seriously the idea that one should treat all 
the objects in an equal way. This means we do not see any reason to try 
to achieve "something more" with smooth objects than we do with objects 
outside this class. And the fact that we treat all the objects in the same way 
provides the third advantage of our theory; it makes the theory much more 
natural than the Colombeau algebra approach. 

• The fourth advantage is that it naturally works with the much more general 



9 As previously mentioned, the current literature lists some ambiguous attempts to incorporate 
the covariant derivative, but no satisfactory theory. 
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(and for the language of distributions natural) concept of piecewise smooth 
manifolds, so the generalization of the physics language into such concep- 
tual framework will give it much higher symmetry. 

• The fifth and last, "small" advantage comes from the fact that by avoiding 
the Colombeau algebra construction we automatically remove the problem 
of how to canonically embed arbitrary tensorial distributions into the alge- 
bra. But one has to acknowledge that this problem was already solved also 
within the Colombeau theory approach [69J. 

The disadvantages compared to usual Colombeau theory What are the disad- 
vantages of this approach compared to Colombeau theory? A conservative person 
might be not satisfied with the fact that we do not have a smooth tensor algebra as 
a subalgebra of our algebra. This means that the classical smooth tensorial fields 
have to be considered to be solutions of equivalence relations only (as opposed 
to the classical, "stronger" view to take them as solutions of the equations). But 
in my view, this is not a problem at all. The equivalence relations contain all the 
classical smooth tensorial solutions (equivalently smooth tensor fields), for the 
smooth initial value problem. So for the smooth initial value problem the equiv- 
alence relations reduce to classical equations. We will suggest how to extend the 
initial value problem for larger classes of distributions and show that it has unique 
solutions. It is true that the equivalence relations might have also many other 
(generally non-linear) solutions in the space constructed]!!] but the situation that 
there exist many physically meaningless solutions is for physicists certainly noth- 
ing new. Our previous considerations suggest that we shall look only for distribu- 
tional solutions (they are the ones having physical relevance) where the solution 
is provided to be unique (up to initial values obviously), if it exists. (This will also 
recover the common, but also many "less" common physical results |173|| .) 



10 



By "solution" we mean here any object fulfilling the particular equivalence relations. 
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3.4 New approach 

3.4.1 The basic concepts/definitions 

Before saying anything about distributional tensor fields, we have to define the 
basic concepts which will be used in all the following mathematical constructions. 
This task is dealt with in this section, so this is the part crucially important for un- 
derstanding all the subsequent theory An attentive reader, having read through 
the introduction and abstract will understand why we are particularly interested 
in defining these concepts. 

Definition of (M,A) 

Definition 3.4.1. By piecewise smooth function we mean a function from an open 
set r^x(C R 4 ) — > W 71 such, that there exists an open set fi 2 (in the usual "open ball" 
topology on M 4 ) on which this function is smooth and Q 2 — fli \ Vl' (where f2' has 
a Lebesgue measure 0). 

Take M as a 4D paracompacO. Hausdorff locally Euclidean space, on which 
there exists a smooth atlas S. Hence (M , S) is a smooth manifold. Now take 
a ordered couple (M, A), where A is the maximal atlas, where all the maps are 
connected by piecewise smooth transformations such that: 

• the transformations and their inverses have on every compact subset of R 4 
all the first derivatives (on the domains where they exist) bounded 
(hence Jacobians, inverse Jacobians are on every compact set bounded), 

• it contains at least one maximal smooth subatlas S C A, 
(coordinate transformations between maps are smooth there). 



11 We will use specifically 4-dimensional manifolds, but one can immediately generalize all the 
following constructions for n-dimensional manifolds. 



CHAPTER 3. MULTIPLICATION OF TENSORIAL DISTRIBUTIONS 126 

Notation. The following notation will be used: 

• By the letter S we will always mean some maximal smooth subatlas of A. 

• Every subset of M on which there exists a chart from our atlas A, we call 
Vl C h- An arbitrary chart on Vl C h from our atlas A is denoted Ch{Vt C h)- 

Notation. Take some set Vt C h- Take some open subset of that set Q! C flch- Then 
Ch(Vtch)\n' is defined simply as Ch(Q,'), which is obtained from Ch{Vt C h) by limit- 
ing the domain to 0! . 

Definition of "continuous to the maximal possible degree" 

Definition 3.4.2. We call a function on M continuous to the maximal possible 
degree, if on arbitrary Vt of Lebesgue measure it is only in such casesll^ 

a) either undefined, 

b) or defined and discontinuous, 

in which there does not exist a way of turning it into a function continuous on Vt 
by 

• in the case of a) extending its domain by the Vt set, 

• in the case of b) re-defining it on Vt. 

Jacobians and algebraic operations with Jacobians 

Now it is obvious that since transformations between maps do not have to be 
everywhere once differentiable, the Jacobian and inverse Jacobian may always be 



12 The expression "Lebesgue measure set" will have in this chapter extended meaning. It refers 
to such subsets of a general manifold M that they have in arbitrary chart Lebesgue measure 0. 
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undefined on a set having Lebesgue measure 0. Now if we understand product in 
the sense of a limit, then the relation 

■K {J^Tu = %, (3-39) 

for example, might hold even at the points where both Jacobian and inverse Jaco- 
bian are undefined. This generally means the following: any algebraic operation 
with tensor fields is understood in such way, that in every chart it gives sets of 
functions continuous to a maximal possible degree. From this follows that the 
matrix product (|3.39|) must be, for jjl — v, equal to 1 and, for /i^v, 0. 



Tensor fields on M 

We understand the tensor field on M to be an object which is: 

• Defined relative to the 1-differentiable subatlas of A everywhere except for 
a set having Lebesgue measure (this set is a function of the given 1-dif- 
ferentiable subatlas). 

• In every chart from A it is given by functions continuous to a maximal pos- 
sible degree. 

• It transforms ^Vt Ch between charts Chi(Vt C h), Ch 2 (^ch) £ Am the tensorial 
way 



Tt.cH, = ^ (J~ l t ■ ■ ■ Tl;; Chl a.e.,13 (3.40) 



Ch) 



where J£ is the Jacobian of the coordinate transformation from Chi(Vt 
to Ch 2 (£lch), and T^'" Chi ,T^"" Ch are tensor field components in charts Chi, 
Ch 2 . As we already mentioned: If T^"" Ch is at some given point undefined, 
so in some chart Chi(Vt C h) the tensor components do not have a defined 



13 This expression means "almost everywhere", that is, "everywhere apart from a set having 
Lebesgue measure 0". 
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limit, this limit can still exist in C/i 2 (^c/i)/ since Jacobians and inverse Jaco- 
bians of the transformation from Chi{Vt C h) to Ch 2 (ft>ch) might be undefined 
at that point as well. This limit then defines T^"" Ch2 at the given point. 

Important classes of test objects 

Notation. The following notation will be used: 

• We denote by C P (M) the class of 4-form fields on M, such that they are 
compactly supported and their support lies within some Vt C h- For such 4- 
form fields we will generally use the symbol to. 

• For the scalar density related to u E C P (M) we use always the symbol u/. 

• By C p {Vtch) we mean a subclass of C P (M), given by 4-form fields hav- 
ing support inside VL C h- Note that only the C p (Vt C h) subclasses form linear 
spaces. 

• Take such maximal atlas S, (3S C S C A), that there exist 4-forms from 
C P (M), such that they are given in this atlas by everywhere smooth scalar 
density ui'. ("Maximal" here means that these 4-forms have in every chart 
outside this atlas non-smooth scalar densities.) By C P ,AM) we mean a class 
of all such elements from C P (M), that they have everywhere smooth scalar 
density in S. 



The letter S will from now on be reserved for maximal atlases defining 
C P ,AM) classes. 

C P (M) is defined as: C P {M) = \J S C P S (M). 

)) means C p 
port inside the given VL C h- 



C s(S) ^h) (or C p (n Ch )) means C p {S) (M) (or C P (M)) element having sup- 
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Topology on C^ S) (Q C h) 

Consider the following topology on each C^.xAQch)'- A sequence from ui n e 
Cg,gJQch) converges to an element u from that set if all the supports of u n lie in 
a single compact set, and in any chart Ch(Vt C h) £ S, for arbitrary k it is true that 

d k u>'(x l ) • c i . d k uj'(x 1 ) 

^rf^jL converges uniformly to g^±. 

3.4.2 Scalars 

This section deals with the definition of scalar distributions as the easiest par- 
ticular example of a generalized tensor field. The explanatory reasons are the 
main ones why we deal with scalars separately, instead of taking more "logical", 
straightforward way to tensor fields of arbitrary rank. 

Definition of D\M), and hence of linear generalized scalar fields 

Definition 3.4.1. We say that B, being a function that maps some subclass of 
C P (M) to W l , is linear, if the following holds: Take such co 1 and u 2 from the do- 
main of B, that they belong to the same class C p (VL C h)- Whenever the domain 
of B contains also their linear combination AiWi + A2W2/ where Ai, A2 £ K, then 

B{\iWi + A2W2) = \iB(u){) + \ 2 B(u 2 ). 

Definition 3.4.2. Now take the space of linear maps F — ¥ R, where F is such set 
that 3S, such that C p '= AM) C F C C P (M). These linear maps are also required 
to be for every Vt C h on C P ,(Vt C h) continuous, (relative to the topology taken in 
section |3.4.1[) . Call set of such maps D'(M), or in words, the set of linear generalized 
scalar fields. 

Important subclasses of D'{M) 

Notation. The following notation will be used: 

• Now take a subset of D'(M) given by regular distributions defined as in- 
tegrals of piecewise continuous functions (everywhere on C P (M), where it 
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converges). We denote it by D' E (M)n 

• Take such subset of D' E (M) that there 3S in which the function under the 
integral is smooth. Call this class D' S (M). 

• Now take subsets of D'(M) such that they have some common set 
U n Cf.= AM) belonging to their domains and are Vf^, Vn continuous on 
C* s Mch)> Denote such subsets by D' , e AM). Obviously T E D' « AM) 
means T E n n L>'~ AM). 

(On) 

• By £>' a x(M) we mean objects such that they belong to £>' * JM) and 
their full domain is given as U„Cf, = AM). 

• If we use the notation D' = (M), we mean objects defined by integrals of 
piecewise continuous functions, with their domain being the class U n Cf , = . (M) . 

• By using D', - AM) we automatically mean subclass of D' - AM), 

J ° S{U n S n o) y > J E(U„S„o) y ' 

such that it is given by an integral of a smooth function in some smooth 
subatlas S C U n <S n . 

D' A (M), hence generalized scalar fields 

Notation. Let us for any arbitrary set D'=. (M]l3 construct, by the use of point- 
wise multiplication of its elements, an algebra. Another way to describe the alge- 
bra is that it is a set of multivariable arbitrary degree polynomials, where different 
variables represent different elements of D'.AM). Call it D', )A (M). 

By pointwise multiplication of linear generalized scalar fields Bi, B 2 E 
D'g. (M) we mean a mapping from u into product of the images (real numbers) 
of the Bi,B 2 mappings: {B x ■ B 2 )(uj) = Bi(oj) ■ B 2 (tu). The domain on which 
the product (and the linear combination as well) is defined is an intersection of 



14 Actually, it holds that if and only if the function is integrable in every chart on every compact 
set in W 1 , then this function defines a regular D'(A1) distribution and is defined at least on the 

whole Cf(M) class - 

15 The union is here trivial, it just means one element S. 
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domains of B x and B 2 (trivially always nonempty containing C^=JM) at least). 
Note also that the resulting arbitrary element of D' - AM) has in general all the 
properties defining D'(M) objects, except that of being necessarily linear. 

Definition 3.4.3. The set of objects obtained by the union U^D'-(M) we denote 
D' A {M), and call them generalized scalar fields (GSF). 

Topology on ^(M) 

If we take objects from D'* AM) and we take a weak (a—) topology on that set, 
we know that any object is a limit of some sequence from D'.AM) objects from 
that set (they form a dense subset of that set). That is known from the classical 
theory. Such a space is complete. 



3.4.3 Generalized tensor fields 

This section is of crucial importance. It provides us with definitions of all the basic 
objects we are interested in, the generalized tensor fields and all their subclasses 
of special importance as well. 

The class D'™(M) of linear generalized tensor fields 

First let us clearly state how to interpret the Jff Jacobian in all the following def- 
initions. It is a matrix of piecewise smooth functions f^ \ Q 2 -^ K, f^ being a 
open subset of M 4 and f2 2 having Lebesgue measure 0. Let it represent transfor- 
mations from Chi(Vt C h) to Ch 2 (Vt C h)- We can map the Jacobian by the inverse of 
the Chi(Vt C h) coordinate mapping to VL C h and it will become a matrix of functions 
Vtch \ fi' — )> IR, Vt' having Lebesgue measure 0. The object J/f • u> is then under- 
stood as a matrix of 4-forms from C p {Vt C h), which also means that outside Vtch we 
trivially define them to be 0. 

Definition 3.4.1. Take some set F, F ± C F C F 2 , where F 1 us such set that 3<S and 
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3<SC<S (S is some maximal smooth atlas) defining it as: 

F 1 = Un Ch {{oo,Ch(Q Ch )) : u E C P { ^(Q Ch ), Ch(Q Ch ) E s} . 
F 2 is defined as: 

f 2 = un Ch {(u,ch(n Ch )) : u e c p (n Ch ), ch(n Ch ) eA)}. 

By a D'™(M) object, the linear generalized tensor Held, we mean a linear mapping 
from F — >■ IR 4 " 1 " for which the following holds: [Jj 

• VQch it is VChkC&ch) ES cS continuous on the class C P . s J£lch)- (Both S, 

0(0) 

S are from the definition of Fi .) 

• This map also ^Vt C h transforms between two charts from its domain, 

Ch^Vtch), Ch 2 (n Ch ),as: 

T^:2{Ch u u) = T-;i {ch 2 , j*...Ji(r x )l....(r x y 6 ■ u) . 

• The following consistency condition holds: If Vt' Ch C ficfc, then T^'f gives 
on u) x C7i(fickW , w E C p (Q' Ch ) the same results^ as on w x Ch(Q C h)- 

We can formally extend this notation also for the case m = n = 0. This means 
scalars, exactly as defined before. So from now on m, n take also the value 0, 
which means the theory in the following sections holds also for the scalar objects. 



16 We could also choose for our basic objects maps taking ordered couples from C (Ctch) x 
Ch(Q! Ch ), (flch 7^ ft'ch)- The linearity condition then automatically determines their values, since 
for u e C p (Qch), whenever it holds that Q' Ch n supp(w) = {0}, they must automatically give 
for any chart argument. Hence these two definitions are trivially connected and choice between 
them is just purely formal (only a matter of "taste"). 

17 By the "same results" we mean that they are defined on the same domains, and by the same 
values. 
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Important subclasses of D'™(M) 

Notation. The following notation will be used: 

• By a complete analogy to scalars we define classes D'^ E (M): On arbitrary 
Vl C h, being fixed in arbitrary chart Chi(Qch) £ A we can express it in an- 
other arbitrary chart Ch 2 (£lch) £ A, as an integral from a multi-index matrix 
of piecewise continuous functions on such subset of C P (M), on which the 
integral is convergent!^. 

• Analogously the class D'^(M) c D'^(M) is defined by objects which can, 
for some maximal smooth atlas <S, in arbitrary charts E3 Chi(Q C h), Ch 2 (Qch) £ 
S, be expressed by an integral from a multi-index matrix of smooth func- 
tions. 

• D'"} - (M) means a class of objects being for every £l Ch in every Ch(Vt C h) £ 
At(Si) continuous on Cg,gJQch)- (Here At(S) stands for a map from atlases 
S to some subatlases of A.) 

• D'™ , . i/t =. ,(M) means a class of objects from D'™ , . xl5 ,AM) having as their 
domain the union 

u Qch u* {(u,ch(n Ch )) : w e cf^^), C7i(ft c/I ) e At(s,)} . 

• If we have classes D /r " ^ C ~JM) and D"? .„*, .(M) where At(SA = S t C 



<S,, we use the simple notation D'^^M), D% lSlo) (M) 



20 



18 Actually we will use the expression "multi-index matrix" also later in the text and it just means 
specifically ordered set of functions. 

19 The first chart is an argument of this generalized tensor field and the second chart is the one 
in which we express the given integral. 

20 We have to realize that the subatlas S n specifies completely the atlas S n , since taking forms 
smooth in S n determines automatically the whole set of charts in which they are still smooth. This 
fact contributes to the simplicity of this notation. 
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Definition of D'™ A (M), hence generalized tensor fields 

Definition 3.4.2. Now define D'^ S , A (M) to be the algebra constructed from the ob- 
jects D'^xs-s (M) by the tensor product, exactly in analogy to the case of scalars (this 
reduces for scalars to the product already defined). The object, being a result of 
the tensor product, is again a mapping V — v IR 4 ™ ", defined in every chart by com- 
ponentwise multiplication. Now denote by D^ A (M) a set given as UsD'^ s ^ A (M), 
meaning a union of all possible D'^ S ^ A (M). Call the objects belonging to this set 
the generalized tensor fields (GTF). 

Notation. Furthermore let us use the same procedure as in the previous defini- 
tion, just instead of constructing the algebras from the classes D'£? S JM), we now 
construct them only from the classes D'™ A( = . . (M) fl D%? s ^ (M), (again by tensor 
product). For the union of such algebras we use the notation D'™ .* . (M). 

Definition of T— objects, their classes and algebras 

Notation. Now let us define the generalized space of objects T l (M). (For ex- 
ample the Christoffel symbol would fall into this class.) These objects are defined 
exactly in the same way as D'™(M) (m + n = I) objects, we just do not require that 
they transform between charts in the tensorial way, (second point in the definition 
of generalized tensor fields). 
Note the following: 

• The definition of T l (M) includes also the case m — 0. Now we see, that the 
scalars can be taken as subclass of T°(M), given by objects that are constants 
with respect to the chart argument. 

• Note also that for a general T m (M) object there is no meaningful differenti- 
ation between "upper" and "lower" indices, but we will still use formally 
the Tj};;; notation (for all cases). 
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Notation. In the same way, (by just not putting requirements on the transforma- 
tion properties), we can generalize the classes 

• D' m - (M) to T m+n - (M), 

• D ':tu tM s t)o) W to r ru - 5 - )o) (M), 

• Kl lSl) W to r«(M), 

• £>^(M) to F£ +n (M), 

• D^(M) to r™ + "(M). 

It is obvious that all the latter classes contain all the former classes as their 
subclasses, (this is a result of what we called a "generalization"). 

Note that when we fix T^(M) objects in arbitrary chart from A, they must be 
expressed by integrals from multi-index matrix of functions integrable on every 
compact set. In the case of the D'^ E (M) subclass it can be required in only one 
chart, since the transformation properties together with boundedness of Jacobians 
and inverse Jacobians, provide that it must hold in any other chart from A. The 
specific subclass of T^(M) is r™ (M), which is a subclass of distributions given 
in any chart from A, (being an argument of the given T— object), by integrals 
from multi-index matrix of smooth functions (when we express the integrals in 
the same chart, as the one taken as the argument). T™, UnSn0 JM) stands again for 
r™(M) objects with domain limited to 

u Qch u n {(u,Ch(tt Ch )) : w e C^ n) (Q C h), Ch(n Ch ) e S n } , 
where S n is given by the condition S n C S n . 
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Notation. Take some arbitrary elements T£- e T r ^(M), u> E C p (Vt C h), and 
Ch k (il C h) £ A. The Tjf-(Ch k) uj) can be always expressed as f Q T^-(Ch k ) ■ u. 
Here T^—(Chk) appearing under the integral denotes some multi-index matrix of 
functions continuous to a maximal possible degree on Qch' For T£- E D'^(M) 
the T£-(Chk) multi-index matrix components can be obtained from a tensor field 
by: 

• expressing the tensor field components in Ch k (ilch) on some subset of R 4 , 

• mapping the tensor field components to Vl C h by the inverse of Ch k (Vt C h)- 
Furthermore tu'(Ch k ) will denote the 4-form scalar density in the chart Ch k (Vtch)- 

Topology on r-^ o) (M) 

If we take the class of TRi s ) (M), and we impose on this class the weak (point or 
a—) topology then the subclass of ^T, nSn )(M) defined as T™, UnS ,(M) is dense 
^ r^ o) (M). The same holds for D% lSl0) (M) and D>™ {UiSio) (M). 

Definition of contraction 

Definition 3.4.3. We define the contraction of a T™(M) object in the expected 
way: It is a map that transforms the object T;£~ E T™(M) to the object T;;;£;; E 

r r 2' 2 {M). 

Now contraction is a mapping D'™(M) -± D^{M) and T m (M) -* T m - 2 (M), 
but it is not in general the mapping D'™ A (M) -± £>^(M), only D'™ A (M) -> 

Interpretation of physical quantities 

The interpretation of physical observables as "amounts" of quantities on the open 
sets is dependent on our notion of volume. So how shall we get the notion of 
volume in the context of our language? First, by volume we mean a volume of 
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an open set. But we will consider only open sets belonging to some VL C h- So take 
some Vtch and some arbitrary Q! c &ch- Let us now assume that we have a 
metric tensor from D'^(M). This induces a (volume) 4-form. Multiply this 4- 
form by a noncontinuous function xci' defined to be 1 inside Q' and everywhere 
else 0. Call it u>q>. Then by volume of an open set Q' we understand: J uqi. Also 
u)q> is object from C P (M) (particularly from C p (Vt C h))- The "amounts" of physical 
quantities on Q! we obtain, when the D'£\(M) objects act on uq>. 

3.4.4 The relation of equivalence («) 

This section now provides us with the fundamental concept of the theory, the con- 
cept of equivalence of generalized tensor fields. Most of the first part is devoted 
to fundamental definitions, the beginning of the second part deals with the basic, 
important theorems, which just generalize some of the basic Colombeau theory 
results to the tensor product of generalized tensor fields. It adds several impor- 
tant conjectures as well. The first part ends with the subsection "some additional 
definitions" and the second part with the subsection "some additional theory". 
They both deal with much less central theoretical results, but they serve very well 
to put light on what equivalence of generalized tensor fields means "physically". 

The necessary concepts to define the equivalence relation 

Notation. Take some subatlas of our atlas, this will be a maximal subatlas of 
charts, which are maps to the whole of R 4 . Such maps exist on each set Vtch and 
they will be denoted as Ch'(Vt C h)- We say that a chart Ch'(VL C h) is centered at the 
point q E Vtch, if this point is mapped by this chart to (in IR 4 ). We will use the 
notation Ch'(q, Vtch)- 



Notation. Take some Vt C h, q E Vt Ch and Ch'(q,£l C h) £ <S. The set of 4-forms 
u € G A n (S, Ch'(q, Qch)) is defined in such way that lu € e C P ( ^Jilch) belongs to 
this class if: 
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a) in the given Ch'(q, Q' Gh ), Ve it holds that: 

Ich'(W )(IL x i l ) u 'e( x ) ^ x = $ko, Y,i k i = k, k<n, fiGN, 

b) the dependence on e is in Ch'{q, Sl' C h) given as e~ 4 u/(|). 

Notation. Take an arbitrary q, Vt C h (q £ ^Ch), Ch'(q, Vt C h) G S and some nat- 
ural number n. For any u e E A n (Ch'(q,Vtch),S) we can, relatively to Ch'(Vtch), 
define a continuous set of maps (depending on the parameter y) 

A n (Cti(q,n C h),S) -> C% (S) (n C h), (3-41) 

such that they are, on Vt C h and in Ch'{Vt C h)r given as u'(^)e~ 4 — ¥ w / (^ £ )e" 4 . 
(To remind the reader u/ is the density expressing u in this chart.) This gives us 
(depending on the parameter y E M 4 ) various C^.^iflch) objects, such that they 
are in the fixed Ch'(flch) expressed by u'i^^-) e~ 4 dx 1 f\ ... f\dx 4 . Denote these 
4-form fields by £o e {y). 

Notation. Now, take any Tjf- E T™ - AM). By applying it in an arbitrary 
fixed chart Ch k {Vt C h) € At{S,&ch) on the 4-form field & e (y), obtained 
from u e E A n (Ch'(q, Q C h),S) through the map (13.41b , we get a function JR 4 ->■ 
M 4 " 1 ".Asa consequence, the resulting function depends on the following objects: 

Tfc (E IT(M)) , oo e (E A n (q, n Ch , Ch\n Ch ))) and Ch k (n Ch ). We denote it by: 
Fit;- {T£:,S, Qch, Ch'(q, Q Ch ),n, u e (y), Ch k (Q Ch )). 

Definition of the equivalence relation 

Definition 3.4.1. Bfc,T£;; E IT(M) are called equivalent (B£y « T£;-), if: 

• they belong to the same classes IT. t( s ))A (M), 

• Vn c/l , Vg (qEn Ch ), VCh'(q,n Ch )ES ( such that B£;;, Tfc E 

r™ t(5 - )M (M) ), \fCh(Q C h) e At(S,Q C h) 3n, such that Vw e G 
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A n (Ch'(q, Vtch), S) and for any compactly supported, smooth function 

R 4 -)- R, (j), it holds: 

lim / \F^ r (B^;,q,n' Ch ,Ch'(n' Ch ),n,Co e (y),Ch(n Ch )) 

e ^° JIR4 [ 

- F^;\T^ : , q ,a ch ,Ch\n' Ch ),nMv),Ch{n Ch ))\ ■ <j>(y)d*y = 0. (3.42) 



Note that for ££*;;, C£J •••, D^'", T^-- having the same domains and being from the 
samel"™ ^(M) classes, it trivially follows that: T£" « -B^'.'.', C£y « £>£- implies 
AiT£;- + \ 2 C£;; « AiB£;; + A 2 £>£;; for Ai, A 2 e R. 

Notation. Now since we have defined an equivalence relation, it divides the ob- 
jects T™(M) naturally into equivalence classes. The set of such equivalence classes 
will be denoted as T A (M). Later we may also use sets of more limited classes of 
equivalence D'^\(M), D'^ A (M) (etc.), which contains equivalence classes (only) 
of the objects belonging to D'™ A (M), D'™ EA (M) (etc.). 

Notation. In some of the following theorems, (also for example in the definition 
of the covariant derivative), we will use some convenient notation: Take some 
object B^- e T^(M). The expression T£["(B%"/u) will be understood in the fol- 
lowing way: Take Chk(Qch) x oj (to E C p {Vtch)) from the domain of T£". Then 
Bg"'(Chk) -w isa multi-index matrix of C p {Vtch) objects. This means that outside 
Vtch set they are defined to be trivially 0. We substitute this multi-index matrix of 
C p (Vt C h) objects to T£", with the chart Ch k (Vt C h) taken as the argument. 

Relation to Colombeau equivalence 

A careful reader now understands the relation between our concept of equiva- 
lence and the Colombeau equivalence relation. It is simple: The previous defi- 
nition just translates the Colombeau equivalence relation (see [47]) into our lan- 
guage and the equivalence classes will naturally preserve all the features of the 
Colombeau equivalence classes (this will be proven in the following theorems). 
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Some additional definitions (concepts of associated field and A class) 

We define the concept (of association) to bring some insight to what our concepts 
mean in the most simple (but most important and useful) cases. It enables us to 
see better the relation between the calculus we defined (concerning equivalence) 
and the classical tensor calculus. It brings us also better understanding of what 
equivalence means in terms of physics (at least in the simple cases). It just means 
that the quantities might differ on the large scales, but take the same small scale 
limit (for the small scales they approach each other). 

Definition 3.4.2. TakeT£;- e FJ(M). Assume that: 

a) V<S, such that T£; G Tf S)A (M), \/Vt Ch and VCh k (VL Ch ) e S 

3 Vt C h \ Q'(Ch k ), (the set Vl'(Ch k ) being in any Lebesgue measure), 
such that Vg G VL Ch \ W(Ch h ), VCh'(q,n Ch ) G S c S 3n, 
such that Vw £ G A n (Cti(q,tt Ch ),S) 

3 ]imTf"(Ch k ,u e ). (3.43) 



b) The limit QM is VCh'(q,n C h) 6 S c S, Vw e G A n (C/i'(g,fich),5) the 
same. 



If both a) and b) hold, then the object defined by the limit (|3.43|) is a mapping: 



Ch{Vt Ch ){e S)) x fi Ch \ fi'(CTi) -> M 4 " l+ ". (3.44) 

We call this map the field associated to T^ ;;• G T™(M), and we use the expression 
A S (T£\"). (It necessarily fulfills the same consistency conditions for &Q h c Vl 2 Ch as 
the r%{M) objects.) 

Definition 3.4.3. Denote by A c r™. - (M) a class of objects, such that each 
TfceA can be Vfi Ch/ Vn, Vw eC^tOc), VS c «S n n Ai(<S n ), 
VChk(£lch) G 5 expressed as a map 

(w,C7i fc )->. / T£;;(Ch h )'U, (3.45) 
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where for T£-(Ch k ) holds the following: In each chart from S for every point z , 
where T£;;(Ch k ) is continuous 3 5 > 0, 3^>0, such that Ve (0 < e < 5)) and 
for arbitrary unit vector n (in the Euclidean metric on M 4 ) 

T»:iCh k , z ) - K£;;e < T£:(Ch k , z + ne) < T^~(Ch k , z„) + K£;;e. (3.46) 



Notation. Take from ([3.4.3]) arbitrary fixed T£;;*, Qch and Ch k (Vt C h)- By the 
notation Q(Ch k ) C fie/! we denote a set (having Lebesgue measure 0) on which is 
T£"(Chk) discontinuous. 

Reproduction of the basic results by the equivalence relation 

Theorem 3.4.1. Any class T!" . , = , (M) contains maximally one linear element. 

Proof. We need to prove that there do not exist such two elements of T m (M), 
which are equivalent. Take two elements B and T from the class T™ .* . . .(M) 
(both with the given domains and continuity). Take arbitrary Vt C h, arbitrary S 
from their domains, and arbitrary Ch'{Vtch) £ S. Map all the C^.gJQch) objects 
to smooth, compact supported functions on R 4 through this fixed chart mapping. 
Now both B and T give, in fixed but arbitrary Ch(Vl Ch ) e At(S) linear, continuous 
maps on the compactly supported smooth functions. (The only difference from 
Colombeau distributions is that it is in general a map to K m , so the difference is 
only "cosmetic".) 

Now after applying this construction, our concept of equivalence reduces for 
every Ch(Vt C h) £ At(S) to Colombeau equivalence from [47]. The same results 
must hold. One of the results says that there are no two distributions being equiv- 
alent. All the parameters are fixed but arbitrary and all the 4-f orms from domains 
of B and T can be mapped to the M 4 functions for some proper fixing of VL C h 
and S. Furthermore, the C$(M) 4-f orms are arguments of B and T only in the 
charts, in which B and T were compared as maps on the spaces of IR 4 functions. 
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So this "arbitrary chart fixing" covers all their domain. As a result B and T must 
be identical and that is what needed to be proven. □ 

Theorem 3.4.2. Any class of equivalence T EA (M) contains maximally one linear 
element. 

Proof. First notice that the elements of Y EA (M) are continuous and defined on 
every C P x. (M) in every chart from A, so they are required to be compared in 
any arbitrary chart from A. By taking this into account, we can repeat the pre- 
vious proof. There is one additional trivial fact one has to notice: the Cg(M) 
domain also uniquely determines how the T E (M) element acts outside Cg(M). 
So if B, T e T EA (M) give the same map on Cg(M), they give the same map 
everywhere. □ 

Theorem 3.4.3. The following statements hold: 

a) Take T£;£ e T a EA (M) such that m Ch , VCh k {Sl Ch ) e A and Vw e C p (Q Ch ) 
T^'f is defined as a map 

(Ch k ,u)^ f Tf„-jCh k ) u ... f T%pJCh k )u. (3.47) 

Then the class of equivalence T EA (M), to which T^- belongs, contains a 
linear element defined (on arbitrary Vt C h) as the map: Vw e C p (Vtch), 
VCh k (n Ch ) e A, 

(Ch k ,u) -> / TftJCh k )...Tfo-JCh k ) u, (3.48) 

if and only if \/Vt Ch , \/Ch k (Q Ch ) e A 3Chi(Q C h) € A such that 

/ Tf v -jCh k )...T^(Ch k ) d 4 x (3.49) 

J Chi{Q C h)\n' 

converges on every compact set Q! C Vtch- 

The same statement holds, if we take instead of T EA (M) its subclass D'£ EA (M) 
and instead of the equivalence class T EA (M), the equivalence class D'£ EA (M). 
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The same statement also holds if we take instead of T^ A (M) and D'^, A {M) 

classes, the classes r^ (UjAt(jfi)o) (Af) and ^(u^^jW/ ( with the exce P" 
tion that the given convergence property shall be considered only for charts 
from UiAt(Si) ). 

b) For any distribution Aft" G F a s{UnSn \(M), and an element T£;; G 
^0j„5„) (-^0/ w e have that A^; ;T^; ■ is equivalent to an element of 
r ri£*)( M )' (and for subclasses £>§ Un5n) (M) c r jj^ (M) and 
A'l(u n 5„ )( M ) e r^ niSi0) (M) it is equivalent to an element of 
D't+us o) (M)). The element is on its domain given as the mapping 

{u,Ch k )^T2;;;{Af;;u). (3.50) 

c) For any tensor distribution Aft- G T a s (M) and an element Tjf- G 
IT, 5 )(M), we have that Ag";T£;" is equivalent to an element of 
^(u + 5 o)(M)> Th e element is on its domain given as mapping 

{u,Ch k )^T£;;{Aft : ;u). (3.51) 

Proof. a) Use exactly the same construction as in the previous proof. For ar- 
bitrary Vl C h and arbitrary Ch k (Vt C h) G A, we see that T£- is for every 
iu G Cg,xJQch) given by (|3.47[) (it is continuous in arbitrary chart on every 



Cg.zAtich))- We can express the map (|3.47[) in some chart Chi{Vt C h) as 

(Ch k ,uj') 
-> / Tf; (Ch k )u'd 4 x... f T^ (Ch k )u'd 4 x. (3.52) 

Jchi{Slch) JChtiQch) 

Then it is a result of Colombeau theory that if 

(C7i*,u/) -> / Tf; (Ch k )...Tjf a $ {Ch k )u'#x 

JCh^Uch) 



(3.53) 



is defined as a linear mapping on compactly supported, smooth R 4 func- 
tions uo' (in our case they are related by Chi{Vt C h) to given C$(M) objects), 
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it is equivalent to (|3.52|) . Now everything was fixed, but arbitrary, so the 
result is proven. From this proof we also see that the simple transforma- 
tion properties of the D'™(M) objects^ are fulfilled by the map (|3.53[) if the 
objects multiplied are from D'£ A (M). So the second result can be proven im- 
mediately. The last two results concerning the classes with limited domains 
trivially follow from the previous proof. 

b) is proven completely in the same way, we just have to understand that be- 
cause of the "limited" domain of the D' b a s , UnS , (M) objects, we can effec- 
tively use the concept of smoothness in this case. 

c) is just the same as b), the only difference is that the domain of the product is 
limited because of the "second" term in the product. 

□ 

Note that this means that tensor product gives, on appropriate subclasses of 
DT EA (M), the mapping D' bEA (M) x D>™ EA -► & b a + n m EA (M). It also means that this 
procedure gives, on appropriate subclasses of T EA (M), the mapping T EA (M) x 
f% A ->■ f^ m (M). The disappointing fact is that this cannot be extended to D'™ A (M). 

Theorem 3.4.4. Take T£~ e IT, Sn o) AM), Bfc E IT, JM) and Lf» E 



r n s( u, Sl o)(M). Then T£; « B£ V implies^ (L ® T)?"£; « (L <g> B) 



\Ot...H- _. IT <r> D\0-M- 



',9...k.. V-^ <y ^ >p...v.. 



Proof. Use the same method as previously. It trivially follows from the results 
of Colombeau theory (especially from the theorem saying that if a Colombeau 
algebra object is equivalent to a distribution, then after multiplying each of them 
by a smooth distribution, they remain equivalent). □ 

Theorem 3.4.5. Contraction (of pi and u index) is always, for such objects T '•;;£;;• E 
D'™ EA {M) that they are equivalent to some linear element, a map to some element 
of the equivalence class from T EA n ~ 2 (M). The equivalence class from 



21 We include also the scalar objects here. 

22 It is obvious that we can extend the definition domains either of T£'" and Bfr", or of Lfr 
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f^ n_2 (M) is such, that it contains (exactly) one element from D'™S XE {M) and this 
element is defined as the map: Vf^, ^Ch k (Vt C h) £ A, Vlu E C p {Vt C h), 

(u,Ch k )^ [ T;;£;;(Ch k ) u. (3.54) 

JOch 
Proof. The proof trivially follows from the fact that contraction commutes with 
the relation of equivalence (this trivially follows from our previous note about 
addition and equivalence). □ 

Some interesting conjectures 

Conjecture 3.4.1. Tensor product gives these two maps: 

. D' bEA {M) x DTea(M) -J- &fi£ A {M). 

. r% A (M) x f b EA (M) -► r%+ A \M). 

Conjecture 3.4.2. Take some B£- G r* , +rS , ,(M). Take an element T**- G 
r|(M), such that V<S C U„S n , Vfi C ft, VCh k (£l C h) G At (5) it holds that Vw G 
C£*.(M) the elements of the multi-index matrix T£~(Ch k ) -u> remain to be from 
the class C^AM) C U„C£ (J (M). Then it holds that Bf;;T£;; is equivalent to 
an element of rjj^ o) (M)" (For subclasses 1%^ o) (M) and ^(M) it is 
equivalent to an element £) /m + a - ^(M).) The element is on its domain given 
as mapping 

(u,Ch k )^BZ;;;(T£;;u). (3.55) 

Some additional theory 

Theorem 3.4.1. Any arbitrary T£;; G A (as defined by 13 .4.3|> defines an A 8 {T£;;) 
object on M. Take any arbitrary S from the domain of T£\" , any arbitrary S C 
At{S) fl 5, any arbitrary f^c/i and any arbitrary Ch k (Vt C h) G 5. Then for fi c/l \ 
VL(Ch k ) it holds that multi-index matrix of functions T^—(Ch k ) can be obtained 
from the tensor components of A s (Tjf-) in Ch k (Vt C h) by the inverse mapping to 
Ch k (Q Ch ). 
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Proof. For MVtch, take fixed but arbitrary Ch k (Q C h) £ S and take T£-(Ch k ). 
Then Vg G fi^fo, \/Cti{n Ch ,q) G 5, and Vw £ 6 A n (^, C7i'(g, fi c/l )), we see that 
u'(Ch') is a delta-sequence. That means we just have to show, that on the set 
where T^-(Ch k ) is continuous in the 13.4.31 sense, the delta-sequencies give the 
value of this multi-index matrix. So write the integral: 



ich'(n ch 
By substitution x = e.z we obtain: 



/ T£;[{Ch k ){x)] i u/(-) ci 4 *. 

Jch'<n nh ) e ve/ 



ici 



7?- [(C7i fc ) (e.z)] w'(z) d 4 ^. 

'Ch'(Sl h) 

But from the properties of Tjf- [(C7ifc)(x)] it follows that 
/ (T£- [(Cfc) (z ) - K;y)} u/(z) d A z 

JCh'iClch) 



< [ Tf"[(Ch k )(zo + m)] u\z)d A z 
Jch'(n Ch ) 

< [ (Tt: [(Ch k )(z )\ + K£;;e) u'(z) d A z, 
Jch'(n nh ) 



>ch'(n Ch ) 
for some e small enough. 

But we are taking the limit e — \ which, considering the fact that w(x) are 
normed to 1, means that the integral must give Tjf- [(C7ifc)(z )]. The set, where it 
is not continuous in the sense of 13.4.31 has Lebesgue measure 0. That means the 
Vt C h part, which is mapped to this set has Lebesgue measure 0. But then the values 
of the multi-index matrix in the given chart at this arbitrary, but fixed point give 
us an associated field (and are independent on delta sequence obviously). □ 

Theorem 3.4.2. The field associated to a T£; G A n D'™ E {M), transforms for 
each Vt C h, for every pair of charts from its domain, Chi(Qch), Ch 2 (VL C h) on 
some M/(£l(Chx) U Vt{Ch2)), as an ordinary tensor field with piecewise smooth 
transformations 23 



23 Of course, some transformations in a generalized sense might be defined also on the Q(Chi) U 
n(Ch 2 ) set. 
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Proof. All this immediately follows from what was done in the previous proof, 
and from the fact that union of sets with Lebesgue measure has Lebesgue mea- 
sure 0. □ 

Note, that if there exists such point that for the object T™{M) we have 

]JmTf"(u e ) = ±oo 

at that point, then the field associated to this object, can be associated to another 
object, which is nonequivalent to this object. This means that the same field can 
be associated to mutually non-equivalent elements of F™(M). This is explicitly 
shown and proven by the next example. 

Theorem 3.4.3. Take 5(q,Ch k (£l C h)) £ D', So JM) being defined as mapping from 
each 4-form C^,*. (M) (S C S) to the value of this form's density at the point q in 
the chart Ch k {Vt C h) £ S, (q G Qch)- Then any power n G N + of 5(q, Ch k {Vt C h)) is 
associated to the function being defined on the domain M \ {q} and everywhere 0. 
Note that this function is associated to any power (n G N + ) (being a nonzero nat- 
ural number) of 5(q), but different powers of 5(q) are mutually nonequivalenQ 

Proof. Contracting powers of S(q, Ch k {Vt C h)) with a sequence of 4-forms from ar- 
bitrary A n (Ch'(q',£l' Ch ),S), (q' ^ q) (they have a support converging to another 
point than q) will give 0. For q = q' (|3.43|) gives 

lim e- 4 a/[(C7ufc)(0)] = ±oo . 

Now explore the equivalence between different powers of S(Ch k , q). 

5 n (q, Ch k {Vt C h)) applied to u € (x) G A n (Ch' k (q, Qch),S) will lead to the expression 

e" 4n w n (-). Then if we want to compute 



24 It is hard to find in our theory a more "natural" definition generalizing the concept of delta 
function from M. n . But there is still another natural generalization: it is an object from r9 y s o JM), 
defined as: S(Ch k (fl C h),q,^) = u' [(Ch k )(q)], Ch k (Q C h) € S n , u £ C^gJClch) ($n C S n ), qis 
image of q given by the chart mapping Ch k (Qch)- So it gives value of the density u' in the chart 
Ch k (Clch)/ at the chart image of the point q. 
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,n /' X \ l , ,m / X \ \ i^f \ 74, 



^i c , M U w "W-^ wm l7"* (x)rfl 



it leads to 



Sa ^ $ (°) / K« - ^y ^ m w) ** 

e ^° e Jch'(n Ch ) V e "■ ; / 

which is for n ^ m, n,m E N + clearly divergent, hence nonzero. D 

Note that despite of the fact that within our algebras we, naturally, have all the 
n E N + powers of the delta distribution, they are for n > 1, unfortunately, not 
equivalent to any distribution. 



Theorem 3.4.4. We see that the map A s is linear (in the sense analogous to 13.4.11) , 
and for arbitrary number of gfr",...,h£" 6 An P, ,= . AM) one has: Take 

\/VL Ch , \/n, \/S c S n nAt(S), \/Ch k {VL Ch ) e S, uA(C7i fc ) to be the union of all 
&i(Ch k ) related to the objects g£y, ..., h%'\ Then 

A s {g%: ® ... <8> K-) = A s {g% : ) ® ... ® A S {K;:) 

on Vl C h \ Ui£li(Chk). Here the first term is a product between A objects and the 
second is the classical tensor product. 

Proof. It is trivially connected with previous proofs: Note that from the definition 
(|3.4.2|) for appropriate 4-form fields co e we have 

A 9 (jgZ~®...®h£")(Ch k )=Km gt"[Ch kl u t )...K"(Ch h ,u 6 ). (3.56) 



But for the objects #?;;;, ..., h%\\\ E A, with respect to the theorem (|3.4.1[) neces- 
sarily 

lim g%"iCh k ,u e )...h£~(Ch k ,u e ) = A s (gf~) ® ... ® A S (K~)- (3-57) 

This proves the theorem. □ 
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This means that the result of tensor multiplication of elements from A (it has 
product of two scalars as a subcase) is always equivalent to some element from A. 
This is a result closely related to the theorem (|3.4.3|) . It tells us that multiplication 
is a mapping between equivalence classes formed of more constrained classes as 



those mentioned in the part a) of the theorem (|3.4.3|) 



Conjecture 3.4.1. If T£- G D'^AM) has an associated field, then it transforms 
on its domains as a tensor field. 

Conjecture 3.4.2. If T£;; G T^ o) (M) has an associated field and Lf~ G T n s {M), 
then A a (T£;;®L%-) = A s {T£;;)®A s {Lf;;). (The <g> signhas again slightly different 
meaning on the different sides of the equation). 

Conjecture 3.4.3. Take C£;;,Dfc,Ffc;;,B%™ G T^(M), such that they belong to 
the same classes TJ^ - (M). Also assume that V«S, such that C£;;,Dfc;,F$;;;, 
Bf- G F™ S {M), and VS C S, each of the elements C£;;,D£;;,Ffc"; Bf- has 
for every At(S) associated fields defined on the whole M. Then F2- w Bo- and 

C »... „ D ,... impHes {C ^ F) M..J : „ {D g, jB) M...a..._ 

3.4.5 Covariant derivative 

The last missing fundamental concept is the covariant derivative operator on 
generalized tensor fields (GTF). This operator is necessary to formulate an appro- 
priate language for physics and generalize physical laws. Such an operator must 
obviously reproduce our concept of the covariant derivative on the smooth tensor 
fields (through the given association relation to the smooth manifold). This is pro- 
vided in the following section. The beginning of the first part is again devoted to 
fundamental definitions. The beginning of the second part gives us fundamental 
theorems, again just generalizing Colombeau results for our case. After these the- 
orems we, (similarly to previous section), formulate conjectures representing the 
very important and natural extensions of our results (bringing a lot of new sig- 
nificance to our results). The last subsection in the second part being again called 
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"some additional theory" brings (analogously to previous section) just physical 
insight to our abstract calculus and is of lower mathematical importance. 

Definition of (^-derivative and connection coefficients 

Definition 3.4.1. We define a map, called the <9-derivative, given by smooth vector 
field U l (smooth in the atlas S) as a mapping TV^JM) — > TT^JM), given on its 
domain VT2c/i and Ch k (Vt C h) as: 

T^; {u) {Ch k ,u) = -T£;;(Ch k , (U a co), a ) u e C p (n Ch ). 

Here (U a tu) ta is understood in the following way: We express U a in the chart 
Ch k (Q C h) and take the derivatives (U a (Ch k ) u'(Ch k )) a in the same charH 
Ch k (Vt C h)- They give us some function in Ch k (Qch), which can be (in this chart) 
taken as expression for density of some object from C p {Vt C h)- This means we 
trivially extend it to M by taking it to be everywhere outside Vt C h- This is the 
object used as an argument in T£;\ 

To make a consistency check: This T^"',™ is an object which is defined at 
least on the domain C P [S AM) for 5C5 (S related to U l in the sense that U l is 

0(0) 

smooth in S), and is continuous on the same domain. This means it belongs to 
the class T n (M). To show this take some arbitrary Qch and some arbitrary chart 
Ch k (Vt C h) G S. We see that within Ch k (Vtch) the expression (U a u'), a is smooth 
and describes 4-forms, which are compactly supported, with their support being 
subset of Vtch- Hence they are from the domain of T£\" in every chart from S. 
In any arbitrary chart from S we trivially observe, (from the theory of distribu- 
tions), that if u n — > u, than ^'" ([/ ,(w„) — >■ T^'" ^Auo). It means that T^'^Joj) is 
continuous. 

Definition 3.4.2. Now, by generalized connection we denote an object from T 3 (M) 
such that: 



25 We will further express that the derivative is taken in Ch k (^lch) by using the notation 

(U a u'(Ch k )) A{Chk)a] . 
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• The set 

U nch Us I (w, C7i(Q c/l )) : w 6 C^(fi C fc), C7i(fi C /») G ^l 
belongs to its domain. 



CTiJ 



• It is V«S, MVtch, yCh k (VL C h) e v4 continuous on Cg.gJQch) with Ch k {Vt 
taken as its argument. 

• It transforms as: 

T%{Ch^e>) = I^(C/n, ((J- 1 )^^ 1 )*^ - -ftOr 1 )^) w). (3.58) 

Definition of covariant derivative 

Definition 3.4.3. By a covariant derivative (on Q,ch) hi the direction of a vector 
field J7*(M), smooth with respect to atlas S, of an object T<f- e 1^(0^), we 
mean: 

D ciu) TZ:iu) = TZ:; t(u) (u) + r^(u)T£-(ir u ) - r: p (u)Tt-(u^). (3.59) 

This definition (|3.59[) automatically defines covariant derivative everywhere 
on VTgJM). This can be easily observed: The <9-derivative still gives us an object 
from TVgJM), and the second term containing generalized connection is from 

T™ } (M) trivially too. 

Definition 3.4.4. Furthermore, let us extend the definition of covariant derivative 
to the class T^ A (M) just by stating that on every nonlinear object (note that every 
such object is constructed by tensor product of linear objects) it is defined by the 
Leibniz rule. (This means it is a standard derivative operator, since it is trivially 
linear as well.) 
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The S' n class 

Notation. Take as V n some n-times continuously differentiable subatlas of A. 
Then the class S' n related to the atlas V n is formed by objects T£;; G r® .= . . (M), 

such that: 

• U m S m = V n and Vm, At(S m ) = V n/ 

• it is given \/VL Ch , \/Ch k {Vt Ch ) G V n , Mm, Moo G C* § A£l C h) as a map 

(oo,Ch k )^ J T»;;;{Ch k )u. 

Here Tj};(Ch k ) is a multi-index matrix of n-times continuously differen- 
tiable functions (if it is being expressed in any arbitrary chart from V n ). 

Basic equivalence relations related to differentiation and some of the interest- 
ing conjectures 

Theorem 3.4.1. The following statements hold: 

a) Take a vector field U\ which is smooth at S C V n+ i c V n for n > 1, (for- 
mally including also n = n + 1 = oo, hence 5), plus a generalized connection 
T» a G r|(M), and T£ ;■ G £>£ a (M) n S n- s 'n ^ here related to the given V n . 
Then it follows that: D C {u)T^;;; is an object from S' n+1 (being related to the 
given V n+1 ). Moreover, the equivalence class S' n+1 of the image contains 
exactly one linear element given for every chart from its domain as integral 
from some multi-index matrix of piecewise continuous functions. Particu- 
larly this element is for \fil C h and arbitrary such oo G C p (Vt C h), Ch k {Vt C h) 
that are from its domain a map: 

(Ch k , oo) -)■ f U a T^ a (Ch k ) uo. (3.60) 

(Here ";" means the classical covariant derivative related to the "classical" 
connection, components of which are in Ch k given by T^ s (Ch k ), and the 
tensor field, components of which are in Ch k given by Tj}-(Ch k ).) 
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b) Take U l being smooth in S and T^ a e r|(M). Then the following holds: 
The covariant derivative is a map D'^LJM) -» f7£? A (M) and the classes 
f™+™ (M) of the image contain (exactly) one element of D'^ So) (M). 

Proof. a) Covariant derivative is, on its domain, given as 

(u, Ch k ) -> l£" >(I0 (u/) + r£,(u/)!Z£;;(r/»u;) - r?»T£-(tf>u;). 

Take the first term T^"', V y Express it on Vt C h in Ch k (Vt C h) £ "Cn as the map: 

w-» / T^-(Ch k )co'(Ch k )d 4 x. 

Jch k (U C h) 

So analogously: 

T^- (Ch k ,u) = - [ T£;;{Ch k ) (U"(Ch k ) u'(Ch k )) i[{Chk)a] d 4 x. 

Jch k (n Ch ) 

Here to' is in arbitrary chart from V n , being n-times continuously differen- 
tiable, (the domain is limited to such objects by the second covariant deriva- 
tive term, which is added to the (^-derivative), and such that the expression 
(U % (Ch k ) u>'(Chk)),[(Ch k )i] is in any Ch k (Vt C h) € V n n-times continuously dif- 
ferentiable. 

Now by using integration by parts, (since all the objects under the integral 
are at least continuously differentiable (n > 1), it can safely be used), and 
considering the compactness of support we obtain: 

T H:::,(u)( Ch k> u >) 

T£;;{Ch k ) (U a (Ch k ) u'(Ch k )) i[{Chk)a] d 4 x 

ch k (n Ch ) 

Tt.iCh k ) t[{ ch k)a ]U a (Ch k ) uj'(Ch k ) d 4 x. (3.61) 

Ch k (n Ch ) 
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Then it holds that: 

T^:(Ch m ) A{Chrn)a] (U a (Ch m ) u/(Chm) d A x 

Ch m (flch) 

{J^(J-y v ...Tl:{Ch k )) j[{ c hh) a}U a (Ch k ) u'(Ch k ) d A x. (3.62) 

Ch k (Q oh ) 

We see that this is defined and continuous in any arbitrary chart from V n+ i 
for every C^.^AM), such that 3S C S and S C T> n+1 . 

Now we see that the second term in the covariant derivative expression is 
equivalent to the map (with the same domain): 

(Ch k ,u) -». J(T» p U a TP;;; - T a up U»T^)(Ch k ) u , 

(see theorem |3.4.3|) , and between charts the objects appearing inside the in- 
tegral transform exactly as their classical analogues. This must hold, since 
T£~ is everywhere continuous (in every chart considered), hence on every 
compact set bounded, so the given object is well defined. This means that 
when we fix this object in chart Ch m {Vt C h), and express it through the chart 
Ch k (Q C h) and Jacobians (with the integral expressed at chart Ch k (Vt C h) ), as 
in the previous case, we discover (exactly as in the classical case), that the 
resulting object under the integral transforms as some object D'^(M), with 
the classical expression for the covariant derivative of a tensor field appear- 
ing under the integral. 

b) The resulting object is defined particularly only on C^.^JM) (S C S). We 
have to realize that T£" can be written as a (N — > oo) weak limit (in every 
chart from S) of T^' N e D'^, So JM). It is an immediate result of previous 
constructions and Colombeau theory, that 

- r?p(0 »£/'.) « -T£7(ry7'.) (3.63) 



CHAPTER 3. MULTIPLICATION OF TENSORIAL DISTRIBUTIONS 155 

Now take Vflch both, Tff'",^ and (|3.63[) fixed in arbitrary Ch k (Vt C h) € 
5. Write both of those objects as limits of integrals of some sequence of 
"smooth" objects in Ch k (Vt C h)- But now we can again use for T^'", V s an 
integration per parts and from the "old" tensorial relations; we get the "ten- 
sorial" transformation properties under the limit. This means that the re- 
sulting object, which is a limit of those objects transforms in the way the 
D'™(M) objects transform. 

□ 



Theorem 3.4.2. Part a) of the theorem (|3.4.1[) can be also formulated through a 



generalized concept of covariant derivative, where we do not require the U % vector 
field to be smooth at some S, but it is enough if it is n + 1 differentiable in V n+ i. 

Proof. We just have to follow our proof and realize that the only reason why we 
used smoothness of U % in S was that it is required by our definition of covariant 
derivative (for another good reasons related to different cases). □ 

This statement has a crucial importance, since it shows that not only all the 
classical calculus of smooth tensor fields with all the basic operations is contained 
in our language (if we take the equivalence instead of equality being the crucial 
part of our theory), but it can be even extended to arbitrary objects from S' n . (If 
the covariant derivative is obtained through connection from the class T%(M).) In 
other words it is more general than the classical tensor calculus. 

We can now think about conjectures extending our results in a very important 
way: 

Conjecture 3.4.1. Take an arbitrary piecewise smooth, and on every compact set 
bounded @ vector field U\ Take also T^ a e T|(M) and such T£;- e D'™ E (M), 



26 To be exact, the expression "covariant derivative" is used in this and the following conjecture 
in a more general way, since we do not put on U l the condition of being smooth in some subatlas 
S. 
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that VQ C h, WCh k (Q C h) e A 3Chi(n C h) e A, inwhidJH 



TZ:iCh k )r% s (ch k ) d A x 

converges on every compact set Q' C Qch- Then the following holds: The covari- 
ant derivative (along U % ) maps this object to an element of some equivalence class 
from T^ +n (Af ). This class contains (exactly) one element from D'™(M). 

Conjecture 3.4.2. Take U l being a piecewise smooth vector field, and T^ a e 
T|(M). Then the following holds: The covariant derivative along this vector field 
isamap: D^ Si o) (M) -> fjg o)A {M), and the classes f™^ o)A (M) of the image 
contain (exactly) one element of D'™(M). 

Theorem 3.4.3. For U l being a smooth tensor field in S with the connection taken 
from r|(M), T£" e r™ o)A (M) and B£y e r™ o) (M), it holds that T£;- « ££•- im- 
plies D%u\T£" « D^jjsBfr;; for arbitrary natural number n. 

Proof. Pick an arbitrary Vl C h and an arbitrary fixed chart Ch'(Vt C h) €E 5. Such a 
chart maps all the 4-forms from the domain of ^T So)A {M) objects to smooth com- 
pact supported functions (given by densities expressed in that chart). The objects 
T£;;((U a u) tCt ) and Bfr;;((U a oj) :a ) are taken as objects of the Colombeau algebra 
(the connection, fixed in that chart, is also an object of the Colombeau algebra) 
and are equivalent to U a (u))T£\\\ a (u)) and U a (u))B£;; a (u;). Here the derivative 
means the "distributional derivative" as used in the Colombeau theory (fulfilling 
the Leibniz rule) and U a (u) is simply a D'£g, So JM) object with the given vec- 
tor field appearing under the integral. But in the Colombeau theory one knows 
that if some object is equivalent to a distributional object, then their derivatives 
of arbitrary degree are also equivalent. It also holds that if any arbitrary object is 
equivalent to a distributional object, then they remain equivalent after being mul- 
tiplied by arbitrary smooth distribution. In the fixed chart we have (still in the 
Colombeau theory sense), 



27 This means we are trivially integrating T^ ;T" on compact sets within subset of R 4 given as 



■06 

image of the given chart mapping. 
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ZA.. ~ V... ' 

But since their ^-derivatives were, in the fixed chart, obtained only by the 
distributional derivatives and multiplication by a smooth function, also their d- 
derivatives must remain equivalent. The same holds about the second covariant 
derivative term (containing connection). So the objects from classical Colombeau 
theory, (classical theory just trivially extended to what we call multi-index matri- 
ces of functions), obtained by the chart mapping of the covariant derivatives of 
T^;; and B£;;, are equivalent in the sense of the Colombeau theory. But the Vt C h 
set was arbitrary and also the chart was an arbitrary chart from the domain of 
T^; , B%"\ So T£- and B£;; are equivalent with respect to our definition. □ 

We can try to extend this statement to a conjecture: 

Conjecture 3.4.3. Take U l being piecewise smooth tensor field, take connection 
from the class r|(M), T£;; G T^ lSl o)A (M), and B£;; G T$ iSi o) (M). Thenit holds 
that T£" ps Bfr" implies D^^T^- ps D^^B^- for arbitrary natural number n, 
if such covariant derivative exists. 

This conjecture in fact means that if we have connection from the class r|(M), 
then the covariant derivative is a map from such element of the class T!™ s \ A (M), 
that it contains some linear element, to T™(M). Note that we can also try to prove 
an extended version of the conjecture, taking the same statement and just extend- 
ing the classes r™^ o)A {M), r^ o) (M) to the classes T? UlSl)A (M), T™ 5;) (M). 

Some additional theory 

Theorem 3.4.4. Take U l to be vector field smooth in some S C V n , with T^ a G 
A and T£;- G S' n n D' b a (M) (n > 1, S' n is related to V n ). Then D C {u)T£;; has 
an associated field which is on M\ Vt{Ch) the classical covariant derivative of 
A S (T£;;). (Q(Ch) is a set on which is T^ a (Ch) continuous and is of Lebesgue 
measure.) It is defined on the whole V n+1 (S C V n+1 ). That means association 
and covariant differentiation in this case commute. 
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Proof. Just take the definition of the classical covariant derivative, and define the 
linear mapping given VT^ and arbitrary Ch k (Vt C h) £ 2\ as 

(w, Cfc*) -> / U v (Ch k ) [A S {TZ : %] (Ch k ) u. (3.64) 

This is an internally consistent definition, since [A S {T^;;;).^\ (Chk) is defined 
everywhere apart of a set having L measure 0. Now from our previous results 
follows that everywhere outside Vt(Ch) 

[A s (Tt;)J(Ch k )=Tt;., ll (Ch k ), 

and so the linear mapping (|3.64|) is equivalent to the object Dc(u)Tj?;;;. Then 
U v A a {T^")., l = A 3 (D C{U) T£;;) everywhere outside the set Q{Ch). □ 



Theorem 3.4.5. An extended analogy of theorem (|3.4.4[) , can be proven, if we use 



generalized concept of covariant derivative, without assuming that vector field is 
smooth in some S, but only n + 1 continuously differentiable within V n+ i. 

Proof. Exactly the same as before. □ 

This means that the aim to define a concept of covariant derivative, "lifted" 
from the smooth manifold and smooth tensor algebra to GTF in sense of associa- 
tion, has been achieved. It completes the required connection with the old tensor 
calculus. 

Conjecture 3.4.4. Take IP smooth in S, T£;; £D'™ {So) (M), such that 3A S (T£;;), 
T L £ r l( M )- Then 3 A S (D C{U) T£;;) and holds that 

A S (D C(U) TZ:;) = U a A s {T^;;), a . 

(Hence, similarly to ®, the covariant derivative operator commutes with associa- 
tion for some significant number of objects.) 

Note that for every class T™ t , S j (M) we can easily define the operator of Lie 
derivative along arbitrary in S smooth vector field V (not along the generalized 
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vector field, but even in the case of covariant derivative we did not prove anything 
about larger classes of vector fields than smooth vector fields). Lie derivative can 
be defined as (L V T, u) = (T, L v oo). This is because the Lie derivative preserves 



ip 

jiracivra uic j^iupci lien ui auv.ii <„ 

holds that S C S. 



n-forms and also preserves the properties of such C;~. (M) classes, for which it 



3.4.6 Basic discussion of previous results and open questions 

We have constructed the algebra of GTFs, being able to incorporate the concept of 
covariant derivative, (with the given conditions on vector fields and connection), 
for a set of algebras constructed from specific distributional objects. The use of 
these ideas in physics is meaningful where the operations of tensor product and 
covariant derivative give a map from appropriate subclass of D'™(M) class to the 
elements of T m (M) containing a D'™(M) element. This is always guaranteed 
to work between appropriate subclasses of piecewise continuous distributional 
objects, but a given physical equivalence might specify a larger set of objects for 
which these operations provide such mapping. Note that the whole problem lies 
in the multiplication of distributions outside the D'^(M) class. (For instance it 
can be easily seen that square of 5(Chk, q) as introduced before is not equivalent 
to any distribution.) This is because the product does not have to be necessarily 
equivalent to a distributional object. Even worse, in case it is not equivalent to 
a distributional object, the product is not necessarily a mapping between equiv- 
alence classes of the given algebra elements (D'£\(M)). The same holds about 
contraction. 

But even in such cases there can be a further hope. For example we can aban- 
don the requirement that certain quantities must be linear, (for example the con- 
nection), and only some results of their multiplication are really physical (mean- 
ing linear). Then it is a question whether they should be constructed (constructed 
from the linear objects as for example metric connection from the metric tensor) 
through the exact equality or only through the equivalence. If we take only the 
weaker (equivalence) condition, then there is a vast number of objects we can 
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choose, and many other important questions can be posed. Even in the case that 
the mathematical operations do depend on particular representatives of the equiv- 
alence classes, there is no necessity to give up; in such situation it might be an 
interesting question if there are any specific "paths" which can be used to solve 
the physical equivalence relations. The other point is that if these operations do 
depend on the class members, then we can reverse this process. It means that 
for example in the case of multiplication of two delta functions we can find their 
nonlinear equivalents first and then take their square, thus obtaining possibly an 
object belonging to an equivalence class of a distribution. 

As I mentioned in the introduction, these are not attempts to deal with phys- 
ical problems in a random, ad hoc way. Rather I want to give the following in- 
terpretation to what is happening: The differential equations in physics should 
be changed into equivalence relations. For that reason they have plentiful solu- 
tions in the given algebrao. (By a "solution" one here means any object fulfilling 
the given relations.) One obtains much "more" solutions than in the case of clas- 
sical partial differential equations (but all the smooth distributions representing 
"classical" solutions of the "classical" initial value problem are there), but what is 
under question is the possibility to formulate the initial value problem for larger 
classes of objects than D'^ E (M) and D'^ So JM) (see the next section). Moreover, if 
this is possible then there remains another question about the physical meaning 
of those solutions. It means that even in the case we get nonlinear objects as so- 
lutions of some general initial value problem formulation, this does not have to 
be necessarily something surprising; the case where physical laws are solved also 
by physically meaningless solutions is nothing new. The set of objects where we 
can typically search for physically meaningful solutions is defined by most of the 
distributional mappings (that is why classical calculus is so successful), but they 
do not have to be necessarily the only ones. 



28 After one for example proves that covariant derivative is a well defined operator on all GTF 
elements, then it is the whole GTF algebra. 



CHAPTER 3. MULTIPLICATION OF TENSORIAL DISTRIBUTIONS 161 

3.4.7 Some notes on the initial value problem within the partial 
differential equivalence relations («) on D'™(M) 

In this section we will suggest how to complete the mathematical structure de- 
veloped and will get some idea how a physical problem can be formulated in our 
language. It is again divided into what we call "basic ideas" and "some additional 
ideas". The first part is of a considerable importance, the second part is less im- 
portant, it just gives a suggestion how to recover the classical geometric concept 
of geodesies in our theory. 

The basic ideas 

The approach giving the definition of the initial value problem Take a hyper- 
surface (this can be obviously generalized to any submanifold of lower dimen- 
sion) iV c M, which is such that it gives in some subatlas An C A a piecewise 
smooth submanifold. In the same time An is such atlas that 3S S C An- 

If we consider space of 3-form fields living on N (we give up on the idea of 
relating them to 4-forms on M), we get two types of important maps: 

• Take such D'^(M) objects, that they have in every chart from An asso- 
ciated (tensor) fields defined everywhere on N, apart from a set having 3 
dimensional Lebesgue measure equal to 0. Such objects can be, in every 
smooth subatlas S C Am, mapped to the class D'^(N) by embedding their 
associated tensor fields^] into N. This defines a tensor field T£\" living on a 
piecewise smooth manifold N. Furthermore if A^d is some largest piece- 
wise smooth atlas on N, the tensor field Tff~ defines on its domain a map 



Vn Ch C N, Ch k (n Ch ) e A 3D and LueC F (n 



Ch 



{uj,Ch k )^ / T^;(Ch k )u. (3.65) 

'n Ch 



29 As previously noted, the given T{f £ D'^(M) we can define in every chart by (Chk,u>) 

j nch [A s {TZ)]{Ch k )u. 
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This is object from the class D'^(N). What remains to be proven is that in 
any smooth subatlas we map the same D'^ E (M) object to D'™(N), otherwise 
this formulation is meaningless. 

• The other case is a map D'^ So) (M) ->■ D'™ {So) (N) (S c An), defined in a 
simple way: The objects from D%g, So JM) are mapped as associated smooth 
tensor fields (in the previous sense). After this step is taken, one maps the 
rest of distributional objects from D'£? So) (M) by using the fact that they are 
weak limits of smooth distributions D'" l s , s AM). (This is coordinate inde- 
pendent for arbitrary tensor distributions.) So in the case of objects outside 
the class D'^,, s , (M) we embed the smooth distributions first, and take the 
limit afterwards (exchanging the order of operations). The basic conjecture 
is that if this limit exists on M, it will exist on N (in the weak topology), by 
using the embedded smooth distributions. 

Now we can say that initial value conditions of, (for example), second-order 
partial differential equations are given by two distributional objects from 
D%So)( N i)> D T { So){ N 2) (DTe(N 1 ) / D'™ e {N 2 )) ontwohypersurfaces N U N 2 (not 
intersecting each other). The solution is a distributional object from the same class, 
which fulfills the ~ equation and is mapped (by the maps introduced in this 
section) to these two distributional objects. 

Useful conjecture related to our approach Note, that we can possibly (if the 
limit commutes) extend this "initial value" approach through the D'^(M) class 
to all the weak topology limits of the sequences formed by the objects from this 
class. This means extension to the class of objects belonging to D'™(M), such that 
for any chart from A they have the full C P (M) domain and D'^(M) is dense in 
this class. 

Some additional ideas 

"Null geodesic solution" conjecture Let us conjecture the following: 



CHAPTER 3. MULTIPLICATION OF TENSORIAL DISTRIBUTIONS 163 

Conjecture 3.4.1. Pick some atlas S. Pick some g^ u G D^ s , So JM), such that it has 
A s (g^u), being a Lorentzian signature metric tensor field. Take some VL C h and two 
spacelike hypersurfaces H 1 , H 2l H x n H 2 = {0}, H x n Vt C h ^ {0}, # 2 n fi C h 7^ 
{0}. Furthermore Hi, H 2 are such that there exist two points q\ G Hi fl Slchr 
q 2 E H 2 (1 Vt C h separated by a null curve geodesies (relatively to A s (g^)), and 
the geodesies lies within Vt C h- Construct such chart Ch k (Vt C h) G S, that both of 
the hypersurfaces are hypersurfaces (they are smooth manifolds relatively to S) 
given by u = const, condition (u is one of the coordinates) and the given geodesies 
is representing w-coordinate curve. 

Take classical free field equation with equivalence: D g $ w (with g^ u as 
previously defined). Then look for the distributional solution of this equation with 
the initial value conditions being 5{Chi k (Vt C h H Hi), qi) G D'£f So \{Hi) on the first 
hypersurf ace and 5(Ch 2k (£lch^H 2 ),q 2 ) G D'^ So AH 2 ) on the second hypersurface. 
(Here Chi k (Vt C h H Hi), Ch 2k (£l C h H H 2 ) are coordinate charts, which are the same 
on the intersection of the given hypersurface and VL C h as the original coordinates 
without u.) Then the solution of this initial value problem is a mapping $, which 
is such, that it can be in chart Ch k (Vt C h) expressed as 



du I \\dx l (5{x l { qi ))u'(Ch k ){u,x 3 )) 



UJ -¥ 



(where x l (qi) is image of qi in chart mapping Ch k (VL C h) )■ 

We can formulate similar conjectures for timelike and spacelike geodesies, we 
just have to: 

• instead of point separation by null curve, consider the separation by timelike 
or spacelike curve, 

• instead of the "massless" equation we have to solve the (D 9 ± m 2 )$ fa 
equation (m being arbitrary nonzero real number). Here ± depends on the 
signature we use and on whether we look for timelike or spacelike geodesies. 



The rest of the conditions are unchanged (see |3.4.1[) . Some insight to our con 



jectures can be brought by calculating the massless case for flat Minkowski space, 
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using modified cartesian coordinates (u = x — ct, x, y, z). We get the expected re- 
sults. 



3.5 What previous results can be recovered, and how? 

As was already mentioned, our approach is in some sense a generalization of 
Colombeau approach from ||46|, which is equivalent to canonical M n approach. 
So for Vl C h, after we pick some Ch'(Vt C h) £ S, (which determines the classes 
A n (M) related to this chart), and by considering only the objects C^ ( ,AVt C h) (S C 
S), (hence considering the D'^ So \ (M) class only), we obtain from our construction 
the mathematical language used in [46J. But all the basic equivalence relations 
from Colombeau approach have been generalized first to the class D'^ S(>)A (M) 
and also to appropriate subclasses of the D'^ A (M) class. 

3.5.1 Generalization of some particular statements 

Now there are certain statements in R", where one has to check whether they are 
just a result of this specific reduction, or not. A good example is a statement 

H n 5 « — 5 . (3.66) 

n + 1 

(H is Heaviside distribution.) What we have to do is to interpret the symbols 
inside this equation geometrically. This is a M 1 relation. H is understood as a 
D', S AM) element and defined on the manifold (one dimensional, so the geom- 
etry would be quite trivial) by integral given by a function (on M ) obtained by 
the inverse coordinate mapping substituted to H. Now take some fixed chart 
C/i^M 1 ). The derivative is a covariant derivative along the smooth vector field 
U, which is constant and unit in the fixed coordinates Chk(M}). Then 5 can be 
reinterpreted as 5(Ch k , q), where q is the point in the chart C^IR 1 ). Then the 
relation can be generalized, since it is obvious that (see the covariant derivative 
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section) D C (u)H = 5(q, Ch k ) and so 






Dc(U)H = D C (u)L(Hf) 






= D C (u)L(HJ ) p 


* D c{u) {H n+1 ) - 


= (n + l)H n D c{u) H. 



(3.67) 

By L we mean here a regular distribution defined by the function in the brackets, 
hence an object from D' E (M). (To be precise and to avoid confusion in the nota- 
tion, we used for the Heaviside function the symbol Hf, while for the Heaviside 
distribution the usual symbol H.) This is nice, but rather trivial illustration. 

This can be generalized to more nontrivial cases. Take the flat R" topological 
manifold. Fix such chart Ch^iW 1 ) covering the whole manifold, that we can 
express Heaviside distribution in this chart through Hf{x\) G D', So AM). This 
means the hypersurface where Hf is discontinuous is given in Ch k (W l ) asx\ = 0. 
Now the derivative will be a covariant derivative taken along a smooth vector 
field being perpendicular (relatively to the flat space metric^) to the hypersurface 
on which is Hf discontinuous. We easily see that the covariant derivative of H 
along such vector field gives a distribution (call it 5 E D', So AM)), which is in the 
chart Ch k (R. n ) expressed as 

5(4>) = S X1 ( / (f)(x 1 , x 2 ... x n ) dx 2 ... dx n J . (3.68) 

This distribution reminds us in some sense the "geodesic" distribution from the 
previous part. Then the following holds: 

H n 5 w -^— 5 . (3.69) 

n + 1 

This generalized form of our previous statement can be used for computations 
with Heaviside functional metrics (computation of connection in fixed coordi- 
nates). 



30 Note that since it is flat space it makes sense to speak about a perpendicular vector field, since 
we can uniquely transport vectors to the hypersurface. 
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3.5.2 Relation to practical computations 

These considerations (for example) imply that the result from canonical Colom- 
beau IR n theory derived by [|46| can be derived in our formalism as well. This is 



also true for the geodesic computation in curved space geometry from |104[ . The 
results derived in special Colombeau algebras (in geometrically nontrivial cases) 
are more complicated, since in such cases the strongest, v4°°(IR n ), version of the 
theory is used. This version is not contained in our chart representations. (This is 
because we are using only A n (R") with finite n.) It is clear that all the equivalence 
relations from our theory must hold in such stronger formulation (since obviously 
A°°(M) C A n (M) Vn E N) and the uniqueness of distribution solution must hold 
as well. This means that at this stage there seems to be no obstacle to reformulate 
our theory by using A°°(Ch k , q, Vt C h) classes (and taking elements from D'.AM) 
at least), if necessary. But it is unclear whether one can transfer all the calculations 
using A°°(Chk, q, &ch) classes to our weaker formulation. 

The strong formulation was used also in the Schwarzschild case |[76|, but there 
is a problem. The fact that the authors of ||76l regularize various functions piece 
by piece does not have to be necessarily a problem in Colombeau theorycJ- But, 
as already mentioned, the problem lies in the use of formula for R„ v , origi- 
nally derived within smooth tensor field algebra. If we want to derive in the 
Schwarzschild case Ricci tensor straight from its definition, we cannot avoid mul- 
tiplications of delta function by a non-smooth function. This is in Colombeau 
theory deeply non-trivial. 

In the cases of Kerr's geometry and conical spacetimes theory this problem 
appears as well. As a consequence of this fact, calculations are mollifier depen- 
dent, not being (in the strict sense) results of our theory anymore. On the other 
hand there is no reasonable mathematical theory in which these calculations make 
sense. This means that a better understanding of these results will be necessary. 
By better understanding of these results provided by our theory we mean their 



31 Although the authors use (in the first part, not being necessarily connected with the results) 
quite problematic embeddings. 
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derivation by a net of equivalence relations, by taking some intermediate quanti- 
ties to be nonlinear. So the results should follow from the principle that the equiv- 
alence relations are the fundamental part of all the mathematical formulation of 
physics. 

3.6 Conclusions 

The main objectives of this work were to build foundations of a mathematical 
language reproducing the old language of smooth tensor calculus and extending 
it at the same time. The reasons for these objectives were given at the beginning 
of this chapter. This work is a first step to such theory but it already achieves its 
basic goals. That means we consider these results as useful independently of how 
successful future work on the topic will turn out to be. On the other hand, the 
territory it opens for further exploration is in my opinion large and significant. It 
offers a large area of possibilities for future work. 

Just to summarize: the result of our work is a theory based purely on equiv- 
alence relations instead of equalities, using a well defined concept of generalized 
tensor field and the covariant derivative operator. This operator is well defined at 
least on the proper subclass of generalized tensor fields. We also defined (using 
some conjectures) the initial value problem for partial differential equivalence re- 
lations. Our theory naturally relates to many results beyond the classical smooth 
tensor calculus, already derived. 



Conclusions 



As we mentioned in the introduction, this thesis is based on three separate re- 
search projects. 

The first project was related to pseudo-Finsler extensions of the general theory 
of relativity, and represented an attempt to find a natural geometric framework 
for possible high energy Lorentz violations. The reason why one was interested 
in such construction is the question of whether it is possible to find a weaker in- 
terpretation of Einstein's equivalence principle consistent with Lorentz symmetry 
violations. The result obtained was, when mathematical simplicity was taken as 
the guidance principle, unfortunately a "no-go" theorem, at least for significant 
number of cases. The cases particularly affected were the bi-metric theories, but 
the analogue model based on bi-refringent crystal optics indicated that the prob- 
lem might affect much larger class of theories. The problem lies in the fact that, 
unlike what one would naturally expect, introducing Lorentzian signature puts 
very tight constraints on Finsler geometry, at least if one wants to keep some of 
the basic geometric concepts well defined and meaningful. This "no-go" result 
we consider to be disappointing, (which is the case of most "no-go" results), but 
certainly very useful. 

The second project focused on the highly damped quasi-normal modes of dif- 
ferent black hole spacetimes. The method of approximation by analytically solv- 
able potentials was used to estimate the highly damped modes for the Schwarz- 
schild and the Schwarzschild-de Sitter (S-dS) black holes. The first served more as 
a consistency check (since the asymptotic formula for the highly damped QNMs 
of the Schwarzschild black hole is well known). But for the Schwarzschild-de Sit- 

168 



ter black hole a lot of new information was extracted from those models, especially 
the link between rational ratios of horizon surface gravities to the periodic be- 
haviour of the QNMs. Also, when periodic, in general the highly damped modes 
do not form only one equi-spaced family as in the case of Schwarzschild black 
hole, but split into multiple families. Strikingly, as we discovered, the same pat- 
terns can be observed in a complementary set of analytic estimations for the highly 
damped modes, those approximations being obtained by monodromy techniques. 
This holds for all types of black hole spacetime so far analysed by those tech- 
niques. That means we were able to significantly generalize our theorems about 
the highly damped mode behaviour to all the presently known analytic results. 
Our results might be interesting also from the viewpoint of the black hole ther- 
modynamics, as the asymptotic QNM behaviour is suspected to be linked to the 
black hole area spectrum |f77[|1101 . 

The third project dealt with the problem of multiplication of tensorial distri- 
butions. Despite the fact that lot has been done in the field in the past ffTTl |173|| , 
the full generalization of the covariant derivative operator was, for example, not 
yet achieved. On the other hand practical results confirm a need for such gen- 
eralization (see 111731 ). We built an alternative construction, which fully operates 
with the Colombeau equivalence relation, but technically avoids Colombeau al- 
gebra construction. It generalizes the concept of covariant derivative to tensorial 
distributions and operates on the much more general, but for the language of 
distributions natural, piecewise smooth manifolds. We are convinced that such 
language might offer conceptual extension of general relativity and might have 
possibly interesting consequences for quantum gravity as well. 



169 



Appendix A 
Bi-refringent crystals 



A.l Basic characteristics of the crystal media 



The basic optics reference we shall use is Born and Wolf, Principles of Optics 11291 . In 
particular we shall focus on Chapter XV, "Optics of crystals", pages 790-818. See 
especially pages 796-798 and pages 808-811. Specific page, chapter, and section 
references below are to the 7th (expanded) edition, 1999/2003. 

The theory of bi-refringent crystal optics is formulated in the preferred inertial 
system, the inertial system of the crystal. The optical medium of the crystal is 
characterized by permeability and permittivity. Permeability fi is taken to be a 
scalar, permittivity e^ a "spatial", (relative to the inertial system of the crystal), 
3x3 tensor. (This is an excellent approximation for all known optically active 
media.) By going to the principal axes we can, without loss of generality, take e^ 
to be diagonal 

" e^ 
e.„ 



(A.l) 



e z 

This fixes the relativistic inertial coordinate system in a unique way. 
We furthermore define "principal velocities" 

c c c 



/fiea 



/flEy 



'fl€z 



(A.2) 
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Note (this is a tricky point that has the potential to cause confusion) that v x is not 
the velocity of light in the x direction — since e x (and so v x ) is related to the prop- 
erties of the electric field in the x direction, the principal velocity v x is instead the 
velocity of a light wave whose electric field is pointing in the x direction. That is, 
for light waves propagating in the y-z plane, one of the polarizations will propa- 
gate with speed v x . 

A.2 Group velocity and ray equation 

The group velocity, v g , in the framework used by Born and Wolf, is identical to 
the "ray velocity", and is controlled by the so-called "ray equation". See (15.2.29), 
page 797. To set some conventions, h will always denote a unit vector in physical 
space — a unit with respect to the usual Euclidean norm, while n is a generic 
position in physical 3-space. In contrast, k will be reserved for a unit wave-vector 
in the dual "wave-vector space". 

Born and Wolf exhibit the ray equation in a form equivalent (Born and Wolf 
use t where we use n) to: 

"' "' j "' -0. (A.3) 



W g - Vvl 1/v 2 - 1/v 2 1/v 2 - 1/v 2 
Here the group velocity (ray velocity) is defined by looking at the energy flux and 



v g = v g h. (A.4) 



We can rewrite this as 



_^M_ + m. + ^M_ = o. (A.5) 

v — v v — v v — v 

9 x 9 y 9 z 

This form of the ray equation encounters awkward "division by zero" problems 
when one looks along the principal axes, so it is advisable to eliminate the denom- 
inators by multiplying through by the common factor (v 2 — v1)(v 2 — v 2 )(v 2 — v 2 ), 
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thereby obtaining: 

~2 2/ 2 2\/ 2 2\ , - 2 2 / 2 2\ / 2 2\ 

"aWp - V y)( V 9 - V z) + n y V y{ V g ~ V z)^g ~ %) 

+h 2 z v 2 (v 2 -v 2 x )(v 2 -v 2 y )=0. (A.6) 

It is this form of the ray equation that, (because it is much better behaved), we shall 
use as our starting point. Now this is clearly a quartic in v g , and by regrouping it 
we can write 

4|~"2 2 , -2 2 , -2 21 
^ L n * U * + Vl/ + n x V z\ 

~v 2 9 [n 2 x v 2 x (v 2 y + v 2 z ) + h 2 y v 2 (v 2 z + v 2 x ) + h 2 z v 2 z (v 2 x + v 2 y )] 

+ [v 2 x v 2 y v 2 ] = 0. (A.7) 

Equivalently 

4 |~"2 -2 -2 , -2 -2 -2 , -2 -2 -21 
V s K u „ U 2 + U y V z V x + n z V x V y \ 

~ v l [™l( v y 2 + v z 2 ) + K^f + v x*) + n&x 2 + v y 2 )] 

+1 = 0. (A.8) 



Now define two quadratics (in terms of the three direction cosines hi) 

- /- s.\ |~"2 -2 -2 i -2 -2 -2 i -2 -2 -21 /a n\ 

<?o(n,n)= [n,^ w 2 + n y v z v x + n z v x v y J, (A.9) 

q 2 (h,h) = -[hl{ v - 2 + v- 2 )+h 2 y {vf + vf)+h 2 z {v x 2 + v y 2 )], (A.10) 
then 



2/ A \_ ga(n, n) ± yg 2 (n, h) 2 - g (h, h) 
9o(n,n) 
But, since h is a unit vector, we could equally well rewrite this as 



2,~x _ g2(n, n) ± Vg 2 (n, h) 2 - g (n, h) (h ■ h) 

9o(n,n) 
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In this form both numerator and denominator are manifestly homogeneous and 
quadratic in the components of n, so for any 3-vector n (now not necessarily of 
unit norm) we can take the further step of writing 



q 2 (n, n) ± ^(n, n) 2 - q (n, n) (n ■ n 



<(n) = ^ 7 ' V ™V ' ™ v ' /v '-■ (A.13) 



g (n,n) 



The function v g (n) so defined is homogeneous of degree zero in the components 
of n: 

Vg(Kn) = v g (n) = v g (h). (A.14) 

The homogeneous degree zero property should remind one of the relevant feature 
exhibited by the Finsler metric. It is also useful to note that 



q 2 (n, n) =F Vg 2 (n, n) 2 - q (n, n) (n • n) 



v g (n) 2 (n ■ n) 



(A.15) 



A.3 Phase velocity and Fresnel equation 

In contrast, the phase velocity, in the framework used by Born and Wolf, is con- 
trolled by the so-called "equation of wave normals", also known as the "Fresnel 
equation". See equation (15.2.24), page 796. The relevant computations are similar 
to, but not quite identical to, those for the group velocity. 

Let us consider a plane wave exp(i[k • x — cot]) and define the phase velocity by 

v p = v p k = — k, (A.16) 

then the Fresnel equation is equivalent (Born and Wolf use s where we use k) to 

P P P 

+ -A^ + -A. = 0. (A.17) 



V 2 — f 2 f 2 — f 2 f 2 — f 2 
u %> u x u p u y u p u z 



This form of the equation exhibits "division by zero" issues if you try to look 
along the principal axes, so it is for many purposes better to multiply through by 
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the common factor (v 2 — v 2 )(v 2 — v 2 )(v 2 — v 2 ) thereby obtaining the equivalent of 
their equation (15.3.1) on page 806: 



^i-v^vi-^)^{vi-vi){vi-vi) 



+ k 2 z (v 2 p - v 2 x )(v 2 - v 2 ) = 0. 



(A.18) 



This is clearly a quartic in v v and by regrouping it, and using k ■ k = 1, we can 
write 



4 2 

V p- V p 



K(vl + vl) + k>(vl + vl) + h>M + iZ) 



+ 



7222,7222,7222 

K v v v z + K V z V x + k z v x v y 



0. 



(A.19) 



Let us now define two quadratics (in terms of the direction cosines hi) 

1 



q 2 (k, k) = ^ k 2 x (v 2 y + v 2 z ) + k 2 y (v 2 + v 2 x ) + k 2 z (v 2 x + v 2 y ) 



and 



7222, 7222, 7222 
K x V y V z + K y V z V x "+" K z V x V y 



?o(k,k) 
so as a function of direction the phase velocity is 



«i(k) = fe (k, k) ± Vg 2 (k, k)^ - g (k, k). 



(A.20) 
(A.21) 

(A.22) 



This is very similar to the equations obtained for the ray velocity. In fact, we can 
naturally extend this formula to arbitrary wave-vector k by writing 



_ g 2 (k,k)± v /g 2 (k,k)2-g (k,k)(k-k) 

' Jp{k) ~ (]Tk) 



(A.23) 



This expression is now homogeneous of order zero in k, so that 

v p (Kk) = v p (k) = v p (k). 



(A.24) 



Again, we begin to see a hint of Finsler structure emerging. 
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A.4 Connecting the ray and the wave vectors 

Connecting the ray-vector n and the wave-vector k in birefringent optics is rather 
tricky — for instance, Born and Wolf provide a rather turgid discussion on page 
798 — see section 15.2.2, equations (34)-(39). The key result is 



v g (n) hi v p {k) ki 



2 — /j)2 „, ( lf\2 — /i.2 



v g [nr — vj v n (kr — v. 



v p (ky 



(A.25) 



which ultimately can be manipulated to calculate h as a rather complicated "ex- 
plicit" function of k — albeit an expression that is so complicated that even Born 
and Wolf do not explicitly write it down. Unfortunately if it comes to pseudo- 
Finsler geometry, all the extra technical machinery provided by Finsler notions of 
norm and distance do not serve to simplify the situation. (The fact that phase and 
group velocities can be used to define quite distinct, and in some situations com- 
pletely unrelated, effective metrics has also been noted in the context of acoustics 
B771H75I.) 

A.5 Optical axes 

To find the ray optical axes we (without loss of generality) take v z > v y > v z , and 
define quantities A± (this is of course the result of considerable hindsight) by: 

"* = t4^ ; ^irfc' (A ' 26) 

Furthermore eliminate h y by using 

\2 i -2 -2 



hi = l-n 2 x -hi, (A.27) 



then 

D = Ql - Qo 
1 

" 4 



1 [(^A + + ^A_) 2 -(A^ + A 2 _)] 



x 



[(h x A + - h z A_) 2 - (A 2 , + A 2 )] . (A.28) 
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Thus the (ray) optical axes are defined by 

(n x A + ± n z A_) 2 = (AJ. + A 2 ). (A.29) 

But thanks to the Cauchy-Schwartz inequality 

(n x A + ± n z A_) 2 < (n 2 + n 2 z )(A 2 + + A 2 _) < (A 2 + + A 2 _). (A.30) 

Therefore on the (ray) optical axis we must have n y = 0, and (n 2 + n 2 ) = 1. So (up 
to irrelevant overall signs) 



A_ 



ei, 2 







A_ 



A^ + A 2 



'A% + A 2 
which we can recast in terms of the principal velocities as 



(A.31) 



e l,2 




i/«2-i/«2 : 



o 



' iM 2 -i/^ 

1/t; 2 - 1/t; 2 



or 



ei,2 



±- 



W 2 hl-V 2 y _ _V^ \ V 2 y ~ Vj 

r (/ y v 2 x -vf 'vyVvl-vl 



(A.32) 



(A.33) 



These are the two ray optical axes. (Compare with equation (15.3.21) on p. 811 of 
Born and Wolf.) 

A similar computation can be carried through for the phase optical axes. We 
again take v z > v y > v z , and now define 



2 2 , a 2 

v x = v y + A + ; 



2 2 a 2 

v, = v,, — A_. 



Eliminate k y by using 

Then 

D 



l-k 2 -ki 



Q2 ~Qi 
1 



4 
x 



(A; :E A + + A; z A_) 2 -(A 2 h + A 2 _; 
\k x A + -k z A_) 2 -(A 2 + + A 2 _ 



(A.34) 

(A.35) 

(A.36) 

(A.37) 
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This tells us that the discriminant factorizes, always. The discriminant vanishes if 

(^A + ±£ 2 A_) 2 = A^ + AL (A.38) 

But by the Cauchy-Schwartz inequality 

(fc x A+ ± fc 2 A_) 2 < (P x + k 2 z )(A 2 + + A 2 _) < (A 2 , + A 2 _). (A.39) 

Thus on the (phase) optical axis we must have k y = and (k 2 + k 2 ) = 1. The two 
unique directions (up to irrelevant overall sign flips) that make the discriminant 
vanish are thus 

e 1>2 = |± f + ; ; f" ) . (A.40) 

y v/A^ + A 2 . ^A 2 . + A 2 

which can be rewritten as 



e, 2 =lW|^i;0;^|. (A.41) 

These are the two phase optical axes. (Compare with equation (15.3.11) on p. 810 
of Born and Wolf.) 



Appendix B 

Special functions: some important 
formulas 



B.l Trigonometric identities 

In the body of the thesis we needed to use some slightly unusual trigonometric 
identities. They can be derived from standard ones without too much difficulty 
but are sufficiently unusual to be worth mentioning explicitly: 

„ cos(A - B) - cos(A + B) 

tan A tan B = K — K — '-; (B.l) 

cos(A -B) + cos(A + BY 

(A + B\ [A-B\ cos B- cos A m x 

tan tan = -; (B.2) 

V 2 J \ 2 J cos^ + cosA' v ; 

and 

cos{A + 2B) + cos A = 2 cos B cos{A + B). (B.3) 



B.2 Gamma function identities and approximations 

The key Gamma function identity we need is 

Y{z) r(l -z) = -^—. (B.4) 
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The Stirling approximation for Gamma function gives 



T(x) 




2 ^ x -Y(i + o( l - 

x \eJ V V x 



Re(x) > 0, |x| — y oo. 



(B.5) 



We also need the following asymptotic estimate based on the Stirling approxima- 
tion 



T(x + 
T(x) 



x 



l + 0[- 

x 



Re(x) > 0, 



\x\ — > oo. 



(B.6) 



B.3 Hypergeometric function identities 



The key hypergeometric function identities we need are Bailey's theorem 

r 



r(f)r(*i; 



2 Fl ( a, 1 ~a,C, - , T^/c+gNp/ c-q+l y 



(B.7) 



which is easily found in many standard references (for example H166I0 , and the 
particular differential identity 

d{ 2 Fi(a,b,c,z)} c-1 



[ 2 Fi (a, b,c-l,z)- 2-Fi («, b, c, z)} 



dz z 

which is easy enough to verify once it has been presented. 



(B.8) 



B.4 Bessel function identities and expansions 

The important differential identity which holds for Bessel functions is the follow- 
ing: 



dJ n (x) 



[j Q _i(x) - J a+1 (x)]. 



dx 2 

What was also needed was the asymptotic expansion (see for example [166]) 



(B.9) 



Ja{x) = \l 




2 

TIX 



r (a, x) cos I x 

2 4 



Q(a,x)sm [ x - — - - 



(B.10) 
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where 

p(a, x) = jt(-iy - ^\t a ^T ] 9 ^ (B - n) 

^ (2x) 2n (2n)\r(± + a-2n) 

and 

«<**> ^ B-'i'i^^M' <B ' 12) 
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D(R n ) 


compactly supported, smooth functions 

(C°°(R n )) 
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D'(R n ) 


space of distributions on the space of com- 
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102 
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£"•/ I e ) 
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M^ n ) 


Moderate functions 


105,107,111 


N(R n ) 


Negligible functions 


105,107,111 


mollifier, 
smoothing kernel 


(usually) e— sequence of smooth, compactly 
supported functions with integral normed to 
1 and with support "stretching" to as e — > 
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C(f) 


embedding of distributions into Colombeau 
algebra by a convolution with a mollifier 


107 


M 


(in the section|3.4[) manifold, on which one can 
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A 


maximal piecewise smooth atlas, where the 
transformation Jacobians are bounded on ev- 
ery compact set 
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maximal smooth subatlas of A 
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^Ch 


subset of a manifold, such that it can be 
mapped to R n by a chart mapping 
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Ch(n Ch ) 


chart from the atlas A on the set VL C h 
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c p (n Ch ) 


class of compactly supported piecewise 
smooth 4-forms, (we work from the be- 
ginning with a 4D manifold), having their 
support inside the set Vt C h 


124 


s 


maximal subatlas of A, such that there exist el- 
ements of the class C p (il C h), that have in this 
subatlas smooth scalar densities 


124 


C S(S)^ Q Ch) 


set of all 4-forms from C p (Vtch), that are in S 
given by smooth scalar densities 
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C s (£lch) 


^c p (S) (n Ch ) 
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C P (M) 
C£(M) 


Un Ch C p (Q Ch ) 
^n Ch C s (D, C h) 


124 


D'(M) 


linear generalized scalar fields 
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D [s)W 


linear generalized scalar fields defined (at 
least) on the class C P ,JM) 
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D \a M ^ 


linear generalized scalar fields defined exclu- 
sively on the class C P ^JM) 
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D' E {M), 
D ' E( . n s n0) ( M ) 


distributional analogue of piecewise smooth 
scalar fields 
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D's(M), 

S(U n S n o)\ > 


objects as close as possible to a distributional 
analogue of smooth scalar fields 
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D'a(M) 


generalized scalar fields 
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D' EA {M), 
D' SA {M) (etc.) 


generalized scalar fields constructed from the 
elements of the classes D' E (M), D' S (M) (etc.) 


129 


DT(M) 


linear generalized tensor fields of rank (m, n) 


127-128 


At(S) 


specific function, which maps atlases to at- 
lases 


129 



183 



Symbol: 


Definition/Explanation: 


Page: 


D' m - (M) 


linear generalized tensor fields of rank (m, n), 
defined (at least) on all the classes C^.AM), 
labeled by I, in the atlases At (Si) 
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D' m - (M) 


linear generalized tensor fields of rank (m,n), 
defined exclusively on the classes Cg, s JM), 
labeled by I, in the atlases At (Si) 


128 


^sw 


distributional analogue of rank (m, n) piece- 
wise smooth tensor fields 


128 


^3W 


objects as close as possible to a distributional 
analogue of rank (m, n) smooth tensor fields 


128 




subclass of D'7 A , S „(M), such that At (Si) 
has as image of Si the atlas Si C Si 


129 


L>^(M) 


generalized tensor fields of rank (m, n) 
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D- A (M), D£ A (M ) (etc.) 


generalized tensor fields of rank (m, n) con- 
structed from the elements of the classes 

DTe(M), DTs(M) 
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r m (M), rf (m), rj?(M), 
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T m - (M), 


Gamma objects (generalized from linear gen- 
eralized tensor fields of rank (a, b), by impos- 
ing no condition on the transformation prop- 
erties) 
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rj(M) 


sets of algebras constructed from the elements 
of the class T m (M) (generalization of D™ A (M)) 
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r^(M),r- A (M)(etc.) 


sets of algebras constructed from the elements 
of the classes T™(M), V r §(M) (etc.) (general- 
ization from D'™ A (M), D'™ SA (M) ) 
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r£:-(Cfc fc ) 


multi-index matrix obtained from the tensor 
components of Tg;;; in the chart Ch k 
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C/i'(g,fi C fc) 


chart mapping from Qch to the whole R 4 , such 
that it maps the point q to 
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A n (5,C/i'(9,«Ch)) 


specific subclass of the class C^xJttch) 


133 



184 



Symbol: 


Definition/Explanation: 


Page: 


u e (y) 


specific one parameter class of elements 
of the class Cf. c -.(M) 


133 


ftj 


the equivalence relation 


133-134 


f m (M), fg(M), 
DT A {M),DT EA {M) 
(etc.) 


sets of equivalence classes constructed 
from the elements of the classes T m (M), 
F%(M), DTa(M), D'™ a (M) (etc.) 


134 


A 


specific subclass of the class 

r m - (M) 

1 E(u n At(S n )o)^ VI 1 


135-136 


^&) 


tensor field associated to T£;- G PJ(M) 


135 


6(C7i fc ) 


set of Lebesgue measure 0, on which a 
particular multi-index matrix T£\"(Chk) 
is discontinous 


136 


(^-derivative 


specific operator acting on the elements 
of the class r™ (M) and generalizing the 
operator of the partial derivative (as a 
part of the covariant derivative defini- 
tion) 
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Dc(u) 


covariant derivative along the vector 
field U* 
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v n 


n-times differentiable subatlas of A 
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S' n related to the 
atlas V n 


specific subclass of r™ (u ^ t( ^ )o)A (M) 
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Jozef Skakala and Matt Visser, The causal structure of spacetime 
is a parametrized Randers geometry 

Abstract: 

There is a by now well-established isomorphism between stationary 4-dimensional 
spacetimes and 3-dimensional purely spatial Randers geometries - these Randers 
geometries being a particular case of the more general class of 3-dimensional 
Finsler geometries. We point out that in stably causal spacetimes, by using the 
(time-dependent) ADM decomposition, this result can be extended to general 
non-stationary spacetimes - the causal structure (conformal structure) of the full 
spacetime is completely encoded in a parameterized (time-dependent) class of 
Randers spaces, which can then be used to define a Fermat principle, and also to 
reconstruct the null cones and causal structure. 
published in Classical and Quantum Gravity 28, (2011), 065007, 
[arXiv: gr-qc/1012.4467]. 
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Jozef Skakala and Matt Visser, Birkhoff-like theorem for rotating 
stars in (2+1) dimensions 

Abstract: 

Consider a rotating and possibly pulsating "star" in (2+1) dimensions. If the star 
is axially symmetric, then in the vacuum region surrounding the star, (a region 
that we assume at most contains a cosmological constant), the Einstein equations 
imply that under physically plausible conditions the geometry is in fact stationary 
Furthermore, the geometry external to the star is then uniquely guaranteed to be 
the (2+1) dimensional analogue of the Kerr-de Sitter spacetime, the BTZ geometry 
This Birkhoff-like theorem is very special to (2+1) dimensions, and fails in (3+1) 
dimensions. Effectively, this is a "no hair" theorem for (2+1) dimensional axially 
symmetric stars: the exterior geometry is completely specified by the mass, angu- 
lar momentum, and cosmological constant. 
arXiv: gr-qc/0903.2128. 



207 



